Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,056,123 A→C 6.1% noncoding (141/183 nt) ssrS → 6S RNA inhibits RNA polymerase promoter binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,056,1230AC 6.1% 424.6 / 16.4 127noncoding (141/183 nt)ssrS6S RNA inhibits RNA polymerase promoter binding
Reads supporting (aligned to +/- strand):  ref base A (67/52);  new base C (6/2);  total (73/54)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 149 total aligned reads displayed.

AGACAAAATTTCTCTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGC  >  NC_000913/3055976‑3056271
                                                                                                                                                   |                                                                                                                                                    
aGACAAAATTTCTCTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaa                                                                                                                                                     <  2:993275/149‑1 (MQ=255)
 gACAAAATTTCTCTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaag                                                                                                                                                    <  4:10844/149‑1 (MQ=255)
    aaaaTTTCTCTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggt                                                                                                                                                 >  4:92648/1‑149 (MQ=255)
    aaaaTTTCTCTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggt                                                                                                                                                 >  1:449563/1‑149 (MQ=255)
          tctcTGAGATGTTCGCACGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggttcaagg                                                                                                                                           <  3:418432/149‑1 (MQ=255)
          tctcTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGTGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggttcaagg                                                                                                                                           <  2:457596/149‑1 (MQ=255)
          tctcTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggttcaagg                                                                                                                                           <  3:1000648/149‑1 (MQ=255)
          tctcTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggttcaagg                                                                                                                                           <  1:1025994/149‑1 (MQ=255)
           ctctGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggttcaaggg                                                                                                                                          <  3:493239/149‑1 (MQ=255)
            tctGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggttcaagggt                                                                                                                                         >  4:197556/1‑149 (MQ=255)
             ctGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACcaagggttcaagggtt                                                                                                                                        <  3:1002777/149‑1 (MQ=255)
              tGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTa                                                                                                                                       >  3:141826/1‑149 (MQ=255)
              tGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTa                                                                                                                                       >  2:861166/1‑149 (MQ=255)
               gagaTGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTAc                                                                                                                                      <  2:482121/149‑1 (MQ=255)
                agaTGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACa                                                                                                                                     >  4:426764/1‑149 (MQ=255)
                         cAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTgcggc                                                                                                                            <  3:197556/149‑1 (MQ=255)
                          aaGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTgcggcg                                                                                                                           >  1:1053978/1‑149 (MQ=255)
                          aaGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTgcggcg                                                                                                                           >  1:568425/1‑149 (MQ=255)
                           aGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGcggcgg                                                                                                                          <  3:846781/149‑1 (MQ=255)
                            gCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCggcggc                                                                                                                         >  3:277673/1‑149 (MQ=255)
                                  cAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGTGGCATCTCg                                                                                                                   >  3:493610/1‑149 (MQ=255)
                                    gTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGga                                                                                                                 <  4:575513/149‑1 (MQ=255)
                                    gTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGga                                                                                                                 <  2:449563/149‑1 (MQ=255)
                                    gTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGga                                                                                                                 <  2:803216/149‑1 (MQ=255)
                                       cccTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGGACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGAt                                                                                                              >  4:692502/1‑149 (MQ=255)
                                       cccTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGAt                                                                                                              >  1:447206/1‑149 (MQ=255)
                                       cccTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGAt                                                                                                              >  1:522292/1‑149 (MQ=255)
                                         cTGAGCCGATATTTCATACCACCAGAATGTGGCGCTCCGTGGTTGGTCAGCATGCTCTGTCCGTCCGAGAAGCCTTAAAACTGCGAGGACACATTCACCGTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTc                                                                                                            <  2:522292/149‑1 (MQ=255)
                                          tGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTcc                                                                                                           >  2:609366/1‑149 (MQ=255)
                                           gAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTccc                                                                                                          <  2:302784/149‑1 (MQ=255)
                                            aGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCct                                                                                                         <  4:969772/149‑1 (MQ=255)
                                                   atTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGTGGCATCTCGGAGATTCCCTTCTTATc                                                                                                  <  4:493610/149‑1 (MQ=255)
                                                        aTACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGAtt                                                                                                             >  4:599492/1‑133 (MQ=255)
                                                            cacaAGAATGTGGCGCTCCGCGGTAGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAg                                                                                         >  1:906776/1‑149 (MQ=255)
                                                              caAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGcc                                                                                       <  1:903943/149‑1 (MQ=255)
                                                               aaGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCa                                                                                      >  4:963369/1‑149 (MQ=255)
                                                               aaGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCa                                                                                      >  3:563671/1‑149 (MQ=255)
                                                                        gcgcTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATc                                                                                                  <  3:813028/128‑1 (MQ=255)
                                                                        gcgcTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAAc                                                                             <  4:563671/149‑1 (MQ=255)
                                                                        gcgcTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAAc                                                                             <  4:407573/149‑1 (MQ=255)
                                                                        gcgcTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCACTTCTTATCTGGCACCAGCCATGACGCAAc                                                                             >  4:569277/1‑149 (MQ=255)
                                                                                cgGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGACACCAGCCATGACGCAACTACCAGaa                                                                     <  2:632247/149‑1 (MQ=255)
                                                                                  gTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACt                                                                   >  4:578177/1‑149 (MQ=255)
                                                                                     ggTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTccc                                                                <  3:231386/149‑1 (MQ=255)
                                                                                       tGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCAc                                                              <  3:92648/149‑1 (MQ=255)
                                                                                          gCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGa                                                           <  4:141826/149‑1 (MQ=255)
                                                                                          gCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGa                                                           >  1:740272/1‑149 (MQ=255)
                                                                                             tGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACAt                                                        >  2:223657/1‑149 (MQ=255)
                                                                                                tCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTAGAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTAt                                                     >  4:74252/1‑149 (MQ=255)
                                                                                                 cGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCAGTGACATTATc                                                    <  4:863139/149‑1 (MQ=255)
                                                                                                  ggTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATcc                                                   <  4:277673/149‑1 (MQ=255)
                                                                                                       gtccgAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACa                                              >  1:288004/1‑149 (MQ=255)
                                                                                                        tccgAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCAc                                                              <  2:723511/132‑1 (MQ=255)
                                                                                                        tccgAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCAc                                                              >  1:723511/1‑132 (MQ=255)
                                                                                                        tccgAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACaa                                             <  3:1025527/149‑1 (MQ=255)
                                                                                                        tccgAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACaa                                             >  2:852988/1‑149 (MQ=255)
                                                                                                           gagaAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGa                                                                                   >  3:753788/1‑108 (MQ=255)
                                                                                                           gagaAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACGAGaa                                          >  1:214007/1‑149 (MQ=255)
                                                                                                                  cTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCcttc                                                                                                       <  3:123876/81‑1 (MQ=255)
                                                                                                                   ttAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATACGCaa                                  >  4:892073/1‑149 (MQ=255)
                                                                                                                          tGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGAtt                           <  1:223657/149‑1 (MQ=255)
                                                                                                                            cgacgaCACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCg                         <  1:852988/149‑1 (MQ=255)
                                                                                                                             gacgacACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCgg                        >  4:822628/1‑149 (MQ=255)
                                                                                                                                gacACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCATCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGcca                     <  4:535864/149‑3 (MQ=255)
                                                                                                                                 acacATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAAc                    >  3:869532/1‑149 (MQ=255)
                                                                                                                                  cacaTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCa                                                               >  4:113916/1‑105 (MQ=255)
                                                                                                                                   acaTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATcc                                      <  1:851731/129‑1 (MQ=255)
                                                                                                                                   acaTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAAcgt                  >  2:820076/1‑149 (MQ=255)
                                                                                                                                     aTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATc                                                                                                  <  3:386277/67‑1 (MQ=255)
                                                                                                                                     aTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATc                                                                                                  >  4:386277/1‑67 (MQ=255)
                                                                                                                                     aTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAAcgtcg                >  3:105833/1‑149 (MQ=255)
                                                                                                                                     aTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAAcgtcg                >  4:705596/1‑149 (MQ=255)
                                                                                                                                     aTTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAAcgtcg                >  4:386870/1‑149 (MQ=255)
                                                                                                                                       tCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAAcgtcgtc              <  3:616490/149‑1 (MQ=255)
                                                                                                                                        cACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAAcgtcgtcg             >  4:758105/1‑149 (MQ=255)
                                                                                                                                          ccTTGAACCCTTGGTTCAAGGGTTACAGCCTGCGGCGGTATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTg           >  3:953813/1‑149 (MQ=255)
                                                                                                                                          ccTTGAACCCTTGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAAc                                                                             >  4:176830/1‑83 (MQ=255)
                                                                                                                                          ccTTGAACCCTTGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTg           >  1:635582/1‑149 (MQ=255)
                                                                                                                                          ccTTGAACCCTTGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTg           >  3:814816/1‑149 (MQ=255)
                                                                                                                                          ccTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTAt                                                     <  4:271208/107‑1 (MQ=255)
                                                                                                                                          ccTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTAt                                                     >  3:271208/1‑107 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACCGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCa                                                                                      <  4:76435/73‑1 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCa                                                                                      >  3:76435/1‑73 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATcc                                      <  3:631068/121‑1 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATcc                                      >  4:631068/1‑121 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGc          <  2:252331/149‑1 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGc          <  1:475062/149‑1 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGc          <  3:705596/149‑1 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGc          >  4:923321/1‑149 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGc          >  2:300026/1‑149 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGc          >  3:426165/1‑149 (MQ=255)
                                                                                                                                           cTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGc          <  1:300026/149‑1 (MQ=255)
                                                                                                                                            ttGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCg         >  4:680222/1‑149 (MQ=255)
                                                                                                                                               aaCCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTa      <  3:578177/149‑1 (MQ=255)
                                                                                                                                                aCCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTaa     >  2:975171/1‑149 (MQ=255)
                                                                                                                                                aCCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTaa     >  2:653643/1‑149 (MQ=255)
                                                                                                                                                 ccAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAAc    <  4:880646/149‑1 (MQ=255)
                                                                                                                                                 ccAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAAc    >  4:128654/1‑149 (MQ=255)
                                                                                                                                                  caagggttcaagggttACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACg   >  3:704592/1‑149 (MQ=255)
                                                                                                                                                   aagggttcaagggttACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGc  >  1:108822/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGACAAAATTTCTCTGAGATGTTCGCAAGCGGGCCAGTCCCCTGAGCCGATATTTCATACCACAAGAATGTGGCGCTCCGCGGTTGGTGAGCATGCTCGGTCCGTCCGAGAAGCCTTAAAACTGCGACGACACATTCACCTTGAACCAAGGGTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGC  >  NC_000913/3055976‑3056271

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: