Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,306 C→G 15.8% 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,3060CG15.8% 293.2 / 42.9 107G73R (GGC→CGC) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (37/53);  new base G (10/7);  total (47/60)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 114 total aligned reads displayed.

CTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGT  >  NC_000913/1950207‑1950404
                                                                                                    |                                                                                                  
cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGcc                                                                                                    <  2:1457097/100‑1 (MQ=255)
cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGcc                                                                                                    >  1:1506951/1‑100 (MQ=255)
 tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGccc                                                                                                   >  2:2694375/1‑100 (MQ=255)
  gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCg                                                                                                  <  2:2522623/100‑1 (MQ=255)
  gcgcGGTTGATGATAGGCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCg                                                                                                  <  1:1595/100‑1 (MQ=255)
    gtggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                                <  1:2532019/98‑1 (MQ=255)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                              <  2:918225/100‑1 (MQ=38)
       gTTGATGATGGTCAGCCAGGACGCCAGCACTTCGATTTCGCCGGTAGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACTTACGGTGccc                                                                                                   <  1:1575767/94‑1 (MQ=255)
       gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                                             <  1:2683835/100‑2 (MQ=38)
        ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcga                                                                                            <  2:1618348/100‑3 (MQ=37)
        ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcga                                                                                            <  2:1444875/100‑3 (MQ=37)
         tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                           <  2:2174131/100‑1 (MQ=255)
         tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                           <  1:1445847/100‑1 (MQ=255)
          gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGa                                                                                          >  2:468983/1‑100 (MQ=255)
           atgatAGTCAGCGAGGACGCCAGCCCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcga                                                                                            <  2:1692668/97‑3 (MQ=37)
           atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGaa                                                                                         >  2:1046772/1‑100 (MQ=255)
           atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGaa                                                                                         >  1:1083296/1‑100 (MQ=255)
            tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAAt                                                                                        <  2:1798733/100‑1 (MQ=255)
            tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAAt                                                                                        <  1:279650/100‑1 (MQ=255)
            tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAAt                                                                                        <  1:1328070/100‑1 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGTCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                   <  1:1602112/100‑1 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                   >  2:1432970/1‑100 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                   >  2:1812993/1‑100 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                   >  1:1238387/1‑100 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                                   <  1:856783/100‑1 (MQ=255)
                   cAGCGAGGACGCCAACACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                 <  1:2457412/100‑1 (MQ=255)
                    aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                 >  2:2102720/1‑99 (MQ=255)
                    aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                                >  1:1169273/1‑100 (MQ=255)
                     gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                               <  1:1184462/100‑1 (MQ=255)
                      cgagtccgcgAGCACTTCGTTTTCGCGGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGGCCTGAATGCAGAAc                                                                                 <  1:2380262/91‑1 (MQ=255)
                       gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                             >  1:1118634/1‑100 (MQ=255)
                          gACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                          <  2:547978/100‑1 (MQ=255)
                          gACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                          <  2:1395997/100‑1 (MQ=255)
                           tctccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                          <  2:1522335/96‑1 (MQ=255)
                            cgcccgccCTTCGATTTCGCCGGTCGCCATATCGCGGTTGATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                        <  1:2199448/100‑1 (MQ=255)
                            cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATCTTTTTTTCGGCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                        <  1:999349/100‑1 (MQ=255)
                            cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                        >  1:1050848/1‑100 (MQ=255)
                            cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                        >  2:5283/1‑100 (MQ=255)
                            cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAACGCAGAACTCATTACGc                                                                        >  1:1702978/1‑100 (MQ=255)
                             gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCa                                                                       <  1:806800/100‑1 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAcagat                                                                      >  1:1522331/1‑96 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                      >  1:518621/1‑100 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCCGAATGCCGAACTCATTACGCAg                                                                      >  1:1393902/1‑100 (MQ=255)
                                   aCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                 >  2:1750930/1‑100 (MQ=255)
                                   aCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                 >  2:267635/1‑100 (MQ=255)
                                      tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                               <  2:1276/99‑1 (MQ=255)
                                      tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                              <  2:2388104/100‑1 (MQ=255)
                                        gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                            <  2:1862205/100‑1 (MQ=255)
                                                cGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                    <  1:2572372/100‑1 (MQ=255)
                                                cGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                    <  2:349664/100‑1 (MQ=255)
                                                    cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                                <  2:1199629/100‑1 (MQ=255)
                                                    cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAATCATTACGCAGTTAACAAGACAGCtttctcacg                                                >  2:1257226/1‑94 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                               >  2:2863802/1‑100 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                               >  2:584351/1‑100 (MQ=255)
                                                      ccATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                              >  1:686252/1‑100 (MQ=255)
                                                      ccATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                              <  1:828281/100‑1 (MQ=255)
                                                      ccATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                              <  1:2165109/100‑1 (MQ=255)
                                                       cATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                              <  1:2820432/99‑1 (MQ=255)
                                                          atCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                          >  1:813014/1‑100 (MQ=255)
                                                           tCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                         <  2:919920/100‑1 (MQ=255)
                                                            cgcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACg                                        >  2:2416748/1‑100 (MQ=255)
                                                             gcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                                       >  2:888087/1‑100 (MQ=255)
                                                                  tAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                  <  1:736622/100‑1 (MQ=255)
                                                                    atatTTTTTTTTTTATTCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                >  1:2388242/1‑100 (MQ=255)
                                                                    atatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                >  2:1377988/1‑100 (MQ=255)
                                                                     tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                               <  2:1651284/100‑1 (MQ=255)
                                                                     tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                               <  2:33652/100‑1 (MQ=255)
                                                                     tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                               <  1:267633/100‑1 (MQ=255)
                                                                      atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                              <  2:2158426/100‑1 (MQ=255)
                                                                      atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                              <  2:696283/100‑1 (MQ=255)
                                                                      atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAGCCGGACCg                              >  2:1264323/1‑100 (MQ=255)
                                                                       tttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                             <  1:1430714/100‑1 (MQ=255)
                                                                         tttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                           <  2:19770/100‑1 (MQ=255)
                                                                             atctCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                       <  1:553523/97‑1 (MQ=38)
                                                                              ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                       >  1:1837793/2‑99 (MQ=37)
                                                                              ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                      >  2:1021586/2‑100 (MQ=38)
                                                                              cGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                      <  2:1622008/100‑1 (MQ=255)
                                                                               gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGACCTCATTACGCAGTTCAGAGGCCAGCTTTAACGCTTCCGCACGATCAGGATCTAAAAATAcc                     >  1:2608915/1‑100 (MQ=25)
                                                                               gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                     >  2:590288/1‑100 (MQ=38)
                                                                               gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                     >  1:1372367/1‑100 (MQ=38)
                                                                               gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                      >  1:2555189/1‑99 (MQ=38)
                                                                               gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                     <  1:1050004/100‑1 (MQ=255)
                                                                               gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                     >  2:671919/1‑100 (MQ=255)
                                                                               gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                     <  1:779425/100‑1 (MQ=255)
                                                                                tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                    >  1:1568113/1‑100 (MQ=38)
                                                                                tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                    >  1:1386248/1‑100 (MQ=38)
                                                                                tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                    >  2:2883551/1‑100 (MQ=38)
                                                                                tCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                    >  2:2158797/1‑100 (MQ=255)
                                                                                tCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                    <  2:1391285/100‑1 (MQ=255)
                                                                                 cACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATaagat                   >  2:553524/1‑96 (MQ=255)
                                                                                 cACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                   >  1:699493/1‑100 (MQ=255)
                                                                                   cgcgCACGTACGGTGCCCGTGCCCTGAATGCATAACTCATTAGGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                   <  1:2345721/98‑1 (MQ=255)
                                                                                   cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa                 >  2:1362633/1‑100 (MQ=255)
                                                                                       cACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGAt             <  2:762667/100‑1 (MQ=255)
                                                                                            tACGGTGCCTGTGACCTGAATTTATAACTCATTACGCAGTTAAGACGCCAGCTTTAACGCGTCCGCACGATCCTGATCGAAAAATACCTGCACGATa            <  1:1292257/97‑1 (MQ=255)
                                                                                               ggTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcg      >  2:847620/1‑100 (MQ=255)
                                                                                               ggTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTc       >  1:954144/1‑99 (MQ=255)
                                                                                                   cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGt  >  2:166454/1‑100 (MQ=255)
                                                                                                   cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGt  <  2:497796/100‑1 (MQ=255)
                                                                                                   cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCAAGCTTTAACGCGTACGCACGATCCGGATGGAAAAAAAAATGCACGATACCTTCGAgg   >  2:2640431/1‑99 (MQ=255)
                                                                                                    |                                                                                                  
CTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGT  >  NC_000913/1950207‑1950404

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: