Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,214,853 C→A 45.5% R603S (CGT→AGT)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,8530CA45.5% 41.3 / 154.7 110R603S (CGT→AGT) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base A (24/26);  total (55/55)
Fisher's exact test for biased strand distribution p-value = 8.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-02

Only 100 of 111 total aligned reads displayed.

TGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGC  >  NC_000913/3214708‑3215000
                                                                                                                                                 |                                                                                                                                                   
tGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGtt                                                                                                                                                  >  2:404015/1‑149 (MQ=255)
 gAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCATTc                                                                                                                                                 >  2:260029/1‑149 (MQ=255)
 gAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTc                                                                                                                                                 >  7:3025/1‑149 (MQ=255)
  aaGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCt                                                                                                                                                <  8:359976/149‑1 (MQ=255)
    gCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTg                                                                                                                                               >  5:66562/1‑148 (MQ=255)
     cAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGaa                                                                                                                                             <  1:404015/149‑1 (MQ=255)
      aaaaGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAg                                                                                                                                            >  3:274342/1‑149 (MQ=255)
      aaaaGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAg                                                                                                                                            >  3:280404/1‑149 (MQ=255)
        aaGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGtg                                                                                                                                          >  4:179087/1‑149 (MQ=255)
         aGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGtgc                                                                                                                                         >  3:102055/1‑149 (MQ=255)
             cTGCGTATGCGTTTCGGTTTCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGtgctgc                                                                                                                                      <  1:39183/148‑1 (MQ=255)
              tGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCg                                                                                                                                     >  1:406142/1‑148 (MQ=255)
                cGTATGCGTTTCGGTATCGATATGAACACCGACTTCACGCTGGAAGAAGTCGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGAGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTa                                                                                                                                   <  5:254670/148‑1 (MQ=255)
                cGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCg                                                                                                                                     <  3:40440/146‑1 (MQ=255)
                cGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTa                                                                                                                                   <  3:179087/148‑1 (MQ=255)
                        tttCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGaa                                                                                                                                             >  3:350372/1‑130 (MQ=255)
                        tttCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGaa                                                                                                                                             <  4:350372/130‑1 (MQ=255)
                             gTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGAt                                                                                                                     >  4:276673/1‑149 (MQ=255)
                              tATCGATATGAACACCGACTACCCGCTGGAAGGAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGAt                                                                                                                     >  6:485691/1‑148 (MQ=255)
                               aTCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGAtt                                                                                                                    >  3:281814/1‑148 (MQ=255)
                                  gATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATc                                                                                                                <  8:3025/149‑1 (MQ=255)
                                  gATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATc                                                                                                                <  5:88219/149‑1 (MQ=255)
                                   atatGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCg                                                                                                               >  1:365428/1‑149 (MQ=255)
                                    tatGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCgg                                                                                                              >  2:327632/1‑149 (MQ=255)
                                    tatGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCgg                                                                                                              <  6:363147/149‑1 (MQ=255)
                                      tGAACACCTACTACACGCTGGAAGAAGTGGGTAAACAGTTCTACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGGAGCTTCCTGGACGATTAATCGGTa                                                                                                            <  1:142989/149‑1 (MQ=255)
                                      tGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTa                                                                                                            <  4:102055/149‑1 (MQ=255)
                                      tGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTATTCGGTa                                                                                                            <  6:254499/149‑1 (MQ=255)
                                       gAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAg                                                                                                           >  6:48075/1‑149 (MQ=255)
                                           acCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTa                                                                                                                                   >  7:167434/1‑121 (MQ=255)
                                           acCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTa                                                                                                                                   <  8:167434/121‑1 (MQ=255)
                                             cGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGa                                                                                                     >  7:220720/1‑149 (MQ=255)
                                              gACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGAt                                                                                                    <  7:282417/149‑1 (MQ=255)
                                                cTACACGCTGGAAGAAGTGTGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCa                                                                                                  <  2:85118/149‑1 (MQ=255)
                                                   cacGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAg                                                                                                 <  4:281814/147‑1 (MQ=255)
                                                   cacGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGc                                                                                               <  1:327632/149‑1 (MQ=255)
                                                          gaagaaGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc                                                                                        >  3:344900/1‑149 (MQ=255)
                                                          gaagaaGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCGTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc                                                                                        <  2:406142/149‑1 (MQ=255)
                                                          gaagaaGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCCTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc                                                                                        >  6:25154/1‑149 (MQ=255)
                                                           aagaagTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgcc                                                                                       >  5:21784/1‑149 (MQ=255)
                                                           aagaagTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgcc                                                                                       <  4:17896/149‑1 (MQ=255)
                                                                gTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCAcc                                                                                  >  5:90268/1‑149 (MQ=255)
                                                                  gggTAAACAGTTCGACGTTACCCTCGAACGTATCCTTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCg                                                                                >  8:280121/1‑149 (MQ=255)
                                                                  gggTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa                                                                                    >  7:465161/1‑145 (MQ=255)
                                                                  gggTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa                                                                                    <  8:465161/145‑1 (MQ=255)
                                                                    gTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGc                                                                              <  5:48075/149‑1 (MQ=255)
                                                                     tAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGaa                                                                                                                                             >  4:143856/1‑85 (MQ=255)
                                                                     tAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGaa                                                                                                                                             <  3:143856/85‑1 (MQ=255)
                                                                     tAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCa                                                                             <  4:424417/149‑1 (MQ=255)
                                                                         cAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAg                                                                         <  3:134590/149‑1 (MQ=255)
                                                                         cAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAg                                                                         <  5:52915/149‑1 (MQ=255)
                                                                          aGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAgg                                                                        >  3:47216/1‑149 (MQ=255)
                                                                           gTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGc                                                                                               >  1:69594/1‑125 (MQ=255)
                                                                            ttCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGcc                                                                      <  4:47216/149‑1 (MQ=255)
                                                                               gACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc                                                                   >  3:448425/1‑149 (MQ=255)
                                                                                  gTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTg                                                                 <  8:121359/148‑1 (MQ=255)
                                                                                  gTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGc                                                                <  4:280404/149‑1 (MQ=255)
                                                                                  gTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGc                                                                <  4:274342/149‑1 (MQ=255)
                                                                                       ccGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaaa                                                           >  7:139934/1‑149 (MQ=255)
                                                                                       ccGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCCTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaaa                                                           >  6:467275/1‑149 (MQ=255)
                                                                                       ccGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaaa                                                           >  7:200737/1‑149 (MQ=255)
                                                                                          cgAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGc                                                        <  8:139934/149‑1 (MQ=255)
                                                                                                 aTCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt                                                 >  8:372302/1‑149 (MQ=255)
                                                                                                 aTCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt                                                 >  2:390926/1‑149 (MQ=255)
                                                                                                  tCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCAc                                                                            <  5:84075/121‑1 (MQ=255)
                                                                                                  tCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCAc                                                                            >  6:84075/1‑121 (MQ=255)
                                                                                                   ccGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtttt                                                <  8:200737/148‑1 (MQ=255)
                                                                                                    cGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCATTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCg                                              <  1:260029/149‑1 (MQ=255)
                                                                                                    cGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCg                                              <  5:485691/149‑1 (MQ=255)
                                                                                                     gTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCa                                                                                                  >  3:201139/1‑96 (MQ=255)
                                                                                                     gTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCa                                                                                                  <  4:201139/96‑1 (MQ=255)
                                                                                                     gTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCgg                                             >  8:144471/1‑149 (MQ=255)
                                                                                                      tCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCCTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGt                                            <  5:25154/149‑1 (MQ=255)
                                                                                                        aGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTgg                                          >  8:127248/1‑149 (MQ=255)
                                                                                                        aGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTg                                           <  7:280121/148‑1 (MQ=255)
                                                                                                           tCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGt                                       >  1:308164/1‑149 (MQ=255)
                                                                                                           tCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGtt                                                                                            <  4:258537/96‑1 (MQ=255)
                                                                                                           tCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGtt                                                                                            >  3:258537/1‑96 (MQ=255)
                                                                                                              aagcgaagGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttt                                    <  4:448425/149‑1 (MQ=255)
                                                                                                              aagcgaagGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttt                                    <  6:90268/149‑1 (MQ=255)
                                                                                                                 cgaagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTa                                 <  2:308164/149‑1 (MQ=255)
                                                                                                                  gaagGCGCTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTAt                                >  8:408965/1‑149 (MQ=255)
                                                                                                                      gcgcTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGcc                            <  1:390926/149‑1 (MQ=255)
                                                                                                                        gcTGCGCAAACTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATGAATCGGTAGGCCGGATCCGGCGTTACGCCGCACCCGGCCCTAGGCCCTCTGCACAAACGCCACCGTTTCGGTGGCGTTTTTTATCGCCCa                          >  5:321000/1‑149 (MQ=255)
                                                                                                                            cgcAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTc                                                                                                                              <  2:177288/45‑1 (MQ=255)
                                                                                                                              cAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAc                     >  7:262249/1‑148 (MQ=255)
                                                                                                                                aaCTGCGTCACCCGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                  <  7:502832/149‑1 (MQ=255)
                                                                                                                                      gTCACCCGAGCCCTTCTGAAGTGCAGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGCATCAGGCGTTACGCCGCACCCCGCACTAGGCCCTCTGCACAAACCCCACCTTTACGGTGGCGTTTATTATCACCCACGCACAACCAgcgc            >  6:295708/1‑149 (MQ=255)
                                                                                                                                       tCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGcc           >  5:88521/1‑149 (MQ=255)
                                                                                                                                       tCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGcc           >  7:57921/1‑149 (MQ=255)
                                                                                                                                       tCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGcc           >  7:50110/1‑149 (MQ=255)
                                                                                                                                       tCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGcc           >  3:67796/1‑149 (MQ=255)
                                                                                                                                         aCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg         <  2:365428/149‑1 (MQ=255)
                                                                                                                                          cccGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTgg        <  2:75629/149‑1 (MQ=255)
                                                                                                                                           ccGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTa                                 <  7:467502/123‑1 (MQ=255)
                                                                                                                                           ccGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTa                                 >  8:467502/1‑123 (MQ=255)
                                                                                                                                           ccGAGCAGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGt       >  7:313759/1‑149 (MQ=255)
                                                                                                                                            cGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTc      <  4:344900/149‑1 (MQ=255)
                                                                                                                                              aGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCa    <  4:67796/149‑1 (MQ=255)
                                                                                                                                                ccGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGc  >  2:228293/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGC  >  NC_000913/3214708‑3215000

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: