Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,391,723 A→T 5.7% L440* (TTA→TAA)  nuoM ← NADH:quinone oxidoreductase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,391,7230AT 5.7% 282.5 / 3.5 90L440* (TTA→TAA) nuoMNADH:quinone oxidoreductase subunit M
Reads supporting (aligned to +/- strand):  ref base A (41/43);  new base T (2/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCGAGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCGCAAA  >  NC_000913/2391574‑2391871
                                                                                                                                                     |                                                                                                                                                    
gCGAGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTa                                                                                                                                                      >  2:161542/1‑150 (MQ=255)
gCGAGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTa                                                                                                                                                      >  1:1181117/1‑150 (MQ=255)
gCGAGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTa                                                                                                                                                      >  1:1312736/1‑150 (MQ=255)
gCGAGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTa                                                                                                                                                      >  2:219060/1‑150 (MQ=255)
 cGAGGTATACAGAATCGGCTGCGGATAGAAGCCAAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATACCTGGCAGTTACTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTaa                                                                                                                                                     >  2:1508615/1‑150 (MQ=255)
  gAGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAAc                                                                                                                                                    <  1:344443/150‑1 (MQ=255)
   aGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACa                                                                                                                                                   <  2:185350/150‑1 (MQ=255)
    ggTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACAt                                                                                                                                                  >  2:848560/1‑150 (MQ=255)
       aTCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGc                                                                                                                                               >  1:772038/1‑150 (MQ=255)
        tCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGcc                                                                                                                                              <  2:287801/150‑1 (MQ=255)
        tCCAGAATCGGCTGCGGATAGAAGCCCAGCACTACCAGCAGCACCACCAGCAACAGGATCATAAAAAGCTCACGCAGCGACATCCCAGGCAGTTACTGTCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGCTGTAACATCGcc                                                                                                                                              >  2:259053/1‑150 (MQ=255)
          cAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAg                                                                                                                                            >  1:1352021/1‑150 (MQ=255)
          cAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAg                                                                                                                                            >  2:1175634/1‑150 (MQ=255)
           aGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTTACATCGcgctt                                                                                                                                           <  2:1239431/150‑5 (MQ=255)
            gAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCg                                                                                                                                          <  2:778315/150‑1 (MQ=255)
            gAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCg                                                                                                                                          <  1:1573198/150‑1 (MQ=255)
             aaTCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTTACATCGcgcttac                                                                                                                                         <  1:167998/150‑7 (MQ=255)
             aaTCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTTACATCGcgcttac                                                                                                                                         <  2:203233/150‑7 (MQ=255)
              aTCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGaa                                                                                                                                        >  2:244235/1‑150 (MQ=255)
                cGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATa                                                                                                                                      <  2:571663/150‑1 (MQ=255)
                 ggCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATaa                                                                                                                                     >  1:73797/1‑150 (MQ=255)
                 ggCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATaa                                                                                                                                     >  2:252703/1‑150 (MQ=255)
                   cTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAAc                                                                                                                                   >  1:1402355/1‑150 (MQ=255)
                        gATAGAAGCCCAGCAGTACCAGCAGCCCCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATg                                                                                                                              >  1:823428/1‑150 (MQ=255)
                        gATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATg                                                                                                                              >  2:236260/1‑150 (MQ=255)
                        gATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATg                                                                                                                              >  2:239651/1‑150 (MQ=255)
                        gATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATg                                                                                                                              >  2:245616/1‑150 (MQ=255)
                         aTAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGc                                                                                                                             <  2:564653/150‑1 (MQ=255)
                          tAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCa                                                                                                                            >  1:1041896/1‑150 (MQ=255)
                                cccAGCAGTACCAGCAGCACCACCCGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGAcc                                                                                                                      >  2:398934/1‑150 (MQ=255)
                                cccAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGAcc                                                                                                                      <  1:1297759/150‑1 (MQ=255)
                                  cagcagTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAg                                                                                                                    <  1:848560/150‑1 (MQ=255)
                                    gcagTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGcc                                                                                                                  <  2:415358/150‑1 (MQ=255)
                                          ccAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGCCCAGCCCAAAgg                                                                                                            >  2:1494947/1‑150 (MQ=255)
                                          ccAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAgg                                                                                                            <  1:137252/150‑1 (MQ=255)
                                              cagcaCCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTaga                                                                                                        <  2:517537/150‑1 (MQ=255)
                                               agcaCCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTagag                                                                                                       <  2:1352021/150‑1 (MQ=255)
                                                  accaccAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATa                                                                                                    >  1:1512075/1‑150 (MQ=255)
                                                     accaGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACg                                                                                                 >  2:1600118/1‑150 (MQ=255)
                                                     accaGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACg                                                                                                 >  1:757031/1‑150 (MQ=255)
                                                     accaGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACg                                                                                                 >  1:229546/1‑150 (MQ=255)
                                                      ccaGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACgg                                                                                                <  2:282951/150‑1 (MQ=255)
                                                           aaCAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATc                                                                                           >  2:143517/1‑150 (MQ=255)
                                                               ggATCATTACCAGCTCCCGCAGCGACATCGCTGGCAGTTCCTGGCAGGCCATCTGGCTTTTCGCTTTAGCGAAGTACGCGCGATGTAACATCGCCAGCGTAGAAACCGATGCACAGAGCCGCCCCTAGGTTGAGAGAACGGTCATCTCCg                                                                                       >  2:859530/1‑150 (MQ=255)
                                                                gATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCgg                                                                                      <  2:1041896/150‑1 (MQ=255)
                                                                gATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCgg                                                                                      >  1:726818/1‑150 (MQ=255)
                                                                gATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCgg                                                                                      <  1:239651/150‑1 (MQ=255)
                                                                gATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCgg                                                                                      <  1:245616/150‑1 (MQ=255)
                                                                       aaCAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAAcc                                                                               >  2:1299139/1‑150 (MQ=255)
                                                                         cAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTg                                                                             <  2:1312736/150‑1 (MQ=255)
                                                                           gCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGa                                                                           <  2:1185139/150‑1 (MQ=255)
                                                                            cTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGaa                                                                          <  2:431333/150‑1 (MQ=255)
                                                                             tCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAg                                                                         <  2:1510851/150‑1 (MQ=255)
                                                                                  cAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGcc                                                                    <  2:772038/150‑1 (MQ=255)
                                                                                  cAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGcc                                                                    <  1:244235/150‑1 (MQ=255)
                                                                                   aGCGACCTCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCg                                                                   <  2:482182/150‑1 (MQ=255)
                                                                                        cATCCCTGGCAGTTCCTGGCTGCCAATCTGGCTTTTCCCTTAACCGAAGTAAGCGCGAAGAGACAACGCCAGCGAATAAACACAAGCAAAGCCCAGCCCAAACGTAGAGATAACGCTAATCCCCCCAACCACCAGGAAGCTGCCCAACac                                                              >  2:1216421/1‑149 (MQ=255)
                                                                                          tCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGaa                                                            >  2:280181/1‑150 (MQ=255)
                                                                                            ccTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATc                                                          <  2:1323908/150‑1 (MQ=255)
                                                                                              tGGCAGTTCCTGCCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCaa                                                        >  2:81356/1‑149 (MQ=255)
                                                                                               ggCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATa                                                       <  1:1175634/150‑1 (MQ=255)
                                                                                                 cAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATaaa                                                     >  1:331846/1‑150 (MQ=255)
                                                                                                   gTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTag                                                                                                         <  2:54765/96‑1 (MQ=255)
                                                                                                   gTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAAtt                                                   >  1:745116/1‑150 (MQ=255)
                                                                                                   gTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAAtt                                                   <  1:161542/150‑1 (MQ=255)
                                                                                                   cttCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTag                                                                                                         >  1:54765/2‑96 (MQ=255)
                                                                                                      ccTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGc                                                >  2:1547579/1‑150 (MQ=255)
                                                                                                       ctggctggCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGcc                                               >  2:133696/1‑150 (MQ=255)
                                                                                                          gctggcAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgac                                            <  1:555/150‑1 (MQ=255)
                                                                                                              gcAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAg                                        <  1:166558/150‑1 (MQ=255)
                                                                                                                aaTCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGtt                                      <  1:133696/150‑1 (MQ=255)
                                                                                                                  tCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTAGCCGACGAAGTTAc                                    >  2:636157/1‑150 (MQ=255)
                                                                                                                   cTGGCTTTTCGCTTTACGGAAGTAAGCGGGATGTAACATCGCCAGCGAATAAACTGATGAAAAGACCAGCGGAAAGGTTGAGATAACGGTAATCACCGGTACAACGTGGAAGCTGCCGAACAGAATCATAAATTAGCCGACGAAGTTAcc                                   <  2:712026/150‑1 (MQ=255)
                                                                                                                   cTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTAcc                                   <  2:1402355/150‑1 (MQ=255)
                                                                                                                     ggCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCgg                                 <  1:1083811/150‑1 (MQ=255)
                                                                                                                       cTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTg                               >  2:1597431/1‑150 (MQ=255)
                                                                                                                         tttCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGcc                             <  2:1512075/150‑1 (MQ=255)
                                                                                                                         tttCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGcc                             <  2:757031/150‑1 (MQ=255)
                                                                                                                           tCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAg                           <  1:1508615/150‑1 (MQ=255)
                                                                                                                            cGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAgg                          >  1:481975/1‑150 (MQ=255)
                                                                                                                             gCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGc                         >  1:96940/1‑150 (MQ=255)
                                                                                                                                ttACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATc                      <  1:219060/150‑1 (MQ=255)
                                                                                                                                 tACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATcc                     >  2:1601028/1‑150 (MQ=255)
                                                                                                                                    cGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCaa                  >  2:936771/1‑150 (MQ=255)
                                                                                                                                       aGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTg                               <  2:276080/134‑1 (MQ=255)
                                                                                                                                       aGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTg                               >  1:276080/1‑134 (MQ=255)
                                                                                                                                       aGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCg               <  2:229546/150‑1 (MQ=255)
                                                                                                                                       aGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCg               >  1:1466386/1‑150 (MQ=255)
                                                                                                                                        gTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGt              >  1:357285/1‑150 (MQ=255)
                                                                                                                                         cgctggcgATGTTACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGtt             >  2:26813/6‑150 (MQ=255)
                                                                                                                                             ggcgATGTTACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCa         >  1:1574749/2‑150 (MQ=255)
                                                                                                                                              gcgATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCAc        <  1:1395798/150‑1 (MQ=255)
                                                                                                                                                gATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCg      <  1:1601028/150‑1 (MQ=255)
                                                                                                                                                 aTGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCGc     >  1:1377637/1‑150 (MQ=255)
                                                                                                                                                    tAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCGCaaa  <  1:236260/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGAGGTATCCAGAATCGGCTGCGGATAGAAGCCCAGCAGTACCAGCAGCACCACCAGCAACAGGATCATAAACAGCTCACGCAGCGACATCCCTGGCAGTTCCTGGCTGGCAATCTGGCTTTTCGCTTTACCGAAGTAAGCGCGATGTAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGAAGTTACCGGTGCCAGGCATCCCAAGCGTTGCCACCGCAAA  >  NC_000913/2391574‑2391871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: