Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,630 A→G T428A (ACT→GCT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,6300AG97.3% 230.4 / NA 73T428A (ACT→GCT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (41/30);  minor base C (0/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAG  >  NC_000913/4186484‑4186774
                                                                                                                                                  |                                                                                                                                                
cGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGct                                                                                                                                                <  1:546762/149‑1 (MQ=255)
  ttCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGctct                                                                                                                                              >  1:1278491/1‑149 (MQ=255)
  ttCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTAGCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGctcg                                                                                                                                              >  1:32846/1‑148 (MQ=255)
    cATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGc                                                                                                                                            <  1:518499/149‑1 (MQ=255)
      tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCAc                                                                                                                                          >  1:796529/1‑149 (MQ=255)
      tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCAc                                                                                                                                          >  1:583123/1‑149 (MQ=255)
      tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCAc                                                                                                                                          >  1:1179402/1‑149 (MQ=255)
       cTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGTTGATCTGTGCACCGGCTCTGCAcc                                                                                                                                         <  2:544227/149‑1 (MQ=255)
        tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCg                                                                                                                                        >  1:200446/1‑149 (MQ=255)
           ggCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCgccg                                                                                                                                     >  2:924631/1‑146 (MQ=255)
                 cTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCCCCGGCTCTGCACCGTCTGGGTAt                                                                                                                               >  2:133691/1‑149 (MQ=255)
                   ggAACTGCGTGGTCTTGCTACAACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATcc                                                                                                                             >  2:1256772/1‑149 (MQ=255)
                    gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCGGGCACTGCTGAACCGTGCACCGGCTCTGCACCGTGCGGGTGTCCa                                                                                                                            >  2:723879/1‑149 (MQ=255)
                    gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCa                                                                                                                            >  2:1051684/1‑149 (MQ=255)
                          cttGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCAtt                                                                                                                      <  2:1348326/147‑1 (MQ=255)
                           gTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTATGCCACGCCCGGGTATCCAGGCAgtt                                                                                                                     >  2:1338127/1‑149 (MQ=255)
                              ggcTTGCTACCACCATTAACGCTGCGGAGGGTATTTTTCAGCGCGACGAAGCCGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                                                                                  <  1:183202/147‑1 (MQ=255)
                               tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                                                                                 >  2:1258412/1‑149 (MQ=255)
                               tCTTGCTACCACCATTAAAGCCGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                                                                                 <  1:944817/149‑1 (MQ=255)
                                cTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAAcc                                                                                                                >  1:381272/1‑149 (MQ=255)
                                cTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAAcc                                                                                                                >  2:1086218/1‑149 (MQ=255)
                                 ttacTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCACTTGAACCg                                                                                                               >  2:616063/4‑149 (MQ=255)
                                       ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                                                         >  1:1307509/1‑149 (MQ=255)
                                       ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                                                         >  1:50769/1‑149 (MQ=255)
                                         accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGAt                                                                                                       >  2:1268066/1‑149 (MQ=255)
                                           caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCg                                                                                                     <  1:575338/149‑1 (MQ=255)
                                            aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                                                                                    >  1:1087379/1‑149 (MQ=255)
                                            aTTAAAGCTGCGAAGAAAAGGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                                                                                    <  1:820464/149‑1 (MQ=255)
                                              tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                                                                                                  <  2:583123/149‑1 (MQ=255)
                                                   cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTGGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                                                                             <  1:1125806/149‑1 (MQ=255)
                                                        aagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATGCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTAt                                                                                        >  1:974119/1‑149 (MQ=255)
                                                              aTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCGCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                                                  <  1:976597/149‑1 (MQ=255)
                                                                      gcgcGAAGAAGCGGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGc                                                                          <  2:5569/149‑1 (MQ=255)
                                                                          gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCCGGCAGTTGAACCGGTACTGAGCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGt                                                                      >  1:245273/1‑149 (MQ=255)
                                                                            agaagCTGTCGTTTGGGATACCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGttt                                                                    <  2:280139/149‑1 (MQ=255)
                                                                                 cggTCGTTTGGGATATCCTGGACGAAGTCATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGCTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCg                                                               <  1:822078/147‑1 (MQ=255)
                                                                                  tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCgg                                                              <  1:1051684/149‑1 (MQ=255)
                                                                                  cgctgTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGGACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCgg                                                              <  1:298507/145‑1 (MQ=255)
                                                                                    tCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCa                                                            >  1:998387/1‑149 (MQ=255)
                                                                                           ggATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCGGGTTTGTGCGGCATATAACg                                                     >  2:769296/1‑149 (MQ=255)
                                                                                             atatCCTCGACGAATTTATCCGCTAACACCCGGTACTGCTGAACCGTGCACCGGCTCCGCACCGTCTGGGTATCCAGGCATTTGAACCTGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGcc                                                   <  1:51548/149‑1 (MQ=255)
                                                                                                 ccTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCCGCACCCGCTGGCTTGTGCGGCATATAACGCCGACt                                               >  2:118725/1‑149 (MQ=255)
                                                                                                  cTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACtt                                              >  2:436989/1‑149 (MQ=255)
                                                                                                  cTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACtt                                              <  1:203877/149‑1 (MQ=255)
                                                                                                   tGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGAAAAGCTATCCAGCTGCACCCGCGGGTTTGGGCGGCAGATAACGCCGACTTc                                             >  2:61998/1‑149 (MQ=255)
                                                                                                      aCGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGAt                                          <  2:1278491/149‑1 (MQ=255)
                                                                                                      aCGAAGTCAGCCGCGACCCCCCGGTACTGCTGAACCGTGCACCGGCTCCGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGAt                                          <  2:452830/149‑1 (MQ=255)
                                                                                                       cGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATg                                         >  1:232714/1‑149 (MQ=255)
                                                                                                         aaGTTACCCCCGAACCCCCGGTACGGCGTAACCGGGCACCGCCCCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGt                                       <  1:133691/149‑1 (MQ=255)
                                                                                                          aggtATCCGCGTACACCCGGTACTGCTGAACCATGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGACGGtt                                      <  2:765909/146‑2 (MQ=255)
                                                                                                             tacccGCGAACACCCGGTACTGCGGAACCGGGCACCGGCTCTGCACCGGCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc                                   <  2:381272/146‑1 (MQ=255)
                                                                                                              aTCCGCGAACACCCGGTACTGCTGAACCGCGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCa                                  <  1:644510/149‑1 (MQ=255)
                                                                                                                 cgcgAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGAt                               >  1:1025492/1‑149 (MQ=255)
                                                                                                                 cgcgAACACCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGCCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGGGACGAGAt                               >  2:696487/1‑149 (MQ=255)
                                                                                                                       ccaccGGTTCTGCTGACCCGGGCCCCGCCTCTGCAGCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGGAAATCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGt                         <  1:844565/146‑1 (MQ=255)
                                                                                                                        acCCGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGAGGCATATAACGCCCACTTCGATGGAGACCAGATGGCTGtt                        >  2:21202/1‑149 (MQ=255)
                                                                                                                         cccGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTc                       <  1:1258412/149‑1 (MQ=255)
                                                                                                                          ccGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGGCTTCGATGGTGACCAGATGGCTGTTCa                      >  1:1088199/1‑149 (MQ=255)
                                                                                                                          ccGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGCTGGTGACCAGATGGCTGTTCa                      >  1:595973/1‑149 (MQ=255)
                                                                                                                          ccGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGGTGGCTGTTCa                      >  1:283745/1‑149 (MQ=255)
                                                                                                                          ccGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCa                      >  2:580206/1‑149 (MQ=255)
                                                                                                                          ccGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGGGACCAGATGGCTGTTCa                      >  1:1088048/1‑149 (MQ=255)
                                                                                                                          ccGGTACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCAGGTTTGTGCGGCGTATAACGCCGGCTTCGATGGTGACCCGATGGCTGTTCa                      >  2:580524/1‑149 (MQ=255)
                                                                                                                              tACTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTa                  <  2:50769/149‑1 (MQ=255)
                                                                                                                               aCTGCTGAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTAc                 >  1:320946/1‑149 (MQ=255)
                                                                                                                                 tgctgAACCGTGCACTGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGCAGGCAAAGCTATCCAGCTGCAACCGCTGGTTTGCGCGGCCTGTAACGCCGACTTCGATGGTGACAACATGGCAGTTCACTGACCg               >  2:141524/1‑149 (MQ=255)
                                                                                                                                  gctgAACCGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTAccac              <  2:1179402/149‑3 (MQ=255)
                                                                                                                                  gctgAACCGTGCACCGGCACTGCACCGTCTGGGTATCCAGGCATTTGAACAGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGc              >  2:877796/1‑149 (MQ=255)
                                                                                                                                        ccGTGCACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGc        <  1:1133214/149‑1 (MQ=255)
                                                                                                                                          gTGCACCGGCTCTGCACCGTCTTGTTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTg      <  1:1283320/149‑1 (MQ=255)
                                                                                                                                             ccccGGCTCCGCACCGTCTGGGTATGCAGGCATTGGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGaa   <  2:99717/147‑1 (MQ=255)
                                                                                                                                             cACCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGaa   >  2:556480/1‑149 (MQ=255)
                                                                                                                                              aCCGGCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAg  <  1:138016/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAG  >  NC_000913/4186484‑4186774

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