Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,840,141 A→G intergenic (‑148/‑74) nlpA ← / → yicS lipoprotein‑28/uncharacterized protein YicS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,840,1410AG99.0% 352.9 / ‑5.4 99intergenic (‑148/‑74)nlpA/yicSlipoprotein‑28/uncharacterized protein YicS
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (54/44);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTACAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGCATCGTATACG  >  NC_000913/3840017‑3840275
                                                                                                                            |                                                                                                                                      
tACCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGc                                                                                                                         >  2:1162461/1‑140 (MQ=255)
 aCCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCg                                                                                                                        >  1:1037549/1‑140 (MQ=255)
 aCCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCg                                                                                                                        >  1:508604/1‑140 (MQ=255)
  ccTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGa                                                                                                                       >  2:431185/1‑140 (MQ=255)
      aCTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGaa                                                                                                                   >  2:86647/1‑140 (MQ=255)
       cTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGaaa                                                                                                                  >  2:983873/1‑140 (MQ=255)
        tataGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAg                                                                                                                 >  2:881022/1‑140 (MQ=255)
          taGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGaaa                                                                                                                  >  2:147953/1‑137 (MQ=255)
          taGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGaaa                                                                                                                  <  1:147953/137‑1 (MQ=255)
           aGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCaa                                                                                                              >  1:1160954/1‑140 (MQ=255)
            gCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTACAGCATGATGAGAGCGATGAAAGCaaa                                                                                                             >  2:148321/1‑140 (MQ=255)
             cGGTAAGGTCGCTTCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAAc                                                                                                            >  2:659245/1‑140 (MQ=255)
             cGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAGAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAAc                                                                                                            >  2:1124304/1‑140 (MQ=255)
             cGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAAc                                                                                                            >  1:491811/1‑140 (MQ=255)
               gTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTa                                                                                                          <  1:955360/140‑1 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAg                                                                                                                 <  2:306007/132‑1 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAg                                                                                                                 >  1:306007/1‑132 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCaaa                                                                                                             >  1:478311/1‑136 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCaaa                                                                                                             <  2:478311/136‑1 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCaa                                                                                                              >  1:24285/1‑135 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCaa                                                                                                              <  2:24285/135‑1 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTac                                                                                                         >  1:740565/1‑140 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTac                                                                                                         >  2:628097/1‑140 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTac                                                                                                         >  1:805767/1‑140 (MQ=255)
                tAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTac                                                                                                         >  1:1126469/1‑140 (MQ=255)
                  aGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTacac                                                                                                       >  2:324074/1‑140 (MQ=255)
                  aGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTacac                                                                                                       >  2:246644/1‑140 (MQ=255)
                  aGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTacac                                                                                                       >  1:974879/1‑140 (MQ=255)
                    gTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACtt                                                                                                     <  1:628097/140‑1 (MQ=255)
                      cGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTAt                                                                                                   <  1:881022/140‑1 (MQ=255)
                       gCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCaaa                                                                                                             <  2:1169087/129‑1 (MQ=255)
                       gCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCaaa                                                                                                             >  1:1169087/1‑129 (MQ=255)
                       gCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACt                                                                                                      <  2:727752/136‑1 (MQ=255)
                       gCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACt                                                                                                      >  1:727752/1‑136 (MQ=255)
                       gCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATg                                                                                                  >  2:170476/1‑140 (MQ=255)
                       gCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATg                                                                                                  <  2:491811/140‑1 (MQ=255)
                        cTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCAATGAAAGCAAACTACACTTATGc                                                                                                 >  1:1159098/1‑140 (MQ=255)
                         tGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCt                                                                                                <  1:1162461/140‑1 (MQ=255)
                              tttAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCAtt                                                                                           >  1:272110/1‑140 (MQ=255)
                              tttAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCAtt                                                                                           <  1:86647/140‑1 (MQ=255)
                                   aaaaCGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATa                                                                                      <  2:1087636/140‑1 (MQ=255)
                                       cgaaacgCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATAcaca                                                                                  >  1:586872/1‑140 (MQ=255)
                                        gaaacgCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACAt                                                                                 <  1:983873/140‑1 (MQ=255)
                                              cTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTaa                                                                           <  1:69978/140‑1 (MQ=255)
                                                 tCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAg                                                                        >  1:176253/1‑140 (MQ=255)
                                                  cGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAgg                                                                       <  2:1160954/140‑1 (MQ=255)
                                                      aaGAATAAAAAGGAATAAAAGTGGAATATAGAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGa                                                                   <  1:1124304/140‑1 (MQ=255)
                                                      aaGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGa                                                                   <  1:659245/140‑1 (MQ=255)
                                                         aaTAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATc                                                                >  2:647615/1‑140 (MQ=255)
                                                            aaaaaGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCCATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCAt                                                             <  2:808348/140‑1 (MQ=255)
                                                               aaGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGaa                                                          <  2:974879/140‑1 (MQ=255)
                                                               aaGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGaa                                                          <  1:324074/140‑1 (MQ=255)
                                                                  gAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAAcc                                                       <  2:958447/140‑1 (MQ=255)
                                                                  gAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCAATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAAcc                                                       <  2:1159098/140‑1 (MQ=255)
                                                                  gAATAAAAGTGGAATATAAAGTCTGGCACTTCATACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAAcc                                                       <  2:958448/140‑1 (MQ=255)
                                                                   aaTAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAAcc                                                       >  1:725919/1‑139 (MQ=255)
                                                                   aaTAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAAcc                                                       <  2:725919/139‑1 (MQ=255)
                                                                   aaTAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCa                                                      >  2:943631/1‑140 (MQ=255)
                                                                       aaaGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGa                                                  <  2:814628/140‑1 (MQ=255)
                                                                        aaGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGac                                                 >  1:1002646/1‑140 (MQ=255)
                                                                         aGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGaca                                                <  1:453534/140‑1 (MQ=255)
                                                                         aGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGaca                                                <  1:647615/140‑1 (MQ=255)
                                                                          gTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGacac                                               <  2:1126469/140‑1 (MQ=255)
                                                                             gAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACactac                                            >  2:730044/1‑140 (MQ=255)
                                                                                  taAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAAcc                                                       <  2:1154478/124‑1 (MQ=255)
                                                                                  taAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAAcc                                                       >  1:1154478/1‑124 (MQ=255)
                                                                                   aaaGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTa                                      >  2:1144659/1‑140 (MQ=255)
                                                                                    aaGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTAt                                     >  2:359626/1‑140 (MQ=255)
                                                                                          ggCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCAc                               >  2:340178/1‑140 (MQ=255)
                                                                                          ggCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCAc                               >  2:528054/1‑140 (MQ=255)
                                                                                          ggCACTTCTTACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCAc                               >  2:1129048/1‑140 (MQ=255)
                                                                                                 cttACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAGGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTc                        >  1:1086788/1‑140 (MQ=255)
                                                                                                 cttACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTc                        >  2:296281/1‑140 (MQ=255)
                                                                                                  ttACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCg                       <  1:431185/140‑1 (MQ=255)
                                                                                                  ttACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCg                       <  1:943631/140‑1 (MQ=255)
                                                                                                  ttACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCg                       <  1:1144659/140‑1 (MQ=255)
                                                                                                  ttACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCg                       <  2:176253/140‑1 (MQ=255)
                                                                                                   tACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACtttt                           >  2:504324/1‑135 (MQ=255)
                                                                                                   tACTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACtttt                           <  1:504324/135‑1 (MQ=255)
                                                                                                    aCTGACCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCt                     <  1:1129048/140‑1 (MQ=255)
                                                                                                        aCCACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGc                 >  1:198803/1‑140 (MQ=255)
                                                                                                         ccACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGcc                >  2:442219/1‑140 (MQ=255)
                                                                                                         ccACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGcc                >  2:266885/1‑140 (MQ=255)
                                                                                                         ccACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGcc                <  1:296281/140‑1 (MQ=255)
                                                                                                         ccACGCAAGCACGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGcc                <  2:130533/140‑1 (MQ=255)
                                                                                                          cACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCa               >  1:1158187/1‑140 (MQ=255)
                                                                                                          cACGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCa               >  2:775410/1‑140 (MQ=255)
                                                                                                           aCGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAg              >  1:189073/1‑140 (MQ=255)
                                                                                                           aCGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAg              <  2:740565/140‑1 (MQ=255)
                                                                                                            cGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCt                     <  2:83711/132‑1 (MQ=255)
                                                                                                            cGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCt                     >  1:83711/1‑132 (MQ=255)
                                                                                                            cGCAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAgg             >  1:1144857/1‑140 (MQ=255)
                                                                                                              cAAGCGCGTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGCa           <  2:805767/140‑1 (MQ=255)
                                                                                                                  cgcgTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGc            >  1:243876/1‑135 (MQ=255)
                                                                                                                  cgcgTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGc            <  2:243876/135‑1 (MQ=255)
                                                                                                                  cgcgTTTCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGCATCGt       >  1:88503/1‑140 (MQ=255)
                                                                                                                      tttCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGCATCGTATAc   <  2:586872/140‑1 (MQ=255)
                                                                                                                       ttCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGcc                >  1:665315/1‑126 (MQ=255)
                                                                                                                       ttCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGcc                <  2:665315/126‑1 (MQ=255)
                                                                                                                       ttCGTGCAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGCATCGTATACg  >  2:64537/1‑140 (MQ=255)
                                                                                                                            |                                                                                                                                      
TACCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAACGCTATCGATAAGAATAAAAAGGAATAAAAGTGGAATATAAAGTCTGGCACTTCTTACTGACCACGCAAGCGCGTTTCGTACAGCATGATGAGAGCGATGAAAGCAAACTACACTTATGCTTCATTCGATACACATTTGTAAAAGGTTGAATCCATGAAACCAACGACACTACTGCTTATTTTCACTTTTTTCGCTATGCCAGGCATCGTATACG  >  NC_000913/3840017‑3840275

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