Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,440,923 C→T G36D (GGC→GAC)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,9230CT98.7% 291.4 / NA 76G36D (GGC→GAC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (43/32);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGA  >  NC_000913/3440776‑3441063
                                                                                                                                                   |                                                                                                                                            
aGGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCa                                                                                                                                            >  1:182806/1‑150 (MQ=255)
     gAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGcc                                                                                                                                       >  2:167992/1‑150 (MQ=255)
      aGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCAc                                                                                                                                     <  1:644248/151‑1 (MQ=255)
       gCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACg                                                                                                                                    <  2:97506/151‑1 (MQ=255)
        cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGc                                                                                                                                   >  2:766393/1‑151 (MQ=255)
         aGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGct                                                                                                                                  >  1:136578/1‑151 (MQ=255)
               tCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTaa                                                                                                                              >  1:301053/1‑149 (MQ=255)
               tCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTaa                                                                                                                              >  1:52101/1‑149 (MQ=255)
                  aGGATATCTTCCTGAACGCCTTCTTTGGAGCTGTACTCATGTAGTACAGCATCAAGCTCAACCTCGGTCACCGCGCAACACGGCATCGATGGGAGCAGAATACGCCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGTCt                                                                                                                         >  2:815045/1‑151 (MQ=255)
                   ggATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTc                                                                                                                        <  2:257644/151‑1 (MQ=255)
                   ggATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTc                                                                                                                        <  1:506006/151‑1 (MQ=255)
                      tatCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCa                                                                                                                       >  2:642470/1‑149 (MQ=255)
                           tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGt                                                                                                                  <  2:374923/149‑1 (MQ=255)
                           tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTc                                                                                                                 <  2:162652/150‑1 (MQ=255)
                                aaCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTg                                                                                                           <  2:964626/151‑1 (MQ=255)
                                 aCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTgg                                                                                                          <  2:346091/151‑1 (MQ=255)
                                  cGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                         >  2:347380/1‑151 (MQ=255)
                                    ccTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgt                                                                                                       <  2:339824/151‑1 (MQ=255)
                                     cttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtg                                                                                                      >  1:339000/1‑151 (MQ=255)
                                      ttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtgc                                                                                                     >  1:123267/1‑151 (MQ=255)
                                      ttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtgc                                                                                                     >  1:568008/1‑151 (MQ=255)
                                       tcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtgcg                                                                                                    >  1:38469/1‑151 (MQ=255)
                                       tcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtgcg                                                                                                    >  1:663798/1‑151 (MQ=255)
                                            ggTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGaa                                                                                               <  2:14377/151‑1 (MQ=255)
                                             gTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAAc                                                                                              <  1:30796/151‑1 (MQ=255)
                                                 tGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAc                                                                                          >  1:666038/1‑151 (MQ=255)
                                                    aCTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                                                                                        >  2:983417/1‑150 (MQ=255)
                                                    aCTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTg                                                                                       >  2:166282/1‑151 (MQ=255)
                                                         tGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCg                                                                                   >  1:482828/1‑150 (MQ=255)
                                                          gtagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                 >  1:475795/1‑151 (MQ=255)
                                                          gtagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                 <  2:163110/151‑1 (MQ=255)
                                                            agtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                               >  2:424048/1‑151 (MQ=255)
                                                            agtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                               <  2:568463/151‑1 (MQ=255)
                                                            agtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                               >  1:992852/1‑151 (MQ=255)
                                                             gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATc                                                                              <  2:292430/151‑1 (MQ=255)
                                                                 acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                          >  2:729845/1‑151 (MQ=255)
                                                                  ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                          >  1:286909/1‑150 (MQ=255)
                                                                     tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgc                                                                      >  1:81113/1‑151 (MQ=255)
                                                                      cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                     >  1:760106/1‑151 (MQ=255)
                                                                      cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                     <  2:156055/151‑1 (MQ=255)
                                                                           tcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtt                                                                >  1:334397/1‑151 (MQ=255)
                                                                             aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtttt                                                              <  2:727042/151‑1 (MQ=255)
                                                                             aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtttt                                                              >  2:344093/1‑151 (MQ=255)
                                                                              aCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTa                                                             <  1:308714/151‑1 (MQ=255)
                                                                                cTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtttt                                                              >  1:177803/1‑148 (MQ=255)
                                                                                 tCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGa                                                           >  2:735351/1‑150 (MQ=255)
                                                                                 tCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGa                                                           >  2:654058/1‑150 (MQ=255)
                                                                                  cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaaa                                                         <  2:756166/151‑1 (MQ=255)
                                                                                     tCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAct                                                       >  2:713624/1‑150 (MQ=255)
                                                                                      cACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                      <  1:924758/150‑1 (MQ=255)
                                                                                         cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTc                                                  >  2:942338/1‑151 (MQ=255)
                                                                                          gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                 <  2:768509/151‑1 (MQ=255)
                                                                                                 ccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAcc                                          >  1:807378/1‑151 (MQ=255)
                                                                                                  cGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAccc                                         >  2:234066/1‑151 (MQ=255)
                                                                                                  cGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAccc                                         >  1:991873/1‑151 (MQ=255)
                                                                                                        cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                   <  1:370850/151‑1 (MQ=255)
                                                                                                        cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAt                                    <  1:538296/150‑1 (MQ=255)
                                                                                                           tGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgtg                                >  2:52256/1‑151 (MQ=255)
                                                                                                             agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                              >  2:77576/1‑151 (MQ=255)
                                                                                                              gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTcc                             >  2:753544/1‑151 (MQ=255)
                                                                                                               agCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                            <  1:841262/151‑1 (MQ=255)
                                                                                                                  aGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                         >  1:560075/1‑151 (MQ=255)
                                                                                                                    aaTACGGCGCAGTGCGTTTCCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                       <  1:897543/151‑1 (MQ=255)
                                                                                                                       aCGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTg                     <  1:889593/150‑1 (MQ=255)
                                                                                                                       aCGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTg                     <  2:103603/150‑1 (MQ=255)
                                                                                                                         ggCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTa                  <  1:333452/151‑1 (MQ=255)
                                                                                                                            gcAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACt                >  2:224427/1‑150 (MQ=255)
                                                                                                                            gcAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTa               <  1:167993/151‑1 (MQ=255)
                                                                                                                                    ttACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTAc       >  1:599910/1‑151 (MQ=255)
                                                                                                                                     tACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt      >  2:687990/1‑151 (MQ=255)
                                                                                                                                     tACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt      >  1:538135/1‑151 (MQ=255)
                                                                                                                                      aCCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt     <  2:913751/151‑1 (MQ=255)
                                                                                                                                       cccAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt     <  1:545052/150‑1 (MQ=255)
                                                                                                                                       cccAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg    >  2:21021/1‑151 (MQ=255)
                                                                                                                                         cAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa  <  2:568000/151‑1 (MQ=255)
                                                                                                                                         cAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa  <  2:483231/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
AGGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGA  >  NC_000913/3440776‑3441063

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