Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 432,193 +G intergenic (‑180/+119) tsx ← / ← yajI nucleoside channel, receptor of phage T6 and colicin K/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913432,1891.G91.2% 261.1 / 20.4 81intergenic (‑176/+123)tsx/yajInucleoside channel, receptor of phage T6 and colicin K/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base G (41/32);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAGCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑G‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAGAT  >  NC_000913/432090‑432277
                                                                                                    |                                                                                            
aggagCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTcgggc                                                                                               <  1:319065/100‑1 (MQ=255)
 ggagCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTcgggc‑g                                                                                             >  2:1301166/1‑100 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑g                                                                                       <  2:89169/100‑1 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑g                                                                                       <  1:30215/100‑1 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑g                                                                                       <  1:1368226/100‑1 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑g                                                                                       >  1:1319997/1‑100 (MQ=255)
    gCCTATTATAAAAATCATTACTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑g                                                                                       <  1:774495/100‑1 (MQ=255)
     ccTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑gg                                                                                      <  1:1031423/100‑2 (MQ=255)
       tattatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑ggat                                                                                    <  1:855041/100‑4 (MQ=255)
       tattatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑ggat                                                                                    >  1:1203428/1‑97 (MQ=255)
       tattatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑ggat                                                                                    <  1:959758/100‑4 (MQ=255)
        cgatctaaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑ggatt                                                                                   <  2:1095770/94‑5 (MQ=255)
          tataAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGt                                                                                 >  1:1311572/1‑100 (MQ=255)
            taAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTAc                                                                               >  1:352267/1‑100 (MQ=255)
             aaaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACt                                                                              >  2:227204/1‑100 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTg                                                                             <  2:1311582/100‑1 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTg                                                                             >  2:140927/1‑100 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑g‑g‑g‑ggattagatcg                                                                             >  1:1095760/1‑90 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑G‑G‑G‑G‑ATTGTACTGa                                                                            <  1:385691/100‑1 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCCATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTg                                                                             >  2:1357599/1‑100 (MQ=255)
               aaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGa                                                                            <  1:1179366/100‑1 (MQ=255)
                   tccgatctacaAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACttt                                                                        <  2:485482/92‑1 (MQ=255)
                   cATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACttt                                                                        >  1:842148/1‑100 (MQ=255)
                    aTTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTc                                                                       >  2:678703/1‑100 (MQ=255)
                      tcttCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGc                                                                     <  1:782447/97‑1 (MQ=255)
                       gCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGcc                                                                    >  1:181954/1‑100 (MQ=255)
                        cTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGccc                                                                   <  2:311263/100‑1 (MQ=255)
                         tgacaAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑G‑GTG‑G‑ATTGTACTGACTTTCGCCCt                                                                  <  1:1720060/98‑1 (MQ=255)
                         tCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCt                                                                  <  2:1022014/100‑1 (MQ=255)
                         tCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACAAATTCGCGCGG‑G‑G‑G‑ATTGTACTGACTTTCGcaga                                                                  >  2:782452/1‑97 (MQ=255)
                          cacaAAAATATGTTTAGTTTCAAAGTTCTATCATTTACATAATCAATTAAGTTCACGTTTACACACAATCGGGCGG‑G‑G‑G‑AATGGACGGACTTTCGCCCt                                                                  >  1:1253287/1‑99 (MQ=255)
                           acaAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACtttagatcgga                                                                >  1:485479/1‑92 (MQ=255)
                             aaaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACCCACATTCGGGC‑G‑G‑GAG‑ATAGTACTAACTTTCTCCCt                                                                  >  2:1616975/1‑96 (MQ=255)
                             aaaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCttgtt                                                              <  2:769940/100‑1 (MQ=255)
                              aaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCttgttg                                                             <  1:961235/100‑1 (MQ=255)
                              aaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCttgttg                                                             <  1:385248/100‑1 (MQ=255)
                               aaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGc                                                            <  2:361544/100‑1 (MQ=255)
                               aaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGc                                                            <  1:85990/100‑1 (MQ=255)
                                     gtttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACt                                                       >  2:1658540/1‑99 (MQ=255)
                                     gtttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCCGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACt                                                       >  1:1068080/1‑99 (MQ=255)
                                      tttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACttt                                                     <  1:254630/100‑1 (MQ=255)
                                      tttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACttt                                                     >  1:934249/1‑100 (MQ=255)
                                       ttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACtttt                                                    <  1:347816/100‑1 (MQ=255)
                                          gtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACt                                                 <  2:1083430/100‑1 (MQ=255)
                                            ttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTAt                                               >  2:1476536/1‑100 (MQ=255)
                                            ttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑GAG‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTAt                                               >  2:27826/1‑100 (MQ=255)
                                                caGTTCTATCATTTACGTAATCGATTGCTTTCACGTTTACACACATTCGGTCTG‑G‑T‑G‑ATTTTACTTACTTTCGCCCTTGTTGCAAACTTTTACTATCAAt                                           <  1:1099539/100‑1 (MQ=255)
                                                caGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAAt                                           >  1:1552600/1‑100 (MQ=255)
                                                caGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAAt                                           >  1:1555711/1‑100 (MQ=255)
                                                  gTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATcc                                         <  2:778990/100‑1 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCTGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      <  1:600393/100‑1 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      >  2:1523261/1‑100 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      >  1:188357/1‑100 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      <  1:1623481/100‑1 (MQ=255)
                                                      tATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      >  2:556655/1‑99 (MQ=255)
                                                       aTCGTTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGtn                                    <  1:1452947/100‑2 (MQ=255)
                                                       aTCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTg                                    >  2:12830/1‑100 (MQ=255)
                                                        tCATTTACTTAATCTATTGCGTTCACGTTTATACACATTCAATC‑G‑G‑GAG‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGa                                   <  2:109418/100‑1 (MQ=255)
                                                        tCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGa                                   >  2:1051759/1‑100 (MQ=255)
                                                        tCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGa                                   <  1:591660/100‑1 (MQ=255)
                                                         cATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGAt                                  >  1:1639505/1‑100 (MQ=255)
                                                         cATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGAt                                  >  2:1386228/1‑100 (MQ=255)
                                                           tttACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGt                                <  2:1403498/100‑1 (MQ=255)
                                                           tttACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGt                                <  2:250251/100‑1 (MQ=255)
                                                            ttACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTa                               <  1:573661/100‑1 (MQ=255)
                                                             tACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTaa                              >  1:1634052/1‑100 (MQ=255)
                                                               cGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGc                            <  1:138452/100‑1 (MQ=255)
                                                                 tAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCag                          >  1:314000/1‑100 (MQ=255)
                                                                 tAATCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCag                          >  1:118295/1‑100 (MQ=255)
                                                                   aTCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCagag                        >  1:1648442/1‑100 (MQ=255)
                                                                   aTCGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCagag                        <  2:1085452/100‑1 (MQ=255)
                                                                     cGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTa                      >  1:816613/1‑100 (MQ=255)
                                                                     cGATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTa                      <  1:1357589/100‑1 (MQ=255)
                                                                      gATTGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTaa                     >  2:31222/1‑100 (MQ=255)
                                                                       aTTTCGTTCACGTGTACACACATTCGGGC‑G‑G‑G‑GTAGTGTACTGACTTGCGCCCTTGTTTCAAACTTTTACTATCAATGCAGGTGATGTAAGCAGATTAAc                    >  1:534920/1‑100 (MQ=255)
                                                                        ttGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACa                   <  2:876431/100‑1 (MQ=255)
                                                                         tGCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaa                  <  2:842154/100‑1 (MQ=255)
                                                                          gCGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaa                 >  2:568658/1‑100 (MQ=255)
                                                                           cGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaaa                >  1:546575/1‑100 (MQ=255)
                                                                           cGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaaa                >  2:92981/1‑100 (MQ=255)
                                                                           cGTTCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaaa                >  1:1796585/1‑100 (MQ=255)
                                                                              tCACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa             >  1:649697/1‑100 (MQ=255)
                                                                               cACGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa             >  1:1109132/1‑99 (MQ=255)
                                                                               cACGTTTACACACATTCGGGC‑G‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACg           >  2:1034716/1‑100 (MQ=255)
                                                                                 cGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg          <  2:352269/100‑1 (MQ=255)
                                                                                 cGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg          <  1:1487720/100‑1 (MQ=255)
                                                                                 cGTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg          <  2:181954/100‑1 (MQ=255)
                                                                                  gTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACggg         <  1:87326/100‑1 (MQ=255)
                                                                                  gTTTACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACggg         <  2:1214680/100‑1 (MQ=255)
                                                                                     tACACACATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGaa      >  1:1366538/1‑100 (MQ=255)
                                                                                      acacacATTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGATGATGTAAGCAGAGTAACAAAATGACGGTGaaa     >  1:526448/1‑100 (MQ=255)
                                                                                       cacacaTTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAg    >  2:1281482/1‑100 (MQ=255)
                                                                                       cacacaTTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCAAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAg    >  1:530373/1‑100 (MQ=255)
                                                                                         cacaTTCGGGCGG‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAGAt  >  2:631744/1‑100 (MQ=255)
                                                                                                    |                                                                                            
AGGAGCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGC‑G‑G‑G‑G‑ATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAGAT  >  NC_000913/432090‑432277

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