Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,306 C→G 18.7% G73R (GGC→CGC)  aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,3060CG18.7% 367.4 / 78.3 151G73R (GGC→CGC) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (56/66);  new base G (21/7);  total (77/74)
Fisher's exact test for biased strand distribution p-value = 6.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 159 total aligned reads displayed.

CTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTC  >  NC_000913/1950207‑1950405
                                                                                                    |                                                                                                   
cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                                     >  1:1886223/1‑100 (MQ=255)
  gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                   <  2:491329/100‑1 (MQ=255)
    gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    <  2:3405596/97‑1 (MQ=255)
    gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 >  1:1037867/1‑100 (MQ=255)
    gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 <  2:2903707/100‑1 (MQ=255)
    gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 <  1:2685801/100‑1 (MQ=255)
     cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                                <  2:721523/100‑1 (MQ=255)
      ggTTGATGATAGTCAGCGAGGCCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGGTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  2:1760657/100‑1 (MQ=255)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  1:1933561/100‑1 (MQ=38)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  1:3255183/100‑1 (MQ=38)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 <  2:2531808/98‑1 (MQ=255)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  2:596708/100‑1 (MQ=255)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  2:1499156/100‑1 (MQ=255)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  1:2078547/100‑1 (MQ=255)
         tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                                            >  2:2096669/1‑100 (MQ=255)
         tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                                            <  1:2491714/100‑1 (MQ=255)
          gttgatAGTCAGCGAGGACGCCGCCGCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                                           <  2:1742274/98‑5 (MQ=25)
          gatgatAGTCAGCGAGGACGTCCGCACTTCGATTTTGCCGGTCGCCATAGCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaa                                                                                            <  2:2545776/99‑4 (MQ=25)
          gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGa                                                                                           >  1:3239512/1‑100 (MQ=255)
          gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGa                                                                                           >  1:2012101/1‑100 (MQ=255)
           atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGaa                                                                                          >  1:1045292/1‑100 (MQ=255)
           atgatAGTCAACGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGTGGTTAATATTTTTGTCGTCACGCGCA‑CGTACGGTGCTCGTGACCTGaa                                                                                          <  2:860050/100‑1 (MQ=255)
              atAGTCGGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                                       <  2:3278303/100‑1 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                    <  1:2922119/100‑1 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                    >  2:3816176/1‑100 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                    >  2:3752629/1‑100 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                    >  2:2454837/1‑100 (MQ=255)
                   cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                  >  1:3472413/1‑100 (MQ=255)
                    aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                  >  2:65501/1‑99 (MQ=255)
                     gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                >  1:2954564/1‑100 (MQ=255)
                     gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                <  2:1578915/100‑1 (MQ=255)
                      cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                               <  1:1358589/100‑1 (MQ=255)
                       gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                              <  2:366225/100‑1 (MQ=255)
                       gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                              <  2:2035290/100‑1 (MQ=255)
                         ggACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                            <  1:1784162/100‑1 (MQ=255)
                          gACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                           <  1:2455492/100‑1 (MQ=255)
                           acaccAGCAACTCGATCTCGCCGGTCGCCATATCGAGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                           <  1:2555943/96‑1 (MQ=255)
                           aCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGGA‑CGTACGGTGCCTGTGTCCTGAATGGAGAACTCATTACg                                                                          <  2:1012886/100‑1 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       >  2:3116657/1‑100 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       >  2:2692052/1‑100 (MQ=255)
                               cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                      <  1:2637659/100‑1 (MQ=255)
                                aGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGtt                                                                     <  1:395855/100‑1 (MQ=255)
                                    cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 >  2:900624/1‑100 (MQ=255)
                                    cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 >  1:3044371/1‑100 (MQ=255)
                                    cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 >  1:81245/1‑100 (MQ=255)
                                        gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                             <  2:2009232/100‑1 (MQ=255)
                                        gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                             >  1:1487079/1‑100 (MQ=255)
                                         aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                            >  1:2821569/1‑100 (MQ=255)
                                         aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGGGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                            >  1:3762086/1‑100 (MQ=255)
                                          tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                                           >  1:2401693/1‑100 (MQ=255)
                                           ttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                                          <  2:1113522/100‑1 (MQ=255)
                                                 ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                    >  1:3332842/1‑100 (MQ=255)
                                                 ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                    <  2:3368156/100‑1 (MQ=255)
                                                   tCGCCATATCGCGGTTAATATTTTTTACGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGATCTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                    <  2:1535272/98‑1 (MQ=255)
                                                   tCGCCATATCGCGGTTAATATCTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgc                                                  <  1:380982/100‑1 (MQ=255)
                                                    cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                                 >  1:1120353/1‑100 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGGATGCAGAAATCATTACGAAGTTCCGAGGCCAGCTTTAAcgag                                                 >  2:3272916/1‑97 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                                >  1:199924/1‑100 (MQ=255)
                                                      ccATATCGCGGTTAATATTTTTTTTTTTTTTTTTA‑CGTACGGTGCCCGTGACCTGAATTAATAATTTATTATGATGTTCAGAGGCCAGCTTTAAcgc                                                  >  2:839945/1‑97 (MQ=255)
                                                          atCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                           <  1:2096680/100‑1 (MQ=255)
                                                             gcgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                                        >  1:861467/1‑100 (MQ=255)
                                                              cgGTTAATATTTTTTCCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATCCAGAACTCATTACGGAGTTCAGAGGCCAGCTTTAACACGTCAGCACGa                                        >  2:3632530/1‑99 (MQ=255)
                                                               ggTTAATATTTTTTTCGTCCCGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGGGGCCAGCTTTAACGCGTCCCCACGacc                                      >  2:955989/1‑98 (MQ=255)
                                                                  tAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATTCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATc                                      >  1:2268514/1‑97 (MQ=255)
                                                                  tAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                   <  2:1411997/100‑1 (MQ=255)
                                                                  tAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                   <  2:1047462/100‑1 (MQ=255)
                                                                    atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                 >  1:1237901/1‑100 (MQ=255)
                                                                      atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                               <  1:187308/100‑1 (MQ=255)
                                                                       tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                              <  2:2306498/100‑1 (MQ=255)
                                                                       tttttttCGCCACGCGCA‑CTTCCGGTGCCCGTGCCCTGAATGCATAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                <  2:445254/98‑1 (MQ=255)
                                                                        ttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaa                             <  1:3014589/100‑1 (MQ=255)
                                                                             aggccACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCCGTTCAGAGGCCAGCTGTAACGCGTCCGCACGATCCGGATCGAAAAATa                        >  1:1776018/5‑100 (MQ=37)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  2:1341458/2‑100 (MQ=38)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  2:1543523/2‑100 (MQ=38)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  2:1974390/2‑100 (MQ=38)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCAGATCggaagagcg                       >  2:2697374/2‑92 (MQ=25)
                                                                              acccACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGACCTCATTACGCAGTTCCCAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaa                           >  2:229978/4‑96 (MQ=25)
                                                                               gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      >  1:2731618/1‑100 (MQ=38)
                                                                               gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      <  1:519119/100‑1 (MQ=255)
                                                                               gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTGTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      >  2:2219976/1‑100 (MQ=255)
                                                                                tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                     >  2:3843539/1‑100 (MQ=38)
                                                                                tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                     >  1:672363/1‑100 (MQ=38)
                                                                                tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                     >  1:2738891/1‑100 (MQ=38)
                                                                                 cACGGGCACCGTAC‑GTGCGCGTGACCTGAGTGCAGAACCCATTACGCAGTGCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    >  1:2288822/1‑100 (MQ=255)
                                                                                 cACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    >  1:1204922/1‑100 (MQ=255)
                                                                                   cgggcaccGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa                  >  2:2050102/8‑100 (MQ=255)
                                                                                   cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc                   >  2:3773883/1‑99 (MQ=255)
                                                                                   cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa                  <  1:3929119/100‑1 (MQ=255)
                                                                                   cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa                  >  1:3717063/1‑100 (MQ=255)
                                                                                     ggcaccGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACg                >  1:3162551/6‑100 (MQ=255)
                                                                                        a‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa             <  1:3651326/100‑1 (MQ=255)
                                                                                          cGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc            >  1:3399502/1‑100 (MQ=255)
                                                                                            tACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCt          >  2:3164741/1‑100 (MQ=255)
                                                                                             aCGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCtt         <  1:2048369/100‑1 (MQ=255)
                                                                                              cGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTc        >  2:2981679/1‑100 (MQ=255)
                                                                                                gTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgc      >  2:459179/1‑100 (MQ=255)
                                                                                                gTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgc      <  1:2346136/100‑1 (MQ=255)
                                                                                                  gCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgcg     <  2:2881431/99‑1 (MQ=255)
                                                                                                  gCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAAGATACCttagagg    >  1:726029/1‑95 (MQ=255)
                                                                                                    ccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTc  >  2:1270148/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTC  >  NC_000913/1950207‑1950405

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