Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,818 C→T 6.5% C476Y (TGC→TAC)  ydiU ← UPF0061 family protein YdiU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,789,8180CT 6.5% 188.6 / 2.5 61C476Y (TGC→TAC) ydiUUPF0061 family protein YdiU
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/2);  total (31/30)
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.63e-01

CGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAAC  >  NC_000913/1789669‑1789967
                                                                                                                                                     |                                                                                                                                                     
cGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGagc                                                                                                                                                       <  2:1777272/150‑1 (MQ=255)
  cTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGagcag                                                                                                                                                     <  2:525779/150‑1 (MQ=255)
    tCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCt                                                                                                                                                   <  2:402479/150‑1 (MQ=255)
     ccccGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTg                                                                                                                                                  >  2:1897407/1‑150 (MQ=255)
            gCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCa                                                                                                                                           <  1:1807495/150‑1 (MQ=255)
             ccTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAg                                                                                                                                          <  2:1496793/150‑1 (MQ=255)
              cTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGc                                                                                                                                         <  1:1073374/150‑1 (MQ=255)
               tATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGcc                                                                                                                                        >  2:1694556/1‑150 (MQ=255)
                   cGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCCAGCAGCTGACTTCCAGCCGttt                                                                                                                                    <  2:151288/150‑1 (MQ=255)
                      cGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAcc                                                                                                                                 <  1:775246/150‑1 (MQ=255)
                       gAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAccc                                                                                                                                <  2:1110698/150‑1 (MQ=255)
                        aaTTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAcccc                                                                                                                               <  2:1984172/150‑1 (MQ=255)
                         aTTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCa                                                                                                                              >  1:303094/1‑150 (MQ=255)
                            tCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTc                                                                                                                           <  1:1049680/150‑1 (MQ=255)
                            tCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGAACATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTc                                                                                                                           >  1:911660/1‑150 (MQ=255)
                             ccagccagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCa                                                                                                                          >  2:768021/1‑150 (MQ=255)
                              cagccagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAg                                                                                                                         >  1:1351132/1‑150 (MQ=255)
                               agccagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAgg                                                                                                                        >  1:831162/1‑150 (MQ=255)
                                gccagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGTAGCTGACTTCCAGCCGTTTACCCCAGTCAGGt                                                                                                                       >  1:1809984/1‑150 (MQ=255)
                                gccagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGTAGCTGACTTCCAGCCGTTTACCCCAGTCAGGt                                                                                                                       >  1:1802688/1‑150 (MQ=255)
                                gccagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGt                                                                                                                       >  1:1503468/1‑150 (MQ=255)
                                 ccagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTg                                                                                                                      <  2:1184167/150‑1 (MQ=255)
                                   agACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTg                                                                                                                                                  >  1:1160739/1‑120 (MQ=255)
                                   agACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTg                                                                                                                                                  <  2:1160739/120‑1 (MQ=255)
                                           ttAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGAc                                                                                                            <  1:1483910/150‑1 (MQ=255)
                                               gTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAg                                                                                                        <  1:805305/150‑1 (MQ=255)
                                                    gATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGtcatc                                                                                                   >  2:716331/1‑150 (MQ=255)
                                                         cAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCAtctc                                                                                              >  1:1343900/1‑150 (MQ=255)
                                                         cAAACTGATTATCGATCATGACGGGAAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGCCTGACATAGTCATCAtctc                                                                                              >  2:1830459/1‑150 (MQ=255)
                                                           aaCTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTg                                                                                            <  1:784013/150‑1 (MQ=255)
                                                             cTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTc                                                                                          >  1:191474/1‑150 (MQ=255)
                                                              tGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCg                                                                                         <  1:1738053/150‑1 (MQ=255)
                                                                  tATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGa                                                                                     >  1:1160613/1‑150 (MQ=255)
                                                                      gATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAcc                                                                                                                                 <  1:1503824/102‑1 (MQ=255)
                                                                      gATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAcc                                                                                                                                 >  2:1503824/1‑102 (MQ=255)
                                                                      gATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAgg                                                                                 >  1:1497403/1‑150 (MQ=255)
                                                                             aCGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTTGCAAGCTGCTGGTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGc                                                                          <  2:1711110/150‑1 (MQ=255)
                                                                             aCGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCACAGTAAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGc                                                                          >  1:1741311/1‑150 (MQ=255)
                                                                                ggCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAAc                                                                       <  2:887969/150‑1 (MQ=255)
                                                                                      aaTATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAGGGATTTCGCAACGCCTCa                                                                 <  2:1183689/150‑1 (MQ=255)
                                                                                              cTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCg                                                         >  2:1685876/1‑150 (MQ=255)
                                                                                                ccGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAcccc                                                                                                                               <  1:1819790/78‑1 (MQ=255)
                                                                                                ccGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAcccc                                                                                                                               >  2:1819790/1‑78 (MQ=255)
                                                                                                  gTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGc                                                     <  2:1351132/150‑1 (MQ=255)
                                                                                                        cTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTcc                                               <  2:191474/150‑1 (MQ=255)
                                                                                                                 gCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCa                                      <  1:1897407/150‑1 (MQ=255)
                                                                                                                  cAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCAc                                     >  1:1944494/1‑150 (MQ=255)
                                                                                                                      ctgctgTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCtt                                 >  2:582904/1‑150 (MQ=255)
                                                                                                                      ctgctgTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCtt                                 >  1:373471/1‑150 (MQ=255)
                                                                                                                         ctgTTGACCAAGTAGCCTTTAACTCGAGTAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTc                              <  2:1802688/150‑1 (MQ=255)
                                                                                                                         ctgTTGACCAAGTAGCCTTTAACTCGAGTAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTc                              <  2:1809984/150‑1 (MQ=255)
                                                                                                                          tgTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCt                             <  2:831162/150‑1 (MQ=255)
                                                                                                                                 cAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCtt                      >  2:1581058/1‑150 (MQ=255)
                                                                                                                                       gCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCg                <  2:1874896/150‑1 (MQ=255)
                                                                                                                                             aaCTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGtt          >  2:1873769/1‑150 (MQ=255)
                                                                                                                                              aCTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCt                       >  2:951334/1‑136 (MQ=255)
                                                                                                                                              aCTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCt                       <  1:951334/136‑1 (MQ=255)
                                                                                                                                               cTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTgc        >  2:25773/1‑150 (MQ=255)
                                                                                                                                               cTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTgc        >  2:1749097/1‑150 (MQ=255)
                                                                                                                                               cTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTgc        >  2:160166/1‑150 (MQ=255)
                                                                                                                                                  gAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTgcgc      <  2:1105490/149‑1 (MQ=255)
                                                                                                                                                  gAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTgcgc      >  1:1105490/1‑149 (MQ=255)
                                                                                                                                                  gAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGcc     <  2:1175599/150‑1 (MQ=255)
                                                                                                                                                     cagcTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAAc  <  1:227718/150‑1 (MQ=255)
                                                                                                                                                     cagcTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTAGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTgcgcaact  >  2:661562/1‑146 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATCACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAAC  >  NC_000913/1789669‑1789967

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: