Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,622,834 A→T 5.4% noncoding (81/98 nt) mgrR ← sRNA antisense regulator down regulates YgdQ and EptB, PhoPQ‑regulated, Hfq‑dependent; affects sensitivity to antimicrobial peptides

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8340AT 5.4% 324.2 / 5.3 103noncoding (81/98 nt)mgrRsRNA antisense regulator down regulates YgdQ and EptB, PhoPQ‑regulated, Hfq‑dependent; affects sensitivity to antimicrobial peptides
Reads supporting (aligned to +/- strand):  ref base A (44/53);  new base T (4/2);  total (48/55)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 103 total aligned reads displayed.

TTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGTTAATCACATT  >  NC_000913/1622689‑1622980
                                                                                                                                                  |                                                                                                                                                  
ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAAcc                                                                                                                                                 >  1:863358/1‑149 (MQ=255)
ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAAcc                                                                                                                                                 >  2:1087994/1‑149 (MQ=255)
   tattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACcggc                                                                                                                                              <  2:899616/149‑1 (MQ=255)
      tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                                                           <  1:646595/149‑1 (MQ=255)
          tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAATCCGGCGGTGAAt                                                                                                                                       >  1:1032069/1‑149 (MQ=255)
           ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATATTCCCTTACCTAACCGGCCCCTTGTTACAGTTGTGTACATGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCGGTTTACTGGCGGtttttt                                                                                                                                      <  1:1134283/149‑6 (MQ=255)
           ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                                      <  1:1013013/149‑1 (MQ=255)
           ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                                      <  2:935530/149‑1 (MQ=255)
           ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                                      >  1:925070/1‑149 (MQ=255)
              tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACACGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                                                   <  1:1369423/149‑1 (MQ=255)
               ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                                                  >  1:1524220/1‑149 (MQ=255)
               ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCCCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                                                  <  2:1356298/149‑1 (MQ=255)
                 ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                                                <  2:869713/149‑1 (MQ=255)
                  tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                                               <  1:299817/149‑1 (MQ=255)
                    gCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                                             >  2:1093816/1‑149 (MQ=255)
                     ccGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                                            <  1:533683/149‑1 (MQ=255)
                       gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa                                                                                                                          >  2:1474858/1‑149 (MQ=255)
                            tttCAGCCACAAATGGGATGACTAATGAACGGATATAATCCCGCACCTACCCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                                     <  2:87101/149‑1 (MQ=255)
                            tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                                     <  2:1419763/149‑1 (MQ=255)
                                  cctcaAATGGGATGACAAATGATCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACCGTTGTGTACAAGGGGCCTGATTTGTATGACGGCG‑ATTAAAAACCGTCAGTAAACCGGCTGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                                               <  1:66217/146‑1 (MQ=255)
                                  ccACAAATGGGATGACTAATGAACGGAGATTATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                                               >  2:184307/1‑149 (MQ=255)
                                  caacaAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                                               >  2:192930/3‑149 (MQ=255)
                                    acaAATGGGATGACTAATGAACGGAGTTAATCCCTCACCTACCCGGCCCCTTTTTACAGTTGTGTACATGGGGCCTGATTTTTATGACGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                                             <  1:101584/149‑1 (MQ=255)
                                      aaaTGGGATGACTAATGAACGGAGATAATCCCTCACCT‑ACCGGCCCCTTGTTTCAGTTGTGAACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                                          >  1:549464/1‑149 (MQ=255)
                                        aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                                         <  1:1151937/149‑1 (MQ=255)
                                        aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                                         <  2:1084659/149‑1 (MQ=255)
                                          gggATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTa                                                                                                       <  2:40957/149‑1 (MQ=255)
                                            gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                                     >  2:930940/1‑149 (MQ=255)
                                                  tAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGATaaa                                                                                                >  2:214978/1‑148 (MQ=255)
                                                  tAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa                                                                                               >  1:1422999/1‑149 (MQ=255)
                                                     tGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc                                                                                            <  2:1575387/149‑1 (MQ=255)
                                                       aaCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                                                          >  2:1599177/1‑149 (MQ=255)
                                                           gagaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                                                      >  2:120340/1‑149 (MQ=255)
                                                           gagaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                                                      >  2:558861/1‑149 (MQ=255)
                                                             gaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                                                    <  1:337771/149‑1 (MQ=255)
                                                              aTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTg                                                                                   <  1:1085760/149‑1 (MQ=255)
                                                              aTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTg                                                                                   <  2:315932/149‑1 (MQ=255)
                                                               tAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                                                  >  2:807462/1‑149 (MQ=255)
                                                               tAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                                                  <  1:1559187/149‑1 (MQ=255)
                                                                aaTCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTTCTTTTACGCTAACAGGCATTTTCCTGCa                                                                                 <  1:858043/149‑1 (MQ=255)
                                                                aaTCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCa                                                                                 >  1:569727/1‑149 (MQ=255)
                                                                aaTCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTACTGCa                                                                                 >  2:1588562/1‑149 (MQ=255)
                                                                  tCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                                                >  1:575397/1‑149 (MQ=255)
                                                                  tCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                                               <  2:275760/149‑1 (MQ=255)
                                                                    ccTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa                                                                             >  2:871116/1‑149 (MQ=255)
                                                                    ccTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCGGAATTTTATGCCGGCG‑AAAAAAAACCGCCAGTAAACCGGCGCTGAAGGCATGAATGGATAGATTTGTGTTTTGCATTTTCGCTAACCGGCAGTTTCCTGCACTGa                                                                             >  2:810825/1‑149 (MQ=255)
                                                                     cTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                                            >  1:1204677/1‑149 (MQ=255)
                                                                     cTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                                            <  1:628857/149‑1 (MQ=255)
                                                                      tCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                                           >  2:1120438/1‑149 (MQ=255)
                                                                       cACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa                                                                          <  1:803729/149‑1 (MQ=255)
                                                                        aCCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                                         <  1:317767/149‑1 (MQ=255)
                                                                              ccGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                                                   >  1:777592/1‑149 (MQ=255)
                                                                               cGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCATTAAATCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGt                                                                  <  2:1591112/149‑1 (MQ=255)
                                                                                ggCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGTGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                                 <  1:853775/149‑1 (MQ=255)
                                                                                ggCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                                 <  1:1599177/149‑1 (MQ=255)
                                                                                  ccccTTGTTACAGTTTTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGa                                                               <  2:715244/149‑1 (MQ=255)
                                                                                  ccccTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGa                                                               <  1:1093816/149‑1 (MQ=255)
                                                                                    ccTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGCATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATAGTTGtca                                                             >  1:223449/1‑149 (MQ=255)
                                                                                     cTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                                            >  2:124022/1‑149 (MQ=255)
                                                                                     cTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                                            <  1:1087994/149‑1 (MQ=255)
                                                                                        gTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                                            <  2:1541691/130‑1 (MQ=255)
                                                                                        gTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                                            >  1:1541691/1‑130 (MQ=255)
                                                                                        gTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                                         >  1:603187/1‑149 (MQ=255)
                                                                                         ttACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTATCAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                                        >  2:245967/1‑149 (MQ=255)
                                                                                         ttACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                                        >  1:1423062/1‑149 (MQ=255)
                                                                                             aGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                                    >  2:1331172/1‑149 (MQ=255)
                                                                                               ttGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAAAGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACCGTAGCATCa                                                  >  2:48865/1‑149 (MQ=255)
                                                                                                  tgtACAAGTGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                                               <  2:604934/149‑1 (MQ=255)
                                                                                                    tACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                                             <  2:1422999/149‑1 (MQ=255)
                                                                                                    tACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGAATCAGtttt                                             >  1:1255954/1‑149 (MQ=255)
                                                                                                    tACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGAGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                                             <  2:863358/149‑1 (MQ=255)
                                                                                                     aCAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc                                            >  1:880983/1‑149 (MQ=255)
                                                                                                      cAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                                           <  2:348144/149‑1 (MQ=255)
                                                                                                      cAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                                           >  1:1357073/1‑149 (MQ=255)
                                                                                                      cAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                                           <  2:606502/149‑1 (MQ=255)
                                                                                                            gCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                                     >  1:492240/1‑149 (MQ=255)
                                                                                                            gCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                                     >  2:947131/1‑149 (MQ=255)
                                                                                                             ccTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa                                    <  2:1561152/149‑1 (MQ=255)
                                                                                                             ccTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa                                    <  2:1085057/149‑1 (MQ=255)
                                                                                                              cTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                                            >  1:577819/1‑108 (MQ=255)
                                                                                                              cTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                                            <  2:577819/108‑1 (MQ=255)
                                                                                                               tGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg                                  >  1:133536/1‑149 (MQ=255)
                                                                                                                 aTTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                                <  2:777592/149‑1 (MQ=255)
                                                                                                                  tttttATGACGGCGAAAAAAAAACCGCCAGTAAACAGGCTGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                                <  1:214978/149‑1 (MQ=255)
                                                                                                                   ttttATGACGGCG‑ATAAAAAACCGCCTGTATACCGGAGGTTATTGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTATCAGGCATTTTCCTGCACTGATAACGATTCGTTGACACAGTAGCATCAGTTTTCTCATTGAATGTTaa                              <  1:184307/149‑1 (MQ=255)
                                                                                                                    tttATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaaa                             <  1:947131/149‑1 (MQ=255)
                                                                                                                      tATGACGGCGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc                            <  2:575397/149‑1 (MQ=255)
                                                                                                                            ggCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCtt                     >  1:875361/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaaCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAAc                 >  1:1450159/1‑149 (MQ=255)
                                                                                                                                  gcattcaCCGCCGGTTTACTGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt                >  1:199179/7‑149 (MQ=255)
                                                                                                                                     ttcaCCGCCGGTTTACTGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCgg             >  1:1337329/4‑149 (MQ=255)
                                                                                                                                     ttcaCCGCCGGTTTACTGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCgg             >  1:1400833/4‑149 (MQ=255)
                                                                                                                                     aaatCCGCCAGTAAACCGGCGGGGAATGCTTGCATGGAAAGATTTGTGGTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAGCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCgg             <  2:223338/149‑1 (MQ=255)
                                                                                                                                     aaaaCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCgg             <  1:120340/149‑1 (MQ=255)
                                                                                                                                      aaaCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGt            <  1:1331172/149‑1 (MQ=255)
                                                                                                                                         ccgccgGTTTACTGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGTTaa         >  2:1392858/1‑149 (MQ=255)
                                                                                                                                          cGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGTTAAt        <  2:1470957/149‑1 (MQ=255)
                                                                                                                                          cGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGTTAAt        <  1:377071/149‑1 (MQ=255)
                                                                                                                                            ccAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGTTAATca      >  1:469681/1‑149 (MQ=255)
                                                                                                                                                tAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGTTAATCACAtt  <  2:1497040/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACTCGGTTAATCACATT  >  NC_000913/1622689‑1622980

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