Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,724,412 A→T 5.7% A95A (GCA→GCT yfiM → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,724,4120AT 5.7% 248.1 / 5.5 70A95A (GCA→GCTyfiMputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base T (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

ACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCA  >  NC_000913/2724265‑2724560
                                                                                                                                                   |                                                                                                                                                    
acaGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCaa                                                                                                                                                     >  3:438256/1‑149 (MQ=255)
      tCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCg                                                                                                                                               <  4:918891/149‑1 (MQ=255)
          gggATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCta                                                                                                                                           >  1:1065892/1‑149 (MQ=255)
               gAGCCGGGACCGCAGTGCCATGTTTGGATTGATGTTGTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCg                                                                                                                                      <  3:553931/149‑1 (MQ=255)
                aGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGt                                                                                                                                     >  2:594037/1‑149 (MQ=255)
                  ccGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                   <  4:110741/149‑1 (MQ=255)
                   cGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTg                                                                                                                                  >  3:1109027/1‑149 (MQ=255)
                       aTCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCaa                                                                                                                              >  1:672302/1‑149 (MQ=255)
                        tCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAAc                                                                                                                             <  1:841991/149‑1 (MQ=255)
                                  aTGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCa                                                                                                                   <  2:1065892/149‑1 (MQ=255)
                                       tGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTaaa                                                                                                              <  3:253748/149‑1 (MQ=255)
                                       tGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTaaa                                                                                                              <  1:809349/149‑1 (MQ=255)
                                                   ctctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCccct                                                                                                  >  1:345491/1‑149 (MQ=255)
                                                    tctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTc                                                                                                                                    <  1:781814/114‑1 (MQ=255)
                                                    tctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTc                                                                                                                                    >  2:781814/1‑114 (MQ=255)
                                                     ctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttc                                                                                                <  1:594037/149‑1 (MQ=255)
                                                      tGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttcc                                                                                               >  3:1009439/1‑149 (MQ=255)
                                                           gTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTg                                                                                          <  4:438256/149‑1 (MQ=255)
                                                                     gTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCa                                                                                <  1:1136363/149‑1 (MQ=255)
                                                                      tCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCAt                                                                               <  4:1009439/149‑1 (MQ=255)
                                                                        aaaaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGc                                                                                                                                             >  3:104584/1‑85 (MQ=255)
                                                                        aaaaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTACAAGCACCGGc                                                                                                                                             <  4:104584/85‑1 (MQ=255)
                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGAc                                                                                                                         >  4:593771/1‑103 (MQ=255)
                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGAc                                                                                                                         <  3:593771/103‑1 (MQ=255)
                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATc                                                                                                      <  3:1090093/122‑1 (MQ=255)
                                                                          aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGAGGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATc                                                                                                      >  4:1090093/1‑122 (MQ=255)
                                                                                tttGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAcc                                                                     >  4:951102/1‑149 (MQ=255)
                                                                                   gggATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCg                                                                                                                                               >  1:871677/1‑72 (MQ=255)
                                                                                   gggATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCg                                                                                                                                               <  2:871677/72‑1 (MQ=255)
                                                                                            gccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccacc                                                         >  4:480718/1‑149 (MQ=255)
                                                                                             ccccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccacca                                                        >  2:61836/1‑149 (MQ=255)
                                                                                             ccccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccacca                                                        >  4:43700/1‑149 (MQ=255)
                                                                                              cccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGa                                                                                       >  3:687980/1‑117 (MQ=255)
                                                                                              cccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGa                                                                                       <  4:687980/117‑1 (MQ=255)
                                                                                              cccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCcaccac                                                       >  3:274971/1‑149 (MQ=255)
                                                                                               ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGcg                                                                                                                                            <  2:486588/63‑2 (MQ=38)
                                                                                               ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGcg                                                                                                                                            >  1:486588/1‑62 (MQ=38)
                                                                                               ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCgg                                                                                                                                              <  2:892158/61‑1 (MQ=255)
                                                                                               ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCgg                                                                                                                                              >  1:892158/1‑61 (MQ=255)
                                                                                               ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACg                                                      >  3:387730/1‑149 (MQ=255)
                                                                                                cGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACgg                                                     >  4:434751/1‑149 (MQ=255)
                                                                                                     ggAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGTCCATc                                                >  4:580410/1‑149 (MQ=255)
                                                                                                        gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAg                                             >  2:513701/1‑149 (MQ=255)
                                                                                                        gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAg                                             >  1:599191/1‑149 (MQ=255)
                                                                                                             tggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGt                                           <  4:70767/146‑1 (MQ=255)
                                                                                                             tggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGt                                           >  3:70767/1‑146 (MQ=255)
                                                                                                                    ggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCag                                 >  1:336637/1‑149 (MQ=255)
                                                                                                                     gaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGt                                                                                    >  3:584673/1‑97 (MQ=255)
                                                                                                                     gaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGt                                                                                    <  4:584673/97‑1 (MQ=255)
                                                                                                                     gaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCaga                                >  1:724178/1‑149 (MQ=255)
                                                                                                                       aGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAc                                                                                                           <  4:508927/72‑1 (MQ=255)
                                                                                                                       aGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAc                                                                                                           >  3:508927/1‑72 (MQ=255)
                                                                                                                          aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                           <  3:43700/149‑1 (MQ=255)
                                                                                                                          aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                           >  4:738158/1‑149 (MQ=255)
                                                                                                                          aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                           <  1:61836/149‑1 (MQ=255)
                                                                                                                             tGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCt                        <  2:1014913/149‑1 (MQ=255)
                                                                                                                               gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACt                                                                                                                            <  1:1101351/47‑1 (MQ=255)
                                                                                                                               gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACt                                                                                                                            >  2:1101351/1‑47 (MQ=255)
                                                                                                                               gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGt                      <  4:1109027/149‑1 (MQ=255)
                                                                                                                                   gggATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTcc                  <  2:672302/149‑1 (MQ=255)
                                                                                                                                     gATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCAt                <  1:513701/149‑1 (MQ=255)
                                                                                                                                       tGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTc              <  3:738158/149‑1 (MQ=255)
                                                                                                                                       tGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTc              >  2:478922/1‑149 (MQ=255)
                                                                                                                                          cGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTa           >  4:600167/1‑149 (MQ=255)
                                                                                                                                          agCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTa           >  4:45636/2‑149 (MQ=255)
                                                                                                                                           gCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCagaa                               >  4:902508/1‑128 (MQ=255)
                                                                                                                                           gCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCagaa                               <  3:902508/128‑1 (MQ=255)
                                                                                                                                           gCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt          >  2:1090920/1‑149 (MQ=255)
                                                                                                                                            ccGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAtt         >  1:189236/1‑149 (MQ=255)
                                                                                                                                              ggTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGa       <  3:45636/149‑1 (MQ=255)
                                                                                                                                                 gCAAGCACCGGCTATTCCGTCTGGCAAATGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAttt    <  3:580410/149‑1 (MQ=255)
                                                                                                                                                   aaGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGCAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCa  >  3:276285/1‑149 (MQ=255)
                                                                                                                                                   aaGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCa  >  1:571079/1‑149 (MQ=255)
                                                                                                                                                   aaGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCa  >  3:603060/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCA  >  NC_000913/2724265‑2724560

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