Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,566,645 T→A 5.1% K201I (AAA→ATA) ‡ yjiN ← DUF445 domain‑containing protein YjiN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,566,6450TA 5.1% 259.5 / 3.0 79K201I (AAA→ATA) ‡yjiNDUF445 domain‑containing protein YjiN
Reads supporting (aligned to +/- strand):  ref base T (31/44);  new base A (2/2);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAG  >  NC_000913/4566497‑4566793
                                                                                                                                                    |                                                                                                                                                    
aCTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGAttt                                                                                                                                                     <  1:423635/150‑1 (MQ=255)
    atatCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATcg                                                                                                                                                 <  1:1141734/150‑1 (MQ=255)
     tatCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATcgc                                                                                                                                                >  1:661048/1‑150 (MQ=255)
      atcatcAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATcgcg                                                                                                                                               >  2:161949/1‑150 (MQ=255)
      atcatcAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATcgcg                                                                                                                                               >  2:720487/1‑150 (MQ=255)
       tcatcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCTTTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCgcgc                                                                                                                                              <  1:1250034/150‑1 (MQ=255)
       tcatcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCgcgc                                                                                                                                              >  1:1627018/1‑150 (MQ=255)
       tcatcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgc                                                                                                                                              >  1:1887379/1‑150 (MQ=255)
       tcatcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCACGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCgcgc                                                                                                                                              >  1:318735/1‑150 (MQ=255)
        catcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCt                                                                                                                                             >  1:1342444/1‑150 (MQ=255)
          tcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTgg                                                                                                                                           <  2:744941/150‑1 (MQ=255)
          tcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaag                                                                                                                                           <  1:1075347/150‑4 (MQ=255)
           caAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGga                                                                                                                                          <  2:1548288/150‑1 (MQ=255)
           caAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaagt                                                                                                                                          <  1:904885/150‑5 (MQ=255)
             aGCAAAGAATTCACCGCGTCAGAAACCACCTCCGCGCTATGTTCGCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGaa                                                                                                                                       <  1:873664/150‑1 (MQ=255)
                      ttCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATc                                                                                                                               >  1:117259/1‑150 (MQ=255)
                      ttCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATc                                                                                                                               <  2:1393349/150‑1 (MQ=255)
                       tCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCa                                                                                                                              <  1:302557/150‑1 (MQ=255)
                        cACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCaa                                                                                                                             >  2:616255/1‑150 (MQ=255)
                        cACCGCGTCAGAAACCAACTCCGCGCTATGCTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCaa                                                                                                                             >  1:536617/1‑150 (MQ=255)
                         aCCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAAc                                                                                                                            <  2:708234/150‑1 (MQ=255)
                          ccGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACt                                                                                                                           >  1:1150514/1‑150 (MQ=255)
                           cgcgTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtg                                                                                                                          <  1:161949/150‑1 (MQ=255)
                           cgcgTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtg                                                                                                                          <  2:1080300/150‑1 (MQ=255)
                             cgTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtgtg                                                                                                                        <  1:492624/150‑1 (MQ=255)
                                cAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGa                                                                                                                     >  1:1675551/1‑150 (MQ=255)
                                   aaaCCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCa                                                                                                                  <  2:1381325/150‑1 (MQ=255)
                                   aaaCCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCa                                                                                                                  >  2:441510/1‑150 (MQ=255)
                                      ccAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCg                                                                                                               <  1:11905/150‑1 (MQ=255)
                                         aCTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTAt                                                                                                            <  2:1627018/150‑1 (MQ=255)
                                         aCTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTAt                                                                                                            >  1:753147/1‑150 (MQ=255)
                                         aCTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCACGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTAt                                                                                                            <  2:318735/150‑1 (MQ=255)
                                           tCCGCGCTATGTTCGCCCAACCATTCGTTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCTAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATca                                                                                                          <  1:1884525/150‑2 (MQ=255)
                                            ccGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCca                                                                                                         >  2:501826/1‑150 (MQ=255)
                                              gcgcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagc                                                                                                       <  1:720487/150‑1 (MQ=255)
                                              gcgcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagc                                                                                                       <  1:1925866/150‑1 (MQ=255)
                                              gcgcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagc                                                                                                       >  2:1251008/1‑150 (MQ=255)
                                               cgcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagca                                                                                                      <  2:69596/150‑1 (MQ=255)
                                                gcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagcag                                                                                                     >  2:1887997/1‑150 (MQ=255)
                                                        cGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgat                                                                                             >  1:1199445/1‑150 (MQ=255)
                                                        cGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgat                                                                                             >  1:903408/1‑150 (MQ=255)
                                                         gCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatg                                                                                            >  2:1365211/1‑150 (MQ=255)
                                                          cccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatga                                                                                           >  1:373055/1‑150 (MQ=255)
                                                          cccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatga                                                                                           >  2:670764/1‑150 (MQ=255)
                                                           ccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGAc                                                                                          <  2:1345037/150‑1 (MQ=255)
                                                              aCCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGAt                                                                                       <  1:671669/150‑1 (MQ=255)
                                                                   tCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGttt                                                                                  <  2:1342444/150‑1 (MQ=255)
                                                                     ggTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGttttt                                                                                <  2:1567929/150‑1 (MQ=255)
                                                                      gTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCTAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTg                                                                               <  1:1681906/150‑1 (MQ=255)
                                                                          ggAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCa                                                                           >  2:649331/1‑150 (MQ=255)
                                                                             gaATTTTGGCTTTCAGTGGATGCTGGCTCTCGAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATAc                                                                        <  2:1829319/150‑1 (MQ=255)
                                                                               aTTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATActc                                                                      <  2:737433/150‑1 (MQ=255)
                                                                                   tGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCa                                                                  <  2:1493016/150‑1 (MQ=255)
                                                                                      cTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCa                                                               <  2:1150514/150‑1 (MQ=255)
                                                                                       tttCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCAt                                                              <  2:739534/150‑1 (MQ=255)
                                                                                          cAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCaa                                                           >  1:1058199/1‑150 (MQ=255)
                                                                                          cAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCaa                                                           >  2:1196692/1‑150 (MQ=255)
                                                                                             tGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGc                                                        <  1:649331/150‑1 (MQ=255)
                                                                                                  gctcgctcTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTgg                                                   <  1:441510/150‑1 (MQ=255)
                                                                                                     cgctcTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGc                                                >  2:1756522/1‑150 (MQ=255)
                                                                                                           ccagccagCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGaaa                                          >  1:437546/1‑150 (MQ=255)
                                                                                                            cagccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAg                                         <  2:1503558/150‑1 (MQ=255)
                                                                                                            cagccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAg                                         <  2:1675551/150‑1 (MQ=255)
                                                                                                            cagccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAg                                         >  2:1778499/1‑150 (MQ=255)
                                                                                                               ccagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATc                                      <  2:727361/150‑1 (MQ=255)
                                                                                                                cagcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTATGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCg                                     <  2:246173/150‑1 (MQ=255)
                                                                                                                  gcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGAc                                   <  1:670764/150‑1 (MQ=255)
                                                                                                                   cGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGAcc                                  <  2:903408/150‑1 (MQ=255)
                                                                                                                   cGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGAcc                                  <  1:1887997/150‑1 (MQ=255)
                                                                                                                     aaCAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCtt                                >  1:361071/1‑150 (MQ=255)
                                                                                                                       cAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTAt                              >  1:1824519/1‑150 (MQ=255)
                                                                                                                              cTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatga                                                                                           <  2:1050887/82‑1 (MQ=255)
                                                                                                                              cTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatga                                                                                           >  1:1050887/1‑82 (MQ=255)
                                                                                                                                  gCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATAAAATGTGGGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGa                   <  2:437546/150‑1 (MQ=255)
                                                                                                                                     aTAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATgg                <  2:1042919/150‑1 (MQ=255)
                                                                                                                                     aTAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATgg                <  1:1196692/150‑1 (MQ=255)
                                                                                                                                     aTAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATgg                <  2:1175224/150‑1 (MQ=255)
                                                                                                                                          ccagcgcGATAAATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcgc           >  1:1470883/4‑150 (MQ=255)
                                                                                                                                            tGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcgcgc         <  1:501826/150‑1 (MQ=255)
                                                                                                                                              cgcgATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCgcgcgc       >  1:252918/1‑150 (MQ=255)
                                                                                                                                                   tttATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAg  <  1:1365211/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAG  >  NC_000913/4566497‑4566793

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