Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,966,386 A→C intergenic (+296/‑31) trxA → / → rho thioredoxin 1/transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,966,3860AC93.8% 209.5 / 5.1 66intergenic (+296/‑31)trxA/rhothioredoxin 1/transcription termination factor
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base C (28/33);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

AACAGATTTCTGGCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTG  >  NC_000913_3_pae_pgi/3966241‑3966525
                                                                                                                                                 |                                                                                                                                           
aaCAGATTTCTGGCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTa                                                                                                                                         >  2:44602/1‑150 (MQ=255)
aaCAGATTTCTGGCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTa                                                                                                                                         >  1:773886/1‑150 (MQ=255)
aaCAGATTTCTGGCTCGTCACTCAATCCATCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGCGCTTTa                                                                                                                                         >  2:394710/1‑150 (MQ=255)
aaCAGATTTCTGGCTCGTCACTCAATCAGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTa                                                                                                                                         >  1:640040/1‑150 (MQ=255)
     aTTTCTGGCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCGGCGTACTCGCCCGTGACCAGGAACCGAAAAGACATAAGCCGCTAACCATAAACAGGCCTGGATGATCCTGCCGTGCCATTCACAACATTAAATTAGAGGTTTGccccg                                                                                                                                    >  4:543749/1‑149 (MQ=255)
       ttCTGGCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAg                                                                                                                                  >  1:188003/1‑150 (MQ=255)
        tCTGGCTCGTCCCTCAATCCGTCTTGTCGTTTCAGGTCGTCGTACTCTCCTGTGACCAGGCCGAGAAAAGACATGAGTCGATGACCGTAAACAGGCAGGGCTGATCCTGCCATACCATTCACAACATTAATTTCGAGCGTTCCCCCTAGt                                                                                                                                 >  2:687435/1‑150 (MQ=255)
          tGGCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGttt                                                                                                                               <  2:548958/150‑1 (MQ=255)
           ggcgaggcACTGAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATTACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTa                                                                                                                              <  2:477134/143‑1 (MQ=255)
            gCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTaa                                                                                                                             <  4:517851/150‑1 (MQ=255)
              tCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGa                                                                                                                           >  3:120440/1‑150 (MQ=255)
                 tCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACt                                                                                                                        <  2:793286/150‑1 (MQ=255)
                 tCACTCAATCCGTCTGGTCGTTTCAGTTCTGCGTACTCTCCTGTGCCCGGGCAGCGAAAAGACATGAGTCGATGACCGTAAACATGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACt                                                                                                                        <  2:128797/150‑1 (MQ=255)
                   aCTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTca                                                                                                                      >  2:104550/1‑150 (MQ=255)
                      cAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAGCGTTAAGAACTCacac                                                                                                                   >  1:720195/1‑150 (MQ=255)
                            gTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTAt                                                                                                             <  4:317423/150‑1 (MQ=255)
                              cTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGa                                                                                                           <  2:482242/150‑1 (MQ=255)
                                tGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAAt                                                                                                         <  2:270366/150‑1 (MQ=255)
                                 gTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACTCCACTATGAATc                                                                                                        <  2:379191/150‑1 (MQ=255)
                                 gTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATc                                                                                                        <  3:387564/150‑1 (MQ=255)
                                  tCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCt                                                                                                       >  4:310255/1‑150 (MQ=255)
                                      ttCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTAcc                                                                                                   <  3:179794/150‑1 (MQ=255)
                                              tGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTaaa                                                                                           <  2:366138/150‑1 (MQ=255)
                                              tGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCGTACCGAATTaaa                                                                                           <  1:856731/150‑1 (MQ=255)
                                                       tcCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGcc                                                                                  >  1:627722/1‑150 (MQ=255)
                                                        cgtgtgACCAGGCAGCGAAAATACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCg                                                                                 <  3:278072/148‑1 (MQ=255)
                                                           gtgACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGtt                                                                              >  4:344209/1‑150 (MQ=255)
                                                            tgACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGttt                                                                             >  2:180263/1‑150 (MQ=255)
                                                              aCCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCt                                                                           <  2:132762/150‑1 (MQ=255)
                                                                cAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGATTTAAGGAATACGCCGGTTTCTGa                                                                         <  2:256085/150‑1 (MQ=255)
                                                                cAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGa                                                                         <  4:433757/150‑1 (MQ=255)
                                                                    cAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTg                                                                     <  2:697686/150‑1 (MQ=255)
                                                                      gCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAt                                                                   <  2:373791/150‑1 (MQ=255)
                                                                      gCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAt                                                                   <  4:257727/150‑1 (MQ=255)
                                                                        gAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCa                                                                 >  2:831436/1‑150 (MQ=255)
                                                                          aaaGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAct                                                               >  3:84605/1‑150 (MQ=255)
                                                                           aaGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTTATCActc                                                              <  2:891269/150‑1 (MQ=255)
                                                                               cATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCgg                                                          <  2:640040/150‑1 (MQ=255)
                                                                                  gAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGa                                                       >  4:374475/1‑150 (MQ=255)
                                                                                   aGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaa                                                      >  2:733640/1‑150 (MQ=255)
                                                                                   aGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAAt                                                                                                         <  3:562037/99‑1 (MQ=255)
                                                                                   aGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAAt                                                                                                         >  4:562037/1‑99 (MQ=255)
                                                                                    gccGATGACCGTAAACAGGCATGGATAATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAACCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaaa                                                     <  1:633384/148‑1 (MQ=255)
                                                                                    gTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaaa                                                     >  2:241706/1‑150 (MQ=255)
                                                                                        aTGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAatat                                                 >  2:821566/1‑150 (MQ=255)
                                                                                        aTGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAatat                                                 >  1:369439/1‑150 (MQ=255)
                                                                                             cGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGc                                            >  2:49491/1‑150 (MQ=255)
                                                                                              gTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCt                                           <  1:287427/150‑1 (MQ=255)
                                                                                                    aGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATgg                                               >  1:249183/1‑140 (MQ=255)
                                                                                                    aGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATgg                                               <  2:249183/140‑1 (MQ=255)
                                                                                                    aGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGaaaa                                     >  2:674353/1‑150 (MQ=255)
                                                                                                     ggCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAAc                                    <  1:49491/150‑1 (MQ=255)
                                                                                                      gCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                   <  2:391207/150‑1 (MQ=255)
                                                                                                        aTGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTg                                 >  4:135388/1‑150 (MQ=255)
                                                                                                          ggATGATCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGc                               >  4:496574/1‑150 (MQ=255)
                                                                                                            atgatCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                             <  2:773886/150‑1 (MQ=255)
                                                                                                              gatCCTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                           <  4:130070/150‑1 (MQ=255)
                                                                                                               atCCTGCCATACCATTCACAACAGTAAGATCGAGCGTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATAAAAAAAAAAAAAAATATGGGGCTGGAAAACCTGGCTCGTa                          <  2:369439/150‑1 (MQ=255)
                                                                                                                 ccTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATg                        <  1:436118/150‑1 (MQ=255)
                                                                                                                  cTGCCATACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc                       >  1:202274/1‑150 (MQ=255)
                                                                                                                       aTACCATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                  <  3:344209/150‑1 (MQ=255)
                                                                                                                           cATTCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg              >  2:298910/1‑150 (MQ=255)
                                                                                                                             ttCACAACATTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc            >  2:735214/1‑150 (MQ=255)
                                                                                                                                acaacaTTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatt         >  2:873596/1‑150 (MQ=255)
                                                                                                                                     aTTAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt    <  1:733640/150‑1 (MQ=255)
                                                                                                                                       tAAGTTCGAGCTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTg  <  4:409642/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
AACAGATTTCTGGCTCGTCACTCAATCCGTCTTGTCGTTTCAGTTCTGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTG  >  NC_000913_3_pae_pgi/3966241‑3966525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: