Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,722,301 A→C 5.2% intergenic (+196/+27) sokA → / ← glyS pseudogene, sRNA antisense regulator blocking mokA/hokA translation/glycine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,722,3010AC 5.2% 2369.9 / 96.1 678intergenic (+196/+27)sokA/glySpseudogene, sRNA antisense regulator blocking mokA/hokA translation/glycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (331/310);  new base C (22/14);  total (353/325)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

Only 100 of 685 total aligned reads displayed.

CCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGAAA  >  NC_000913/3722153‑3722448
                                                                                                                                                    |                                                                                                                                                   
ccGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGa                                                                                                                                                     >  2:327219/1‑149 (MQ=255)
  gTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATg                                                                                                                                                   <  2:846970/149‑1 (MQ=255)
     gCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCa                                                                                                                                                >  2:1393913/1‑149 (MQ=255)
     gCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCa                                                                                                                                                <  3:886213/149‑1 (MQ=255)
      cTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAg                                                                                                                                               >  4:1513856/1‑149 (MQ=255)
       tATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAgg                                                                                                                                              >  4:678568/1‑149 (MQ=255)
         ttCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAGAAGCCTGCCAGATGGCAGGCt                                                                                                                                            >  1:1590423/1‑149 (MQ=255)
           cTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAt                                                                                                                                          >  3:485651/1‑149 (MQ=255)
            tAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCtttt                                                                                                                                         <  1:183539/149‑1 (MQ=255)
            tAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                                                                         <  1:398723/149‑1 (MQ=255)
            tAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                                                                         <  1:918221/149‑1 (MQ=255)
            tAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                                                                         <  3:140935/149‑1 (MQ=255)
             aGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTtttt                                                                                                                                        <  4:4468/149‑1 (MQ=255)
             aGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAttt                                                                                                                                        >  1:859906/1‑149 (MQ=255)
               cTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                                                                      <  1:1179476/149‑1 (MQ=255)
                  aataaGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGAt                                                                                                                                                    <  1:1369674/132‑1 (MQ=255)
                  aataaGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaataa                                                                                                                                   <  3:91242/149‑1 (MQ=255)
                       gtttgtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                                                                      >  3:264085/1‑141 (MQ=255)
                        tttgtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATgg                                                                                                                                                  <  1:553361/128‑1 (MQ=255)
                        tttgtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGt                                                                                                                             >  1:254499/1‑149 (MQ=255)
                          tgtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGtta                                                                                                                           <  3:77746/149‑1 (MQ=255)
                            tttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatt                                                                                                                         >  1:1224093/1‑149 (MQ=255)
                               cttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCt                                                                                                                                            <  1:947502/127‑1 (MQ=255)
                                   tcCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCATAAGATGATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCaa                                                                                                                  >  2:1265466/1‑149 (MQ=255)
                                   tcCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCaa                                                                                                                  >  1:1409570/1‑149 (MQ=255)
                                   tcCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCaa                                                                                                                  <  4:366925/149‑1 (MQ=255)
                                     cTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                >  3:1376507/1‑149 (MQ=255)
                                      ttGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAg                                                                                                               >  4:1408616/1‑149 (MQ=255)
                                        gATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaat                                                                                                                                     <  2:472221/125‑1 (MQ=255)
                                         atatAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGa                                                                                                            >  1:631951/1‑149 (MQ=255)
                                          tataATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGaa                                                                                                           <  1:696798/149‑1 (MQ=255)
                                           ataATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGaaa                                                                                                          <  1:327219/149‑1 (MQ=255)
                                             aaTCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAata                                                                                                        >  3:948475/1‑149 (MQ=255)
                                               tCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATc                                                                                                      <  2:631951/149‑1 (MQ=255)
                                                  tAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGc                                                                                                                               >  4:247305/1‑121 (MQ=255)
                                                   aaaaaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaata                                                                                                                    <  2:1584185/131‑5 (MQ=255)
                                                     aaaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTattatt                                                                                                                      <  1:749963/127‑1 (MQ=255)
                                                     aaaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaata                                                                                                                    >  2:1002185/1‑125 (MQ=255)
                                                        aaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTTTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAcgcg                                                                                             <  2:258728/149‑1 (MQ=255)
                                                        aaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAcgcg                                                                                             <  2:1452516/149‑1 (MQ=255)
                                                          tctcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCa                                                                                           <  2:1346470/149‑1 (MQ=255)
                                                              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGACAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAg                                                                                          >  1:590920/1‑146 (MQ=255)
                                                              aTAAAATTAATATATGAGAGAATCTTAATTCAGCAGAAGATTATTAAAGGTTGCTGTGTTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTACTAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGTAGGaa                                                                                       >  2:109000/1‑149 (MQ=255)
                                                                aaaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACa                                                                                     <  4:74553/149‑1 (MQ=255)
                                                                   aTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGtt                                                                                  >  1:1497672/1‑149 (MQ=255)
                                                                    ttAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTc                                                                                 >  1:1443677/1‑149 (MQ=255)
                                                                          tatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAg                                                                           >  4:160506/1‑149 (MQ=255)
                                                                           atGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGt                                                                          <  2:762988/149‑1 (MQ=255)
                                                                            tGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGtt                                                                         <  3:848648/149‑1 (MQ=255)
                                                                             gagaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGttt                                                                        >  1:446660/1‑149 (MQ=255)
                                                                              agaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtc                                                                       >  4:1250270/1‑149 (MQ=255)
                                                                               gaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtct                                                                      >  3:381742/1‑149 (MQ=255)
                                                                                aTAATCTTTATTCAGCATAAGATGATTATAGGTTGCTGTATTATTTATCGATAAAAAAAGACTGCGAGATGGTAGGCTATTTAATAACGGCGTTATTATTGCAACAGAGAAATATCCGCAACGCGCAGAACTAGTTCGAGCAGTTtctc                                                                     >  2:584496/1‑149 (MQ=255)
                                                                                aTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtctc                                                                     <  2:410958/149‑1 (MQ=255)
                                                                                 tAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTcc                                                                    >  2:1196480/1‑149 (MQ=255)
                                                                                  aaTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCa                                                                   <  4:485651/149‑1 (MQ=255)
                                                                                   aTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTattatt                                                                                                                      <  3:351301/97‑1 (MQ=255)
                                                                                    tCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGtt                                                                                  <  3:1502519/132‑1 (MQ=255)
                                                                                      tttATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatta                                                                                                                        >  4:906315/1‑92 (MQ=255)
                                                                                              gcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc                                                       >  2:978942/1‑149 (MQ=255)
                                                                                                  aagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATg                                                              <  2:1038404/138‑1 (MQ=255)
                                                                                                  aagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGtt                                                   >  3:579853/1‑149 (MQ=255)
                                                                                                     attattAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGAt                                                >  4:1174468/1‑149 (MQ=255)
                                                                                                      ttattaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATa                                               <  1:1131939/149‑1 (MQ=255)
                                                                                                         ttaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTattatt                                                                                                                      >  4:642908/1‑75 (MQ=255)
                                                                                                         ttaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGc                                            <  1:797013/149‑1 (MQ=255)
                                                                                                         ttaAAGGTTGATGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTattatt                                                                                                                      <  3:642908/75‑1 (MQ=255)
                                                                                                           aaaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCaa                                          <  3:214928/149‑1 (MQ=255)
                                                                                                           aaaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCaa                                          >  1:1001048/1‑149 (MQ=255)
                                                                                                               gTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGCCAGACGGTTGATACGCAATTCt                                      <  2:1133147/149‑1 (MQ=255)
                                                                                                                   cTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGt                                  <  4:839599/149‑1 (MQ=255)
                                                                                                                    tGTATTATTTAGCGATAAAAAAAGCGTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtctc                                                                     <  3:159261/113‑1 (MQ=255)
                                                                                                                      tattatTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTg                                   >  2:661409/1‑145 (MQ=255)
                                                                                                                      tattatTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtcat                               >  3:806665/1‑149 (MQ=255)
                                                                                                                      tattatTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtcat                               >  4:787494/1‑149 (MQ=255)
                                                                                                                               gCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCttt                                    >  1:1078334/1‑135 (MQ=255)
                                                                                                                                cGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGAc                     <  3:692049/149‑1 (MQ=255)
                                                                                                                                 gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGAcc                    >  4:344237/1‑149 (MQ=255)
                                                                                                                                 gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGAcc                    >  2:696416/1‑149 (MQ=255)
                                                                                                                                  aTAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGt                                                    >  4:95409/1‑116 (MQ=255)
                                                                                                                                    aaaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTc                                                                                 <  4:1078435/85‑1 (MQ=255)
                                                                                                                                     aaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCa                >  2:525937/1‑149 (MQ=255)
                                                                                                                                      aaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCAc               >  4:1329213/1‑149 (MQ=255)
                                                                                                                                       aataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACt              >  1:1293328/4‑149 (MQ=255)
                                                                                                                                       aataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACt              >  4:294736/4‑149 (MQ=255)
                                                                                                                                       aaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtctc                                                                     <  4:1024493/94‑1 (MQ=255)
                                                                                                                                       aaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtctc                                                                     >  3:1024493/1‑94 (MQ=255)
                                                                                                                                       aaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACt              <  2:1534572/149‑1 (MQ=255)
                                                                                                                                        ataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACtt             >  1:476641/3‑149 (MQ=255)
                                                                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCttt                                    >  2:919629/1‑125 (MQ=255)
                                                                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCttt                                    >  3:1096323/1‑125 (MQ=255)
                                                                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACttt            >  4:417229/1‑149 (MQ=255)
                                                                                                                                           gCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTAt          >  4:593545/1‑149 (MQ=255)
                                                                                                                                            ccTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATc         <  1:1560634/149‑1 (MQ=255)
                                                                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa      >  2:1306156/1‑149 (MQ=255)
                                                                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa      <  3:1310577/149‑1 (MQ=255)
                                                                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa      <  1:1025274/149‑1 (MQ=255)
                                                                                                                                                ccAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaag     >  3:428392/1‑149 (MQ=255)
                                                                                                                                                   gATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGaaa  <  2:213949/149‑1 (MQ=255)
                                                                                                                                                    aTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCttt                                    <  3:1187903/114‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGAAA  >  NC_000913/3722153‑3722448

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