Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,836,136 G→C 5.4% D22H (GAC→CAC)  ynjB → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,836,1360GC 5.4% 401.8 / 11.5 112D22H (GAC→CAC) ynjBputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (53/53);  new base C (3/3);  total (56/56)
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.38e-01

Only 100 of 113 total aligned reads displayed.

TGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGAAGATTGTCCGTCTGGC  >  NC_000913/1835990‑1836284
                                                                                                                                                  |                                                                                                                                                    
   cAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTgg                                                                                                                                                 <  1:318937/149‑1 (MQ=255)
     gCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCa                                                                                                                                               >  4:30281/1‑149 (MQ=255)
      cccGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCGATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCaa                                                                                                                                              >  3:756033/1‑149 (MQ=255)
      cccGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCaa                                                                                                                                              >  1:1020182/1‑149 (MQ=255)
      cccGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCaa                                                                                                                                              >  3:890737/1‑149 (MQ=255)
       ccGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAg                                                                                                                                             >  1:913938/1‑149 (MQ=255)
         gTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCatgtg                                                                                                                                           <  3:575893/149‑5 (MQ=255)
          ttttAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaac                                                                                                                                          >  1:1110605/1‑148 (MQ=255)
          ttttAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaaa                                                                                                                                          <  2:637980/149‑1 (MQ=255)
          ttttAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaaa                                                                                                                                          <  2:1110605/149‑1 (MQ=255)
           tttAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACTTGCCAGTCACTGGCatgtgtt                                                                                                                                         >  3:428758/1‑143 (MQ=255)
           tttAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatgtgtt                                                                                                                                         <  1:597830/149‑7 (MQ=255)
              aaTTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTaa                                                                                                                                      <  3:30281/149‑1 (MQ=255)
                ttCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTaaaa                                                                                                                                    >  3:863934/1‑149 (MQ=255)
                   tGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATg                                                                                                                                 >  1:235987/1‑149 (MQ=255)
                      cATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAgg                                                                                                                              <  2:17233/149‑1 (MQ=255)
                       aTGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAg                                                                                                                               <  3:38846/147‑1 (MQ=255)
                       aTGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAg                                                                                                                               >  4:38846/1‑147 (MQ=255)
                       aTGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGc                                                                                                                             >  1:37889/1‑149 (MQ=255)
                           aCATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGa                                                                                                                                >  3:501880/1‑142 (MQ=255)
                           aCATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGa                                                                                                                                <  4:501880/142‑1 (MQ=255)
                           aCATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCaaa                                                                                                                         >  4:617670/1‑149 (MQ=255)
                             aTTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGa                                                                                                                                >  3:680392/1‑140 (MQ=255)
                             aTTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGa                                                                                                                                <  3:962382/140‑1 (MQ=255)
                             aTTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGa                                                                                                                                <  4:680392/140‑1 (MQ=255)
                             aTTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAg                                                                                                                               <  4:8418/141‑1 (MQ=255)
                                       aGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTc                                                                                                             >  3:189006/1‑149 (MQ=255)
                                        ggggCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCt                                                                                                            >  4:564282/1‑149 (MQ=255)
                                        ggggCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCat                                                                                                                                              <  3:720919/115‑2 (MQ=255)
                                        ggggCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATAGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACCCATGCCAGTCACTGGCat                                                                                                                                              >  4:720919/1‑114 (MQ=255)
                                         gggCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTg                                                                                                           <  4:863934/149‑1 (MQ=255)
                                         gggCTGTGCCGTATACCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTg                                                                                                           >  3:81175/1‑149 (MQ=255)
                                           gCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGa                                                                                                                                <  1:294374/126‑1 (MQ=255)
                                           gCTGTGCCGTATTCCTGAAAAACGATCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTATGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGa                                                                                                                                >  2:294374/1‑126 (MQ=255)
                                            cTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTTTGGtt                                                                                                        <  4:86812/149‑1 (MQ=255)
                                              gtgCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTa                                                                                                      >  4:375764/1‑149 (MQ=255)
                                                    tATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCt                                                                                                <  1:963836/149‑1 (MQ=255)
                                                      ttCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGcc                                                                                                                            >  1:927969/1‑119 (MQ=255)
                                                      ttCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGcc                                                                                                                            <  2:927969/119‑1 (MQ=255)
                                                          tGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcg                                                                                          >  4:830911/1‑149 (MQ=255)
                                                          tGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcg                                                                                          >  1:555115/1‑149 (MQ=255)
                                                            aaaaaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcggc                                                                                        <  3:121862/149‑1 (MQ=255)
                                                            aaaaaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcggc                                                                                        <  1:993300/149‑1 (MQ=255)
                                                            aaaaaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcggc                                                                                        >  4:325289/1‑149 (MQ=255)
                                                                      ttAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCaa                                                                              >  2:730987/1‑149 (MQ=255)
                                                                        aaGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTgg                                                                                              <  4:490316/131‑1 (MQ=255)
                                                                        aaGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTgg                                                                                              >  3:490316/1‑131 (MQ=255)
                                                                         aGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGCTTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTa                                                                           >  4:834982/1‑149 (MQ=255)
                                                                          ggAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTaa                                                                          <  2:555115/149‑1 (MQ=255)
                                                                          gaagagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTaa                                                                          >  1:910195/3‑149 (MQ=255)
                                                                           gagagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGc                                                                                <  4:368975/142‑1 (MQ=255)
                                                                           gagagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGc                                                                                >  3:368975/1‑142 (MQ=255)
                                                                            agagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAAcc                                                                        >  3:326114/1‑149 (MQ=255)
                                                                               gatgatGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcgg                                                                                         >  4:1026382/1‑129 (MQ=255)
                                                                               gatgatGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcgg                                                                                         <  3:1026382/129‑1 (MQ=255)
                                                                               gatgatGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCt                                                                     >  2:431360/1‑149 (MQ=255)
                                                                                  gatgCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTAtc                                                                  >  1:214054/1‑149 (MQ=255)
                                                                                   atgCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTAtct                                                                 <  3:375764/149‑1 (MQ=255)
                                                                                   atgCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTAtct                                                                 >  2:96088/1‑149 (MQ=255)
                                                                                    tgCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTAtctc                                                                >  3:639532/1‑149 (MQ=255)
                                                                                       gcCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCTAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGAc                                                             <  2:801716/149‑1 (MQ=255)
                                                                                          aTTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAATAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTgg                                                          <  3:1006237/149‑1 (MQ=255)
                                                                                          aTTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTgg                                                          <  2:519460/149‑1 (MQ=255)
                                                                                           ttGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTggg                                                         <  2:214054/149‑1 (MQ=255)
                                                                                              tgggtggTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGt                                                        >  1:590191/1‑147 (MQ=255)
                                                                                              tgggtggTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGt                                                        <  2:590191/147‑1 (MQ=255)
                                                                                               gggtggTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAg                                                     >  4:803306/1‑149 (MQ=255)
                                                                                               gggtggTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAg                                                     >  4:588801/1‑149 (MQ=255)
                                                                                                 gtggtTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTgg                                                                                                          >  2:435085/1‑94 (MQ=255)
                                                                                                 gtggtTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTgg                                                                                                          <  1:435085/94‑1 (MQ=255)
                                                                                                 gtggtTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGcc                                                                                                 >  3:644936/1‑103 (MQ=255)
                                                                                                 gtggtTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGcc                                                                                                 <  4:644936/103‑1 (MQ=255)
                                                                                                    gtTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAgcggcg                                                >  3:289228/1‑149 (MQ=255)
                                                                                                      tGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCgag                                              <  3:325289/149‑1 (MQ=255)
                                                                                                       gCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCgaga                                             <  4:335567/149‑1 (MQ=255)
                                                                                                       gCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCgaga                                             <  4:890737/149‑1 (MQ=255)
                                                                                                       gCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCgaga                                             <  2:910195/149‑1 (MQ=255)
                                                                                                        cTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGACGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGAt                                            >  3:945395/1‑149 (MQ=255)
                                                                                                           ggATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGaa                                         >  1:351574/1‑149 (MQ=255)
                                                                                                              ttGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAAc                                      <  2:913938/149‑1 (MQ=255)
                                                                                                              ttGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAAc                                      <  1:96088/149‑1 (MQ=255)
                                                                                                                 ttATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAAccaa                                   >  3:17418/1‑147 (MQ=255)
                                                                                                                 ttATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCa                                   >  3:1031717/1‑149 (MQ=255)
                                                                                                                  tATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCAt                                  <  1:634684/149‑1 (MQ=255)
                                                                                                                    tCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTa                                >  1:1041618/1‑149 (MQ=255)
                                                                                                                       cTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGc                             <  3:617670/149‑1 (MQ=255)
                                                                                                                        tGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCt                            <  4:189006/149‑1 (MQ=255)
                                                                                                                          ttttCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGc                                                                      >  4:1138791/1‑105 (MQ=255)
                                                                                                                          ttttCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCtat                          >  1:279399/1‑149 (MQ=255)
                                                                                                                              ctctGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGc                                                                                                                     >  1:123602/1‑54 (MQ=255)
                                                                                                                              ctctGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGc                                                                                                                     <  2:123602/54‑1 (MQ=255)
                                                                                                                                   gCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGaa                 >  2:717828/1‑149 (MQ=255)
                                                                                                                                   gCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGaa                 <  4:1031717/149‑1 (MQ=255)
                                                                                                                                   gCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGaa                 <  4:17418/149‑1 (MQ=255)
                                                                                                                                       ttcttGCCAGTCACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGAAGAtt             >  3:791470/5‑149 (MQ=255)
                                                                                                                                       cacaTGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGAAGAtt             <  2:1020182/149‑1 (MQ=255)
                                                                                                                                         caTGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGAAGATTGt           <  1:431360/149‑1 (MQ=255)
                                                                                                                                          aTGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTaa                                                                          >  2:30110/1‑85 (MQ=255)
                                                                                                                                           tGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGAAGATTGTcc         <  2:279399/149‑1 (MQ=255)
                                                                                                                                                  gACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGAAGATTGTCCGTCTGGc  <  3:564282/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGGGTGAGCGGCGAGATGAAAACCCATTACGCTATAAACCTGAAGATTGTCCGTCTGGC  >  NC_000913/1835990‑1836284

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