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

AGGAGCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAGAT  >  NC_000913/432090‑432277
                                                                                                    |                                                                                            
aggagCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTcgggc                                                                                               <  1:319065/100‑1 (MQ=255)
 ggagCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTcgggcg                                                                                             >  2:1301166/1‑100 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggg                                                                                       <  2:89169/100‑1 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggg                                                                                       <  1:30215/100‑1 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggg                                                                                       <  1:1368226/100‑1 (MQ=255)
    gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggg                                                                                       >  1:1319997/1‑100 (MQ=255)
    gCCTATTATAAAAATCATTACTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggg                                                                                       <  1:774495/100‑1 (MQ=255)
     ccTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCggggg                                                                                      <  1:1031423/100‑2 (MQ=255)
       tattatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggggat                                                                                    <  1:855041/100‑4 (MQ=255)
       tattatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggggat                                                                                    >  1:1203428/1‑97 (MQ=255)
       tattatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggggat                                                                                    <  1:959758/100‑4 (MQ=255)
        cgatctaaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggggatt                                                                                   <  2:1095770/94‑5 (MQ=255)
          tataAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGt                                                                                 >  1:1311572/1‑100 (MQ=255)
            taAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTAc                                                                               >  1:352267/1‑100 (MQ=255)
             aaaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACt                                                                              >  2:227204/1‑100 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTg                                                                             <  2:1311582/100‑1 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTg                                                                             >  2:140927/1‑100 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCgggggattagatcg                                                                             >  1:1095760/1‑90 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGATTGTACTGa                                                                            <  1:385691/100‑1 (MQ=255)
              aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCCATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTg                                                                             >  2:1357599/1‑100 (MQ=255)
               aaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGa                                                                            <  1:1179366/100‑1 (MQ=255)
                   tccgatctacaAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACttt                                                                        <  2:485482/92‑1 (MQ=255)
                   cATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACttt                                                                        >  1:842148/1‑100 (MQ=255)
                    aTTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTc                                                                       >  2:678703/1‑100 (MQ=255)
                      tcttCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGc                                                                     <  1:782447/97‑1 (MQ=255)
                       gCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGcc                                                                    >  1:181954/1‑100 (MQ=255)
                        cTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGccc                                                                   <  2:311263/100‑1 (MQ=255)
                         tgacaAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGTGGATTGTACTGACTTTCGCCCt                                                                  <  1:1720060/98‑1 (MQ=255)
                         tCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCt                                                                  <  2:1022014/100‑1 (MQ=255)
                         tCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACAAATTCGCGCGGGGGATTGTACTGACTTTCGcaga                                                                  >  2:782452/1‑97 (MQ=255)
                          cacaAAAATATGTTTAGTTTCAAAGTTCTATCATTTACATAATCAATTAAGTTCACGTTTACACACAATCGGGCGGGGGAATGGACGGACTTTCGCCCt                                                                  >  1:1253287/1‑99 (MQ=255)
                           acaAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACtttagatcgga                                                                >  1:485479/1‑92 (MQ=255)
                             aaaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACCCACATTCGGGCGGGAGATAGTACTAACTTTCTCCCt                                                                  >  2:1616975/1‑96 (MQ=255)
                             aaaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCttgtt                                                              <  2:769940/100‑1 (MQ=255)
                              aaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCttgttg                                                             <  1:961235/100‑1 (MQ=255)
                              aaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCttgttg                                                             <  1:385248/100‑1 (MQ=255)
                               aaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGc                                                            <  2:361544/100‑1 (MQ=255)
                               aaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGc                                                            <  1:85990/100‑1 (MQ=255)
                                     gtttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACt                                                       >  2:1658540/1‑99 (MQ=255)
                                     gtttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCCGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACt                                                       >  1:1068080/1‑99 (MQ=255)
                                      tttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACttt                                                     <  1:254630/100‑1 (MQ=255)
                                      tttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACttt                                                     >  1:934249/1‑100 (MQ=255)
                                       ttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACtttt                                                    <  1:347816/100‑1 (MQ=255)
                                          gtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACt                                                 <  2:1083430/100‑1 (MQ=255)
                                            ttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTAt                                               >  2:1476536/1‑100 (MQ=255)
                                            ttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGAGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTAt                                               >  2:27826/1‑100 (MQ=255)
                                                caGTTCTATCATTTACGTAATCGATTGCTTTCACGTTTACACACATTCGGTCTGGTGATTTTACTTACTTTCGCCCTTGTTGCAAACTTTTACTATCAAt                                           <  1:1099539/100‑1 (MQ=255)
                                                caGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAAt                                           >  1:1552600/1‑100 (MQ=255)
                                                caGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAAt                                           >  1:1555711/1‑100 (MQ=255)
                                                  gTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATcc                                         <  2:778990/100‑1 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCTGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      <  1:600393/100‑1 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      >  2:1523261/1‑100 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      >  1:188357/1‑100 (MQ=255)
                                                     cTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      <  1:1623481/100‑1 (MQ=255)
                                                      tATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                                      >  2:556655/1‑99 (MQ=255)
                                                       aTCGTTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGtn                                    <  1:1452947/100‑2 (MQ=255)
                                                       aTCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTg                                    >  2:12830/1‑100 (MQ=255)
                                                        tCATTTACTTAATCTATTGCGTTCACGTTTATACACATTCAATCGGGAGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGa                                   <  2:109418/100‑1 (MQ=255)
                                                        tCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGa                                   >  2:1051759/1‑100 (MQ=255)
                                                        tCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGa                                   <  1:591660/100‑1 (MQ=255)
                                                         cATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGAt                                  >  1:1639505/1‑100 (MQ=255)
                                                         cATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGAt                                  >  2:1386228/1‑100 (MQ=255)
                                                           tttACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGt                                <  2:1403498/100‑1 (MQ=255)
                                                           tttACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGt                                <  2:250251/100‑1 (MQ=255)
                                                            ttACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTa                               <  1:573661/100‑1 (MQ=255)
                                                             tACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTaa                              >  1:1634052/1‑100 (MQ=255)
                                                               cGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGc                            <  1:138452/100‑1 (MQ=255)
                                                                 tAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCag                          >  1:314000/1‑100 (MQ=255)
                                                                 tAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCag                          >  1:118295/1‑100 (MQ=255)
                                                                   aTCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCagag                        >  1:1648442/1‑100 (MQ=255)
                                                                   aTCGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCagag                        <  2:1085452/100‑1 (MQ=255)
                                                                     cGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTa                      >  1:816613/1‑100 (MQ=255)
                                                                     cGATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTa                      <  1:1357589/100‑1 (MQ=255)
                                                                      gATTGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTaa                     >  2:31222/1‑100 (MQ=255)
                                                                       aTTTCGTTCACGTGTACACACATTCGGGCGGGGTAGTGTACTGACTTGCGCCCTTGTTTCAAACTTTTACTATCAATGCAGGTGATGTAAGCAGATTAAc                    >  1:534920/1‑100 (MQ=255)
                                                                        ttGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACa                   <  2:876431/100‑1 (MQ=255)
                                                                         tGCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaa                  <  2:842154/100‑1 (MQ=255)
                                                                          gCGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaa                 >  2:568658/1‑100 (MQ=255)
                                                                           cGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaaa                >  1:546575/1‑100 (MQ=255)
                                                                           cGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaaa                >  2:92981/1‑100 (MQ=255)
                                                                           cGTTCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaaa                >  1:1796585/1‑100 (MQ=255)
                                                                              tCACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa             >  1:649697/1‑100 (MQ=255)
                                                                               cACGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa             >  1:1109132/1‑99 (MQ=255)
                                                                               cACGTTTACACACATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACg           >  2:1034716/1‑100 (MQ=255)
                                                                                 cGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg          <  2:352269/100‑1 (MQ=255)
                                                                                 cGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg          <  1:1487720/100‑1 (MQ=255)
                                                                                 cGTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg          <  2:181954/100‑1 (MQ=255)
                                                                                  gTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACggg         <  1:87326/100‑1 (MQ=255)
                                                                                  gTTTACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACggg         <  2:1214680/100‑1 (MQ=255)
                                                                                     tACACACATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGaa      >  1:1366538/1‑100 (MQ=255)
                                                                                      acacacATTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGATGATGTAAGCAGAGTAACAAAATGACGGTGaaa     >  1:526448/1‑100 (MQ=255)
                                                                                       cacacaTTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAg    >  2:1281482/1‑100 (MQ=255)
                                                                                       cacacaTTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCAAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAg    >  1:530373/1‑100 (MQ=255)
                                                                                         cacaTTCGGGCGGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAGAt  >  2:631744/1‑100 (MQ=255)
                                                                                                    |                                                                                            
AGGAGCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGGGAAAGAT  >  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: