Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,998 C→A 6.4% intergenic (+19/+20) nudI → / ← ais nucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,364,9980CA 6.4% 323.0 / 10.0 94intergenic (+19/+20)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base A (3/3);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

TGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGA  >  NC_000913/2364853‑2365144
                                                                                                                                                 |                                                                                                                                                  
tGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                                                                                 <  1:751713/149‑1 (MQ=255)
 gCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                                                                <  3:227760/149‑1 (MQ=255)
 gCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                                                                >  2:334320/1‑149 (MQ=255)
   cAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTg                                                                                                                                              >  3:463550/1‑149 (MQ=255)
   cAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGAAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTTGCCACCCGAAAAACGTTGCGTTTGAAAGGTCTTCTGTAACAATGATGGCAGCCATTACTGGCTg                                                                                                                                              >  4:717599/1‑149 (MQ=255)
    aaCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGt                                                                                                                                             >  1:630912/1‑149 (MQ=255)
     aCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTc                                                                                                                                            <  3:410468/149‑1 (MQ=255)
      ccGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCa                                                                                                                                           <  4:636101/149‑1 (MQ=255)
      ccGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCa                                                                                                                                           <  3:591440/149‑1 (MQ=255)
        gagaAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGa                                                                                                                                         >  1:781107/1‑149 (MQ=255)
         agaAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                                                                        <  3:802972/149‑4 (MQ=255)
             gTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCa                                                                                                                                           <  2:267358/142‑1 (MQ=255)
             gTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCa                                                                                                                                           >  1:267358/1‑142 (MQ=255)
              tGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTCGGTTTGAAGGTCTTTTTTAACAATGATGACCGCCAATACTTGCTTTTAGGTTGgaaa                                                                                                                                  >  4:222929/1‑149 (MQ=255)
              tGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaa                                                                                                                                   >  3:780317/1‑149 (MQ=255)
                   aTAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTa                                                                                                                              <  3:551298/149‑1 (MQ=255)
                    tAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAGGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                                                                             <  2:345095/149‑1 (MQ=255)
                    tAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                                                                             <  1:732765/149‑1 (MQ=255)
                      aaCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                                                                           <  1:481128/149‑1 (MQ=255)
                         gAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                                                                        >  4:293445/1‑149 (MQ=255)
                            gagTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACCATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                                                                     <  1:625923/149‑1 (MQ=255)
                            gagTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTTGTTAAc                                                                                                                     <  4:602039/149‑1 (MQ=255)
                            gagTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                                                                     <  4:679925/149‑1 (MQ=255)
                                    ggACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTccc                                                                                                             <  2:414425/149‑1 (MQ=255)
                                      aCTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                                                                           <  3:180177/149‑1 (MQ=255)
                                            cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGa                                                                                                     <  1:348895/149‑1 (MQ=255)
                                               gggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATaa                                                                                                  <  1:573771/149‑1 (MQ=255)
                                               gggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATaa                                                                                                  >  1:615811/1‑149 (MQ=255)
                                                 gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                                <  1:187766/149‑1 (MQ=255)
                                                 gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                                >  1:21448/1‑149 (MQ=255)
                                                 gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGACACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                                <  4:320495/149‑1 (MQ=255)
                                                  tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACt                                                                                               >  2:498822/1‑149 (MQ=255)
                                                    aaaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACttt                                                                                             <  1:846361/149‑1 (MQ=255)
                                                    aaaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACttt                                                                                             <  4:780317/149‑1 (MQ=255)
                                                        cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTCGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                                                                         <  1:735258/149‑1 (MQ=255)
                                                        cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                                                                             >  1:489706/1‑113 (MQ=255)
                                                        cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                                                                             <  2:489706/113‑1 (MQ=255)
                                                        cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                                                                         >  3:382699/1‑149 (MQ=255)
                                                         tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                                                                        >  2:202363/1‑149 (MQ=255)
                                                          gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttt                                                                                       <  3:586357/149‑1 (MQ=255)
                                                           aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttt                                                                                      <  3:293445/149‑1 (MQ=255)
                                                           aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttt                                                                                      <  1:459832/149‑1 (MQ=255)
                                                             gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                                                                    <  3:772120/149‑1 (MQ=255)
                                                              aTCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGCAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                                                                   >  3:395359/1‑149 (MQ=255)
                                                               tCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCca                                                                                  <  2:711711/149‑1 (MQ=255)
                                                                  ggTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcat                                                                               >  3:18636/1‑149 (MQ=255)
                                                                       aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                                                                          >  4:864019/1‑149 (MQ=255)
                                                                            gATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaa                                                                     <  1:249081/149‑1 (MQ=255)
                                                                            gATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaa                                                                     >  3:196366/1‑149 (MQ=255)
                                                                             aTTTGAATGTCGCCACCCGAAAAACGTTGCGTTTGAAAGGTCTACTGTAACAATGATGACAGCCATTACTGGCTATCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGTCTTTTTCCACATGCATGACTaaa                                                                    <  2:791309/149‑1 (MQ=255)
                                                                             aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                                                                    <  4:159682/149‑1 (MQ=255)
                                                                              tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                   >  2:90450/1‑149 (MQ=255)
                                                                              tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                   >  4:100270/1‑149 (MQ=255)
                                                                              tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                   >  1:777964/1‑149 (MQ=255)
                                                                              tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAAGGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                   >  2:669580/1‑149 (MQ=255)
                                                                                tGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAATAAACCa                                                                 <  2:488514/149‑1 (MQ=255)
                                                                                     gTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAg                                                            >  4:772125/1‑149 (MQ=255)
                                                                                        gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATa                                                         >  3:866326/1‑149 (MQ=255)
                                                                                          cACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                                                                          <  2:474900/82‑1 (MQ=255)
                                                                                          cACCCGAAAAACGTGACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                                                                          >  1:474900/1‑82 (MQ=255)
                                                                                             ccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                                                                    <  3:716529/133‑1 (MQ=255)
                                                                                             ccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                                                                    >  4:716529/1‑133 (MQ=255)
                                                                                              cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAgg                                                   <  1:836645/149‑1 (MQ=255)
                                                                                                aaaaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt                                                 >  2:892904/1‑149 (MQ=255)
                                                                                                   aaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTaa                                              >  1:321063/1‑149 (MQ=255)
                                                                                                     cgttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                            <  1:597331/149‑1 (MQ=255)
                                                                                                       ttacgttTGAGAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGt                                          >  1:652405/1‑149 (MQ=255)
                                                                                                        tacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTc                                         <  2:717943/149‑1 (MQ=255)
                                                                                                         acgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcg                                        >  4:684254/1‑149 (MQ=255)
                                                                                                          cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgc                                       <  3:222929/149‑1 (MQ=255)
                                                                                                           gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgcg                                      >  3:115489/1‑149 (MQ=255)
                                                                                                                 aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGc                                >  2:681900/1‑149 (MQ=255)
                                                                                                                 aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGc                                >  3:387320/1‑149 (MQ=255)
                                                                                                                    gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTa                             >  1:16165/1‑149 (MQ=255)
                                                                                                                      cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATc                           <  3:630586/149‑1 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                            <  1:159764/128‑1 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                            >  2:159764/1‑128 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                       <  2:781107/149‑1 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                       <  4:611829/149‑1 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                       >  1:756476/1‑149 (MQ=255)
                                                                                                                             aaCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCa                    >  1:26202/1‑149 (MQ=255)
                                                                                                                                aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAata                 >  4:115405/1‑149 (MQ=255)
                                                                                                                                aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAata                 >  2:667533/1‑149 (MQ=255)
                                                                                                                                     tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGt            <  2:756476/149‑1 (MQ=255)
                                                                                                                                      gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                       <  3:853899/137‑1 (MQ=255)
                                                                                                                                      gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                       >  4:853899/1‑137 (MQ=255)
                                                                                                                                      gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTc           >  4:467309/1‑149 (MQ=255)
                                                                                                                                       aCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCa          <  4:866326/149‑1 (MQ=255)
                                                                                                                                          gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGc       >  4:568943/1‑149 (MQ=255)
                                                                                                                                            cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCaa     >  4:373352/1‑149 (MQ=255)
                                                                                                                                             aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAAt    <  3:100270/149‑1 (MQ=255)
                                                                                                                                              ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCatga   <  1:613583/149‑4 (MQ=255)
                                                                                                                                              ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATg   >  1:440698/1‑149 (MQ=255)
                                                                                                                                               tACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGa  <  4:18636/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGA  >  NC_000913/2364853‑2365144

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