Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,7210GT45.6% 108.1 / 43.9 79Q72K (CAA→AAA) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base G (11/32);  new base T (36/0);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 5.73e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATG  >  NC_000913/300579‑300863
                                                                                                                                              |                                                                                                                                                     
gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTTTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCCAACGCCCCGCGCCTTTTTTGCGGCGt                                                                                                                                               >  2:298314/1‑151 (MQ=255)
     ggCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCCTTTTTCGCGGCGGcccga                                                                                                                                          >  2:922202/1‑148 (MQ=255)
      gCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCACTACCCCCCGCGCGCTTTTTGCGGCGGCCGCAt                                                                                                                                       >  1:983301/1‑151 (MQ=255)
       ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATc                                                                                                                                     <  1:1086875/151‑1 (MQ=255)
        cACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCa                                                                                                                                   <  2:60748/151‑1 (MQ=255)
          cTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAACCCCTACCCCCCCCCCTTTTTTTGCGGCGGGCGTAttcaa                                                                                                                                   >  1:335596/1‑147 (MQ=255)
          cTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAACCCGCACGCCCCGCGCCTTTTTCGCCGGGGCCGTATCCAg                                                                                                                                 >  2:350145/1‑151 (MQ=255)
           tACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAgg                                                                                                                                <  2:1055375/151‑1 (MQ=255)
           tACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAgg                                                                                                                                <  2:1120188/151‑1 (MQ=255)
            aCCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAggg                                                                                                                               <  1:478363/151‑1 (MQ=255)
            aCCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCAGGAAGCCCCGCGCCTTTTTCGCGGCGGCCCGATTAAggg                                                                                                                               >  2:1233161/1‑151 (MQ=255)
                  gcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACGCCCCGCGCTTTTTTTGCGGGGGCCGGATTCAAgggggggg                                                                                                                          >  2:774636/1‑149 (MQ=255)
                   cAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGa                                                                                                                         <  1:187594/150‑1 (MQ=255)
                       ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCCTTTTTCGCGGGGTCCGTATCCAAGGCGGGGGGGCgg                                                                                                                   >  1:338284/1‑150 (MQ=255)
                       ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGATACCCAGCAACCCCCGCGCCTTTTTCGCGGGGCCCGCGTCTAAGGCGGGGGGGCgg                                                                                                                   >  1:807305/1‑150 (MQ=255)
                        gCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACCCCCCGCCCCTTTTTCGCGGCGGCCCGATCAAAGgggggggggggg                                                                                                                    >  2:90406/1‑142 (MQ=255)
                           tgatgaTAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAgg                                                                                                                      <  1:1035801/144‑1 (MQ=255)
                           tgatgaTAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACCGCCAGTACCCCCCGCGCCTTTTTCGCGGCGGCCGTATCCAggggggggggggtgtgt                                                                                                                 >  2:1036309/1‑139 (MQ=255)
                            gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGGGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTGAATATAGCCAATACGCCCCGCGCTCTTTGCGCGGCGTCCGGGTAAAGGGCGGGGAgggggcccgt                                                                                                              >  1:2464/1‑143 (MQ=255)
                              tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCCTTTTTCGCGGGGGCCGTTTCCAGGGGGGGGGGGGGTCCttttg                                                                                                          >  2:847495/1‑149 (MQ=255)
                              tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAACCAGTACGCCCCGCGCCTTTTTTGCGGCGGCCGGATCAAAGGGGGGGGGGGGGCCCTTgg                                                                                                          >  2:504584/1‑151 (MQ=255)
                              tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGCGTTTTTGTCGCCTTTGCCGAGCGCCCCCGCGTTACTGGCGGTCATTACAGCCCCTATGCCCCGTCCTTTTTTCGCGGCGCCCGGCTCCAGGGTGGGGCAGCGTCCctttg                                                                                                          >  1:1313212/1‑149 (MQ=255)
                                atAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGAAAACCCGGCCGCCCCGCGCCTTTTTTGCGGGGGCCGGATCCAGGGGGGGGGGGGGTCCCTTTgga                                                                                                         >  2:643760/1‑150 (MQ=255)
                                   gTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACCCCCCGCGCTTTTTTGGCGGGGGCCCGATTAAGGGGGGGGGGGGGGCCTTTTGGCaag                                                                                                      >  2:920402/1‑150 (MQ=255)
                                   gTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGAAGAACGCCCGTACACCCCGCGCCTTTTTTCCGGGGGCCGTATCCAAGGGGGGGGGGGGTCCCTTTgggatg                                                                                                      >  1:1406619/1‑147 (MQ=255)
                                   gTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAGCCCGGACGCCCCGCGCCTTTTTTGCGGGGGCCGCGTCCAAGGGGGGGGGGGGGCCCTTTGGCatg                                                                                                      >  2:923595/1‑149 (MQ=255)
                                   gTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGAAAACCCATACGCCCCGCGCTTTTTTGGCGGGGGCCCTATCAAAGGGGGGGGGGGGGCCCTTTGCCaag                                                                                                      >  2:1078857/1‑150 (MQ=255)
                                      cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                                                  <  2:119702/151‑1 (MQ=255)
                                        cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAggg                                                                                                                               <  2:249024/123‑1 (MQ=255)
                                        cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATAACAACCCGTACCCCCCGCGCCTTTTTCGCCGCGGCCGGATCAAGGGGGGGGGGGGGGCCCTTTTCCATGGCggg                                                                                                >  1:485846/1‑150 (MQ=255)
                                          aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACCCCCCGCGCTTTTTTCGCGGCGCCCGGCTCAAAGGGGGGGGGGGGGCCttttggcatgggggga                                                                                               >  2:599474/1‑137 (MQ=255)
                                          aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAc                                                                                              <  1:274112/151‑1 (MQ=255)
                                          aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTACAGGCAGTACGCCCCGCGCTTTTTTTGCGGGGGCCGCATTAAAGGGGGGGGGGGGGCCCTTTGGGAAGGCggga                                                                                               >  2:444932/1‑148 (MQ=255)
                                          aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTACAAGCAGTACCCCCCGCGCCTTTTTTGCGGGGGCCGGATCAAGGGGGGGGGGGCGGCCttttgggaagggggga                                                                                               >  2:668648/1‑137 (MQ=255)
                                          aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGTAAACCAGTACGCCCCGCGCCTTTTTGGCGGCGCCCCGATCCACGGGGGGGGGGCGTGCCTTTGGAAAGGGGGGaa                                                                                               >  1:1028186/1‑151 (MQ=255)
                                           tAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTACAGCCCGTACCCCCCGCGCCTTTTTCGCGGCGGCCGTCTCCAGGGCGGGGGGGCGGCCTTTTGccaatgcgggaa                                                                                              >  1:1076836/1‑139 (MQ=255)
                                                  gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGat                                                                                       <  2:1090864/150‑1 (MQ=255)
                                                   ccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat                                                                                     <  2:830841/151‑1 (MQ=255)
                                                   ccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAACCAGGACGCCCCGCGCTTTTTTCGCGGGGTCCCGATCCAAGGGGGGGGGGGGGCCCTTTGCCATGGCGGGACCCACGatat                                                                                     >  2:734723/1‑151 (MQ=255)
                                                     gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTACAGCCCGTAAGCCCGGCGCTTTTTTCGCGGCGGCCGTCTCCAGGGGGGGGGGGGGGCCTTTTGGCAAGGGGGGAACGGCGATATAg                                                                                   >  2:1188272/1‑151 (MQ=255)
                                                     gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAACCCGTACCCCCCGGGCTTTTTTCGCGGGGGCCGGATCAAAGGGGGGGGGGGGTCCCTTTGCCATGgggggaaccccgaaaaaa                                                                                   >  2:662892/1‑137 (MQ=255)
                                                     gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGTATGAAAGCCCGGACGCCCCGCGCCCTTTTTGCGGCGGCCGGTTCCAAGGGGGGGTGGGGGCCCTTTGGCATGGGCGGGAccgccggataa                                                                                    >  1:1268968/1‑146 (MQ=255)
                                                     gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTTCAACCACGACGCCCCGCGCTTTTTTTGCGGGGTCCCGCTCAAAGGCGGGGAGGGGTCCCTTGGCGAAGGCGGGAACCACGGtataa                                                                                   >  1:1315744/1‑150 (MQ=255)
                                                     gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGTCAACCACTACGCCCCGCGCTTTTTTCGCGGCGGCCGCCTCAAGGGCGGGGCGGCGGCCCTTGGCATTGTCGGTAGCGACGCTACAg                                                                                   >  1:1127741/1‑151 (MQ=255)
                                                      ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGc                                                                                  <  1:297006/151‑1 (MQ=255)
                                                        tAACAGCCTGGCGGTGTTTTTGTCGCTTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGAGAAACAAGACCCCCCGCGCTCTTTTCCCGGCGCCCCAAGTAAGGGGGGTGAGGGCTCCCTTTGGCACGGGGGAACCGACGATATAAcc                                                                                 >  2:982742/1‑150 (MQ=255)
                                                             gCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCGGCGCTTTTTTCGCGGCGTCCGTATCAAGGGGGGGGGGGCGGCCTTTGGGCAAGGCGGGACCGGCGATATAGCCCTaagg                                                                           >  2:39204/1‑148 (MQ=255)
                                                                   cGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGt                                                                      <  1:64276/150‑1 (MQ=255)
                                                                        tttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgggg                                                                 <  2:1006400/150‑1 (MQ=255)
                                                                            tGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttc                                                            <  2:1081655/151‑1 (MQ=255)
                                                                              tCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttctt                                                          <  1:1443498/151‑1 (MQ=255)
                                                                               cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCttcttc                                                         <  1:871128/151‑1 (MQ=255)
                                                                                 ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGTTATGGCAGTCACTACCCCCCGCCCCTTTTTCGCGGGCTCCGCATCAAAGGCGGTGATGCGCCCTTTTGGAATTGCGGAACCCACGAAATAGCCAGAAGCGGCGGTGGGGGctgtttt                                                         >  1:187808/1‑144 (MQ=255)
                                                                                 ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATAATAGCCAGTACGCCCCGCGCTTTTTGCGCGGGGGCCGCACCAAGGGGGGGGGGGCGTCCCTTGGCAATGGCGGGAACGCCCGAATATCCCTAGGCGGCCGCGGGGGGgtcttc                                                         >  2:26279/1‑149 (MQ=255)
                                                                                 ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCAGGACGCCCGGCGCTTTTTGCGCGGCGTCCGTATTAAGGGGGGGGGGGGGGCCTTTGGCAATGGGGGGAACGACGATATAGCCATAGGGGGGGGTGGGGGggtcttt                                                         >  2:699595/1‑148 (MQ=255)
                                                                                   tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGGCCGGATCAAGGGGGGGGGGGGGTCCTTTGGCAATGGGGGGACCGACGATATAGCCATAGGGGGCGGTGGGGGGTTCTTCTTg                                                      >  2:1487488/1‑150 (MQ=255)
                                                                                            tGCCCCGGCGTTGCTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGt                                              <  2:1177771/149‑1 (MQ=255)
                                                                                             gCCCCGGCGTTACTGGCGGTAATTACAACCCCGACGCCCCGCGCCTTTTTCGCGGCGTCCCCATCAAAGGCGGGGAGGGGTCCTTTGGCAATGGCGGGAACCACGATATAGACCTAGGCGGCGGTGGGGGCGTCTTCATGCCCTTCGTAgg                                           >  1:190702/1‑151 (MQ=255)
                                                                                              ccccGGCGTTGCTGGCTGTAATGGCAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGGCGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                           <  1:1481566/150‑1 (MQ=255)
                                                                                              ccccGGCGTTGCTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGGGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                           <  2:907502/150‑1 (MQ=255)
                                                                                              ccccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                           <  1:302656/150‑1 (MQ=255)
                                                                                                    cGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcg                                    <  2:878823/151‑1 (MQ=255)
                                                                                                    cGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcg                                    <  1:594516/151‑1 (MQ=255)
                                                                                                    cGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcg                                    <  1:704991/151‑1 (MQ=255)
                                                                                                      ttACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgt                                   <  1:1465719/150‑1 (MQ=255)
                                                                                                         cTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGGAGcgtgccg                               >  1:1305289/1‑151 (MQ=255)
                                                                                                           ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccgtg                             >  2:609823/1‑151 (MQ=255)
                                                                                                             cGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccgtgc                            <  2:1168288/150‑1 (MQ=255)
                                                                                                                tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCggg                        >  1:410990/1‑150 (MQ=255)
                                                                                                                  aTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGGAg                      >  2:877137/1‑151 (MQ=255)
                                                                                                                       aGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGtttt                  >  2:863241/1‑150 (MQ=255)
                                                                                                                       aGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGtttt                  >  2:755496/1‑150 (MQ=255)
                                                                                                                         ccAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTc                 <  2:106916/149‑1 (MQ=255)
                                                                                                                            gTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCgg            >  2:1437813/1‑151 (MQ=255)
                                                                                                                            gTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCgg            >  2:625896/1‑151 (MQ=255)
                                                                                                                             tACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGc              <  1:1349251/148‑1 (MQ=255)
                                                                                                                             tACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGc              <  1:612237/148‑1 (MQ=255)
                                                                                                                                  ccGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTc      <  2:516786/151‑1 (MQ=255)
                                                                                                                                      cgcTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATg  <  2:1394754/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATG  >  NC_000913/300579‑300863

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