Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,292 C→G 5.6% A77A (GCG→GCC aspS ← aspartate‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,2920CG 5.6% 472.5 / 14.5 144A77A (GCG→GCCaspSaspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (68/68);  new base G (5/3);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

Only 100 of 145 total aligned reads displayed.

GCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACG  >  NC_000913/1950148‑1950441
                                                                                                                                                |                                                                                                                                                      
gCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACg                                                                                                                                                  <  1:542699/150‑1 (MQ=255)
  ggTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTa                                                                                                                                                >  2:77493/1‑150 (MQ=255)
  ggTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTa                                                                                                                                                >  1:221528/1‑150 (MQ=255)
     aTTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACgg                                                                                                                                             <  1:1630835/150‑1 (MQ=255)
     aTTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACgg                                                                                                                                             >  2:619740/1‑150 (MQ=255)
     aTTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATAGCGGTTAATATTTTTTTCGTCACGCGCACGTACgg                                                                                                                                             >  2:1048355/1‑150 (MQ=255)
       ttCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTg                                                                                                                                           >  1:1516423/1‑150 (MQ=255)
       ttCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTg                                                                                                                                           <  2:1355994/150‑1 (MQ=255)
         cAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGcc                                                                                                                                         <  2:1033314/150‑1 (MQ=255)
           gACGCTATTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTATCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCTCACGTACGGTGCCCg                                                                                                                                       <  2:1371547/150‑1 (MQ=255)
           gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCg                                                                                                                                       >  2:732199/1‑150 (MQ=255)
           gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCg                                                                                                                                       <  2:279997/150‑1 (MQ=255)
           gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCg                                                                                                                                       <  1:58654/150‑1 (MQ=255)
             cgcgCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                                                                     >  2:1187979/1‑150 (MQ=255)
             cgcgCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                                                                     <  2:948533/150‑1 (MQ=255)
              gcgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCTTGa                                                                                                                                    <  1:477756/150‑1 (MQ=255)
              gcgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGa                                                                                                                                    <  2:1311608/150‑1 (MQ=255)
                gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                                                                                  <  1:875807/150‑2 (MQ=255)
                gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAcc                                                                                                                                  >  1:622698/1‑150 (MQ=255)
                      tcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATg                                                                                                                            <  2:1519529/150‑1 (MQ=255)
                      tcGGGGTTGACGGGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcgaaaaaa                                                                                                                            <  1:861725/150‑8 (MQ=255)
                       cGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGc                                                                                                                           >  2:1473815/1‑150 (MQ=255)
                        ggTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATAGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                                                                          >  2:750566/1‑150 (MQ=255)
                          tgtTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                                                        >  1:680288/1‑150 (MQ=255)
                            ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                                                      <  1:1146502/150‑1 (MQ=255)
                              gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGTCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATg                                                                                                                            >  2:1436739/1‑142 (MQ=255)
                              gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGTCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATg                                                                                                                            <  1:1436739/142‑1 (MQ=255)
                               aCGTGGTTAGAGTCAAGCGGCAGAACATGTTGGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                                                                   <  1:719691/150‑1 (MQ=255)
                               aCGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                                                                   <  1:159229/150‑1 (MQ=255)
                                cGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                                                                  <  2:1599542/150‑1 (MQ=255)
                                   ggTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                                                               >  1:941720/1‑150 (MQ=255)
                                        gagTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                                                          <  1:260820/150‑1 (MQ=255)
                                          gTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTc                                                                                                        >  2:339310/1‑150 (MQ=255)
                                          gTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTc                                                                                                        >  1:1482519/1‑150 (MQ=255)
                                            cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                                                        >  1:120160/1‑132 (MQ=255)
                                            cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                                                        <  2:120160/132‑1 (MQ=255)
                                            cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                                                      >  2:470777/1‑150 (MQ=255)
                                            cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                                                      >  2:421885/1‑150 (MQ=255)
                                             aaGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                                                     <  1:1514060/150‑1 (MQ=255)
                                             aaGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                                                     >  1:361741/1‑150 (MQ=255)
                                              aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                    >  2:543611/1‑150 (MQ=255)
                                              aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                    <  2:376718/150‑1 (MQ=255)
                                               gCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                                                                   <  1:77493/150‑1 (MQ=255)
                                                 ggCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                                                                 >  1:660740/1‑150 (MQ=255)
                                                 ggCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                                                                 >  2:1073583/1‑150 (MQ=255)
                                                  gCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                                                                >  1:1337126/1‑150 (MQ=255)
                                                      aaCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                                                            <  2:75055/150‑1 (MQ=255)
                                                             gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                                                                     <  1:901989/150‑1 (MQ=255)
                                                                   ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                                                               <  1:1473815/150‑1 (MQ=255)
                                                                    tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                                                              <  1:732199/150‑1 (MQ=255)
                                                                     gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCTACACGCGCACGTACGGTGCCCGTGACATGAATGCAGAACTCATTACGCAGTAAAGAGGCCAGCTTAAACGCGTCCGCACg                                                                             >  1:1069570/1‑150 (MQ=255)
                                                                       tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                                                                           >  1:406290/1‑150 (MQ=255)
                                                                         atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                                                                         >  1:1651764/1‑150 (MQ=255)
                                                                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                                >  2:1587969/1‑106 (MQ=255)
                                                                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                    >  1:488120/1‑118 (MQ=255)
                                                                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                    <  2:488120/118‑1 (MQ=255)
                                                                                gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                                                                  >  1:1540671/1‑150 (MQ=255)
                                                                                gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                                                                  <  1:1593435/150‑1 (MQ=255)
                                                                                 cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaa                                                                 <  1:1240218/150‑1 (MQ=255)
                                                                                  gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                                                                <  1:543611/150‑1 (MQ=255)
                                                                                  gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                                                                <  1:339310/150‑1 (MQ=255)
                                                                                    ggACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                                                              >  1:1037835/1‑150 (MQ=255)
                                                                                        gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                                                          <  1:1187979/150‑1 (MQ=255)
                                                                                         ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                                                         >  1:1085915/1‑150 (MQ=255)
                                                                                          catcACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                                                        <  2:1351929/147‑1 (MQ=255)
                                                                                                ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAACTACCTGCACGAt                                                  >  2:711537/1‑150 (MQ=255)
                                                                                                ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGAt                                                  <  1:1109388/150‑1 (MQ=255)
                                                                                                 tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa                                                 >  2:1574811/1‑150 (MQ=255)
                                                                                                   gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAcc                                               <  2:509102/150‑1 (MQ=255)
                                                                                                   gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAcc                                               <  1:1578941/150‑1 (MQ=255)
                                                                                                   gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAcc                                               >  2:803380/1‑150 (MQ=255)
                                                                                                       tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCCGCTTTAACGCGACCGCACGATCCGGATCGAAAAATACCAGCACGATACCTTcg                                           >  1:1017564/1‑150 (MQ=255)
                                                                                                          ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCgg                                        <  1:77693/150‑1 (MQ=255)
                                                                                                          ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCgg                                        >  1:205393/1‑150 (MQ=255)
                                                                                                                gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCgcgc                                  <  1:413893/150‑1 (MQ=255)
                                                                                                                    tatCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCatat                              <  1:121749/150‑1 (MQ=255)
                                                                                                                     atCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc                                                     >  2:104991/1‑126 (MQ=255)
                                                                                                                     atCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATc                             <  1:1574811/150‑1 (MQ=255)
                                                                                                                       cgcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGa                           <  1:3340/150‑1 (MQ=255)
                                                                                                                         gggTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATg                         <  2:1069570/149‑1 (MQ=255)
                                                                                                                               atatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATc                   <  2:34766/150‑1 (MQ=255)
                                                                                                                               atatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATc                   <  1:359019/150‑1 (MQ=255)
                                                                                                                                  tttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAgg                <  1:750566/150‑1 (MQ=255)
                                                                                                                                   ttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc               <  1:1048355/150‑1 (MQ=255)
                                                                                                                                   ttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc               >  1:25444/1‑150 (MQ=255)
                                                                                                                                    tttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCt              <  1:1062179/150‑1 (MQ=255)
                                                                                                                                       ttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAcc           <  2:361741/150‑1 (MQ=255)
                                                                                                                                        tCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCa          >  2:83847/1‑150 (MQ=255)
                                                                                                                                         cGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCaa         >  1:1317157/1‑150 (MQ=255)
                                                                                                                                          gTCACGGGCACCGTACGTGCGCGTGTCCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGAACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg        >  1:751849/1‑150 (MQ=255)
                                                                                                                                          gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg        >  1:1285436/1‑150 (MQ=255)
                                                                                                                                          gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc               >  1:1219294/1‑143 (MQ=255)
                                                                                                                                          gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc               <  2:1219294/143‑1 (MQ=255)
                                                                                                                                          gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg        >  1:1232460/1‑150 (MQ=255)
                                                                                                                                          gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg        <  2:1540671/150‑1 (MQ=255)
                                                                                                                                           tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGa       >  1:677115/1‑150 (MQ=255)
                                                                                                                                            cACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACATGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGAt      >  2:595291/1‑150 (MQ=255)
                                                                                                                                            cACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGAATATCGATGAAGATCAGGCTACCAAGAt      >  1:182937/1‑150 (MQ=255)
                                                                                                                                             aCGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATc     <  2:1516423/150‑1 (MQ=255)
                                                                                                                                                cgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacg  <  1:483887/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACG  >  NC_000913/1950148‑1950441

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: