Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,788 T→C 17.5% intergenic (+20/+33) oppF → / ← yciU oligopeptide ABC transporter ATPase/UPF0263 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,306,7880TC17.5% 299.3 / 58.7 137intergenic (+20/+33)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base T (61/52);  new base C (10/14);  total (71/66)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 137 total aligned reads displayed.

GAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAA  >  NC_000913/1306691‑1306886
                                                                                                 |                                                                                                  
gtgTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTg                                                                                                  <  1:1878155/98‑1 (MQ=255)
     cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                             >  2:2214331/1‑100 (MQ=255)
     cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGc                                                                                             <  2:233482/100‑1 (MQ=255)
      gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGcc                                                                                            <  1:2712930/100‑1 (MQ=255)
       ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                           >  1:3137416/1‑100 (MQ=255)
        cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCt                                                                                          <  1:1598619/100‑1 (MQ=255)
        cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                          <  1:1077892/100‑1 (MQ=255)
        cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                          <  2:1297218/100‑1 (MQ=255)
        cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                          <  1:1927027/100‑1 (MQ=255)
         aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtt                                                                                         <  1:118863/100‑1 (MQ=255)
         aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtn                                                                                         <  2:2722206/100‑2 (MQ=255)
          aaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCttt                                                                                        >  1:2072235/1‑100 (MQ=255)
           aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                       <  2:260310/100‑3 (MQ=255)
           aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                       <  2:166447/100‑3 (MQ=255)
           aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                       <  2:3791958/100‑3 (MQ=255)
             cacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattg                                                                                     <  2:1243200/100‑2 (MQ=39)
              acGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTTAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCa                                                                                    <  1:2474002/100‑1 (MQ=255)
              acGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCa                                                                                    <  1:2562550/100‑1 (MQ=255)
                tctcTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACg                                                                                  <  1:1035585/97‑1 (MQ=255)
                gTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACg                                                                                  >  2:286487/1‑100 (MQ=255)
                 tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                 >  2:645551/1‑100 (MQ=255)
                     gTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaa                                                                             <  2:740875/100‑1 (MQ=255)
                         tGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCg                                                                         >  2:2049801/1‑100 (MQ=255)
                          ggAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTGGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                        >  1:1503146/1‑100 (MQ=255)
                           gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                       <  2:159954/100‑1 (MQ=255)
                           gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                       <  2:2799651/100‑1 (MQ=255)
                            aGGGGAGTTTCGGACACTCCGCTTTTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                       <  1:1782107/99‑1 (MQ=255)
                             ggggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                     <  1:1090848/100‑1 (MQ=255)
                             ggggAGTGTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                     >  2:2708157/1‑100 (MQ=255)
                                gAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattt                                                                  <  2:40791/100‑1 (MQ=255)
                                 aGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattta                                                                 <  1:115788/100‑1 (MQ=255)
                                  gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                >  2:1768238/1‑100 (MQ=255)
                                  gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                >  1:1643312/1‑100 (MQ=255)
                                   tttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATtatttatt                                                               >  2:1773858/1‑100 (MQ=255)
                                     tCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTcc                                                             <  2:3070488/100‑1 (MQ=255)
                                       aGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTccc                                                            >  1:2928563/1‑99 (MQ=255)
                                        gACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGc                                                          <  2:1909085/100‑1 (MQ=255)
                                        gACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGc                                                          >  2:3404325/1‑100 (MQ=255)
                                          cacTCCGTTTCTTGTCTGCAAGTCGTTGCGCTTTAAAAACAATAAGGGCTGACAGTTTTCAGCACTTTTTCACGCTATCAGTGATgatttatt                                                               <  1:2000589/93‑1 (MQ=255)
                                            cTCCGTTTGTTGCCTGAAAGTCGATCCGCTTTAACAACAATAAGGGGTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                         <  1:983533/97‑1 (MQ=255)
                                            cTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAga                                                      >  1:2775080/1‑100 (MQ=255)
                                             tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                     >  2:3583628/1‑100 (MQ=255)
                                             tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                     <  1:3151888/100‑1 (MQ=255)
                                             tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                     <  2:1328244/100‑1 (MQ=255)
                                              ccGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgatg                                                    >  1:623551/1‑100 (MQ=255)
                                                 tttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGata                                                 <  1:3449801/100‑1 (MQ=255)
                                                 tttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGata                                                 <  2:691784/100‑1 (MQ=255)
                                                    cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                              <  1:1201783/100‑1 (MQ=255)
                                                    cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                              <  2:2550935/100‑1 (MQ=255)
                                                     ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGc                                             <  2:2091/100‑1 (MQ=255)
                                                     gtgtcTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGc                                             <  1:2165070/96‑1 (MQ=255)
                                                      tGCCTTAAAGTCTATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCa                                            <  1:3417018/100‑1 (MQ=255)
                                                      tGCCTGAGAGTCGATCCGCTTTACCAACAATAAGGGCTGACAGCTGCAAGACATTTTTCAAGCTCAACGCGATTATTTACTCCCGCCAGATGATATGGCa                                            >  1:2667008/1‑100 (MQ=255)
                                                        ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaaa                                          >  2:3497553/1‑100 (MQ=255)
                                                        ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCAAGATGATATGGCaaa                                          >  1:1811580/1‑100 (MQ=255)
                                                         ntGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAg                                         >  1:3680340/2‑100 (MQ=255)
                                                         cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAg                                         >  2:1528243/1‑100 (MQ=255)
                                                            aaaGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGttt                                      >  2:1212405/1‑100 (MQ=255)
                                                              aGTCGATCCGCTTTAAAAACAATAAGGGCTGAGAGTTGTAAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                    <  1:2679355/100‑1 (MQ=255)
                                                              aGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                    <  2:2239162/100‑1 (MQ=255)
                                                               gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGTAAAGTTtgtg                                   >  1:2201180/1‑100 (MQ=255)
                                                               gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                   <  2:676044/100‑1 (MQ=255)
                                                                tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGa                                  <  1:148368/100‑1 (MQ=255)
                                                                  gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                <  1:2448512/100‑1 (MQ=255)
                                                                  gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                >  1:2979552/1‑100 (MQ=255)
                                                                  gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                >  2:1174592/1‑100 (MQ=255)
                                                                   aaaaaaaatAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                   <  1:2600905/88‑1 (MQ=255)
                                                                      cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGCTGATATGGCAAAGTTTGTGATCtttt                            >  2:1372941/1‑100 (MQ=255)
                                                                      cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                            >  1:2275328/1‑100 (MQ=255)
                                                                      cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                            <  1:2990313/100‑1 (MQ=255)
                                                                      cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAACCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                            >  1:3304152/1‑100 (MQ=255)
                                                                       gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                           >  2:1982190/1‑100 (MQ=255)
                                                                       gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                           >  1:441166/1‑100 (MQ=255)
                                                                       gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                           >  2:28261/1‑100 (MQ=255)
                                                                          ttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgc                        >  2:1886733/1‑100 (MQ=255)
                                                                           tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                       >  1:249749/1‑100 (MQ=255)
                                                                           tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                       >  1:1398387/1‑100 (MQ=255)
                                                                            aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                      <  2:3472791/100‑1 (MQ=255)
                                                                             aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                     <  2:3111130/100‑1 (MQ=255)
                                                                             aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                     <  2:3112343/100‑1 (MQ=255)
                                                                              aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                    >  2:1772216/1‑100 (MQ=255)
                                                                              aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                    <  2:671541/100‑1 (MQ=255)
                                                                               gtgaaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                   >  2:623057/7‑100 (MQ=38)
                                                                                aCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATa                  <  1:3258414/100‑1 (MQ=255)
                                                                                 cAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATaa                 >  1:2986361/1‑100 (MQ=255)
                                                                                   aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa               >  2:527136/3‑100 (MQ=255)
                                                                                    tAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg              <  2:154436/100‑1 (MQ=255)
                                                                                     aaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAACGTGTGTGATCTTTTTCGCGGCATAACaaa             >  2:1474317/1‑98 (MQ=255)
                                                                                      aGGGCTGACAGTTGTCAGCCCTTTTTTTTTTTTTTATCGATTATTTATTGCCGCCAGATGATTTGTCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  1:548860/1‑100 (MQ=255)
                                                                                      aGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:577675/1‑100 (MQ=255)
                                                                                      aGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:1564178/1‑100 (MQ=255)
                                                                                      aGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  1:915148/1‑100 (MQ=255)
                                                                                      aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:1527970/1‑100 (MQ=255)
                                                                                      aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:3357870/1‑100 (MQ=255)
                                                                                      aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:1588630/1‑100 (MQ=255)
                                                                                        ggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATg          >  2:2430842/1‑100 (MQ=255)
                                                                                         gCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGc         >  1:3163752/1‑100 (MQ=255)
                                                                                          ntGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGCTGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg        >  2:453676/2‑100 (MQ=255)
                                                                                           tGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg       <  2:723571/100‑1 (MQ=255)
                                                                                            gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg      <  2:439629/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAA  >  NC_000913/1306691‑1306886

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