Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,845,901 G→A 5.2% H354Y (CAT→TAT)  topB ← DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,845,9010GA 5.2% 245.0 / 2.4 77H354Y (CAT→TAT) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base A (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

TAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAA  >  NC_000913/1845753‑1846049
                                                                                                                                                    |                                                                                                                                                    
taGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATg                                                                                                                                                      <  4:30862/149‑1 (MQ=255)
     cttcgcttcgTTCTCCGTCAGGTTGATCGCATAAATCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGCTCCGGTGCATGAacac                                                                                                                                                 >  3:266491/1‑149 (MQ=255)
      ttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACt                                                                                                                                                >  4:345603/1‑149 (MQ=255)
        cgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACt                                                                                                                                                <  2:271775/147‑1 (MQ=255)
          cttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATg                                                                                                                                            <  8:112934/149‑1 (MQ=255)
              gTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGTAATGATTGAGTGGTGCGCATCGAACTTTTTTTCATACCAACAGAGGTTGCGTATATCTGGAGCTACCACTGCGTGCGGCAACAGATCCGTGGCATGCACACTGATGGcct                                                                                                                                        >  6:236318/1‑147 (MQ=255)
                 ctcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtc                                                                                                                                      <  2:142044/148‑1 (MQ=255)
                  tcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtca                                                                                                                                     >  8:415957/1‑148 (MQ=255)
                          gTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGa                                                                                                                                              >  5:211700/1‑131 (MQ=255)
                          gTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGa                                                                                                                                              <  6:211700/131‑1 (MQ=255)
                           ttGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTg                                                                                                                           <  5:307731/149‑1 (MQ=255)
                            tGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTgg                                                                                                                          >  4:240152/1‑149 (MQ=255)
                                cGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTc                                                                                                                      >  1:362783/1‑149 (MQ=255)
                                  cAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCg                                                                                                                    <  7:166701/149‑1 (MQ=255)
                                  cAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCg                                                                                                                    <  4:171659/149‑1 (MQ=255)
                                    gAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGc                                                                                                                  <  7:75604/149‑1 (MQ=255)
                                        tCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa                                                                                                              >  3:17082/1‑149 (MQ=255)
                                          cGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAAt                                                                                                             <  7:174901/148‑1 (MQ=255)
                                                gTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGttcttc                                                                                                      <  8:54083/149‑1 (MQ=255)
                                                gTCGGAATGATGGCGTGGTGCGCATCGAACTTTTTGTCATCCCAACATCTGTTTAGTATATCTGGATATAACACTGGCTGCGGCAACAGATCCGGTTCATGAACACTGATGGCATTCATCTACTCGTGGCTTCATGCAAAATGttctta                                                                                                      <  8:102222/149‑2 (MQ=255)
                                                  cGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTg                                                                                                    >  7:147438/1‑149 (MQ=255)
                                                    gAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTgg                                                                                                   >  7:131259/1‑148 (MQ=255)
                                                            gCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGc                                                                                           <  8:93085/148‑1 (MQ=255)
                                                             cGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         <  8:153916/149‑1 (MQ=255)
                                                              gtggtgCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         >  7:75081/1‑148 (MQ=255)
                                                              ggggggcgcATCGACCTTTTTGTCATCCCAACAGCGGTAGCGTATATCTGGAGCTCCCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc                                                                                        <  5:328899/144‑1 (MQ=255)
                                                                  tgCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTTCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATc                                                                                    <  3:329855/149‑1 (MQ=255)
                                                                   gcgcATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAAt                                                                                     >  1:281830/1‑147 (MQ=255)
                                                                    cgcATCGAACTTTTTTTCATACCAACAGCGGTTGCGTATATCTGGATCTACCACTGTCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCAATCATCACCGCGTGGCGTCCTGCAAAATGTTCTTCTGGCAAATAGCGACAATCa                                                                                   <  4:266491/148‑1 (MQ=255)
                                                                      cATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa                                                                                                              <  2:281473/119‑1 (MQ=255)
                                                                      cATCGACCTTTTTGTCATCCCAACAGCGGTTGCGAATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa                                                                                                              >  1:281473/1‑119 (MQ=255)
                                                                             cTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGAt                                                                         <  8:108477/149‑1 (MQ=255)
                                                                             cTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGAt                                                                         <  3:345603/149‑1 (MQ=255)
                                                                                  tGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGt                                                                     <  7:2004/148‑1 (MQ=255)
                                                                                   gTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGa                                                                   >  5:341569/1‑149 (MQ=255)
                                                                                    tCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGa                                                                   <  8:131259/148‑1 (MQ=255)
                                                                                       tCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAg                                                               >  3:278222/1‑149 (MQ=255)
                                                                                        cccAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGc                                                              >  7:401076/1‑149 (MQ=255)
                                                                                          cAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatc                                                                                                                                   >  8:319796/1‑78 (MQ=255)
                                                                                          cAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatc                                                                                                                                   <  7:319796/78‑1 (MQ=255)
                                                                                          cAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCtt                                                            >  6:180052/1‑149 (MQ=255)
                                                                                           aaCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGGACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTtg                                                           >  7:381765/1‑149 (MQ=255)
                                                                                           aaCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         >  1:138742/1‑119 (MQ=255)
                                                                                           aaCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         <  2:138742/119‑1 (MQ=255)
                                                                                            aCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTtgt                                                          <  8:401076/149‑1 (MQ=255)
                                                                                            aCAGCGGTTGCGTAAATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTtgt                                                          <  6:327586/149‑1 (MQ=255)
                                                                                              aGCGGTTGCGAATATCTGGATCTACCACTGGCTGCGGCAACAGATCAGGTGCATGAACACTTATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGAGACAATCAGAACGCGGATAAGTGATTAGCTTGTGc                                                        <  6:297109/149‑1 (MQ=255)
                                                                                                   ttGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTc                                                   <  7:415957/149‑1 (MQ=255)
                                                                                                     gCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCg                                                  <  2:362783/148‑1 (MQ=255)
                                                                                                     gCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCg                                                  <  8:75081/148‑1 (MQ=255)
                                                                                                      cGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTa                                                <  3:240152/149‑1 (MQ=255)
                                                                                                        tataTCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTa                                                >  7:55939/1‑147 (MQ=255)
                                                                                                            tCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGttt                                          >  4:360517/1‑149 (MQ=255)
                                                                                                                gATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTgg                                      >  4:83143/1‑149 (MQ=255)
                                                                                                                    tACCACTGTCTGCGGCAACAGATCCGGTGCATTAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATATCGACAATCATAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGa                                  <  5:236318/149‑1 (MQ=255)
                                                                                                                    tACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGc                                                              <  1:409977/121‑1 (MQ=255)
                                                                                                                    tACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGc                                                              >  2:409977/1‑121 (MQ=255)
                                                                                                                      ccACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGata                                >  8:220718/1‑149 (MQ=255)
                                                                                                                      ccACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGat                                 <  6:341569/148‑1 (MQ=255)
                                                                                                                      ccACTGGCGGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGa                                  <  1:314269/147‑1 (MQ=255)
                                                                                                                       cACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGatat                               >  3:151187/1‑149 (MQ=255)
                                                                                                                        aCTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGatat                               >  2:7091/1‑148 (MQ=255)
                                                                                                                            gCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAg                           >  2:23968/1‑148 (MQ=255)
                                                                                                                             cTGCGGCAACAGATCCGGTGCATGGACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCa                         <  8:381765/149‑1 (MQ=255)
                                                                                                                             cTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCa                         <  7:330841/149‑1 (MQ=255)
                                                                                                                              tGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCAc                        >  4:27459/1‑149 (MQ=255)
                                                                                                                                    aaCAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTg                  >  2:62991/1‑149 (MQ=255)
                                                                                                                                     aCAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTgc                 <  3:27459/149‑1 (MQ=255)
                                                                                                                                       aGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTgcgc               >  4:139080/1‑149 (MQ=255)
                                                                                                                                         aTCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGat                                 <  4:243902/129‑1 (MQ=255)
                                                                                                                                         aTCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGat                                 >  3:243902/1‑129 (MQ=255)
                                                                                                                                         aTCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCAc             >  3:192383/1‑149 (MQ=255)
                                                                                                                                          tCCGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACt            >  7:163899/1‑149 (MQ=255)
                                                                                                                                           ccGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc                                                                                        >  7:372448/1‑72 (MQ=255)
                                                                                                                                           ccGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc                                                                                        <  8:372448/72‑1 (MQ=255)
                                                                                                                                            cGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCa          <  8:163899/149‑1 (MQ=255)
                                                                                                                                               tGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGAc       <  1:7091/149‑1 (MQ=255)
                                                                                                                                                  aTGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCa     >  8:23813/1‑148 (MQ=255)
                                                                                                                                                    gAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCaaaa  >  2:88528/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAA  >  NC_000913/1845753‑1846049

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