Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,107,469 +C intergenic (‑296/‑101) b2965 ← / → yqgA ornithine decarboxylase isozyme/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,107,4691.C99.1% 681.4 / ‑1.4 215intergenic (‑296/‑101)b2965/yqgAornithine decarboxylase isozyme/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (77/136);  total (77/138)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 234 total aligned reads displayed.

GAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGC  >  NC_000913/3107377‑3107555
                                                                                               |                                                                                      
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         <  2:1115070/93‑1 (MQ=255)
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         <  1:6177683/93‑1 (MQ=255)
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         <  1:1628180/93‑1 (MQ=255)
 aaaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                                        <  2:4718957/93‑2 (MQ=255)
 aaaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                                        >  1:5559948/1‑92 (MQ=255)
  aaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcg                                                                                       >  2:5846415/1‑91 (MQ=255)
   aCGGCGATAAGCGAGAGAATGTAAGGAGATGGGGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcgc                                                                                      >  1:1205285/1‑90 (MQ=255)
      gCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaa                                                                                   <  1:1078459/93‑1 (MQ=255)
      gCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaa                                                                                   >  2:1315096/1‑93 (MQ=255)
      gCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaa                                                                                   <  1:4748682/93‑1 (MQ=255)
       cGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaaa                                                                                  <  1:6119406/93‑1 (MQ=255)
       cGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaaa                                                                                  <  2:4406939/93‑1 (MQ=255)
        gATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                                 <  1:5880829/93‑1 (MQ=255)
        gATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                                 >  1:5301049/1‑93 (MQ=255)
            aGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTg                                                                             <  1:1311930/93‑1 (MQ=255)
            aGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTg                                                                             <  2:4854318/93‑1 (MQ=255)
             gCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGc                                                                            >  1:5249812/1‑93 (MQ=255)
             gCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGc                                                                            <  2:5908417/93‑1 (MQ=255)
               gagagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                                          >  2:149781/1‑93 (MQ=255)
                 gagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGc                                                                        <  2:3312796/93‑1 (MQ=255)
                  agaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGcc                                                                       <  1:3784719/93‑1 (MQ=255)
                   gaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                                      <  2:1086959/93‑1 (MQ=255)
                   gaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                                      <  2:4586184/93‑1 (MQ=255)
                    aaTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAg                                                                     <  2:1417399/93‑1 (MQ=255)
                      tGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGa                                                                   <  2:3950125/93‑1 (MQ=255)
                       gTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAc                                                                  >  2:5505315/1‑93 (MQ=255)
                         aaGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACtt                                                                <  1:100954/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               <  1:4347915/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               >  2:1029633/1‑93 (MQ=255)
                               aTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaa                                                          <  1:3456301/93‑1 (MQ=255)
                                 ggCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaaca                                                        <  2:30662/93‑1 (MQ=255)
                                 ggCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaaca                                                        >  1:179464/1‑93 (MQ=255)
                                    gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaac                                                     <  1:1861256/93‑1 (MQ=255)
                                    gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaac                                                     >  2:3052300/1‑93 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCTCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    <  1:3835457/93‑1 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    >  1:1484464/1‑93 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    >  1:3566331/1‑93 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    <  1:2173700/93‑1 (MQ=255)
                                      gCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAacaaacaa                                                   <  2:2710226/93‑1 (MQ=255)
                                       cATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAg                                                  <  2:2383069/93‑1 (MQ=255)
                                       cATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAg                                                  <  2:1676859/93‑1 (MQ=255)
                                        aTCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGc                                                 >  2:3916100/1‑93 (MQ=255)
                                         tCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                                                >  1:1713248/1‑93 (MQ=255)
                                         agtggCAACCCCCGTACAGAGAGGAGAAAATTCGCGAAATCAGGGCCCCTTTCGCAAAGGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                                                <  2:1861368/90‑1 (MQ=255)
                                          cGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCa                                               <  2:1275273/93‑1 (MQ=255)
                                          cGGGAAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCa                                               <  1:78057/93‑1 (MQ=255)
                                           gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                                              >  1:5147549/1‑93 (MQ=255)
                                             gCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCg                                            <  2:2443627/93‑1 (MQ=255)
                                               aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                                          <  1:1477145/93‑1 (MQ=255)
                                               aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                                          >  1:2520683/1‑93 (MQ=255)
                                               aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                                          <  2:5465201/93‑1 (MQ=255)
                                               aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                                          <  1:4952305/93‑1 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         <  1:592227/93‑1 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         >  2:2579403/1‑93 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         <  1:340552/93‑1 (MQ=255)
                                                 aCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACgg                                        <  2:3127520/93‑1 (MQ=255)
                                                  cTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                                       <  1:2147240/93‑1 (MQ=255)
                                                   tCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCa                                      >  1:3045182/1‑93 (MQ=255)
                                                      gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCa                                   <  2:2227842/93‑1 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaaa                                 >  1:4676776/1‑93 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaaa                                 >  1:5434328/1‑93 (MQ=255)
                                                         cAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc                                <  1:5505315/93‑1 (MQ=255)
                                                          agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACg                               <  1:4974004/93‑1 (MQ=255)
                                                           gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg                              >  2:1268805/1‑93 (MQ=255)
                                                           gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg                              <  2:2764914/93‑1 (MQ=255)
                                                           gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg                              >  1:1092237/1‑93 (MQ=255)
                                                            agagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGt                             >  2:1771887/1‑93 (MQ=255)
                                                             gaggagAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGta                            >  1:4857497/1‑92 (MQ=255)
                                                                gagaAAATTCGCGCAATCATGGCCCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata                         <  1:4280450/93‑1 (MQ=255)
                                                                gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata                         >  1:417599/1‑93 (MQ=255)
                                                                gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata                         <  1:397451/93‑1 (MQ=255)
                                                                gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata                         >  1:804823/1‑93 (MQ=255)
                                                                gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata                         <  1:1266842/93‑1 (MQ=255)
                                                                 agaAAATACGCGCAATCAAGGCTCCTGTCGCTAAAGCGTGCAAGCCAGGACTTGGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataa                        <  2:3852709/93‑1 (MQ=255)
                                                                   aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataata                      >  2:5339166/1‑93 (MQ=255)
                                                                   aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataata                      <  2:610257/93‑1 (MQ=255)
                                                                   aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataata                      <  2:5133983/93‑1 (MQ=255)
                                                                    aaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa                     <  1:211889/93‑1 (MQ=255)
                                                                     aaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAg                    <  1:3695765/93‑1 (MQ=255)
                                                                      atgcgcgCAATCATGGCACCTTCCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGa                   <  2:4618794/90‑1 (MQ=255)
                                                                        tCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaaa                 <  2:1735082/93‑1 (MQ=255)
                                                                         cgcgCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaaaa                <  2:2099045/93‑1 (MQ=255)
                                                                         cgcgCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaaaa                <  2:4451652/93‑1 (MQ=255)
                                                                          gcgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAAt               <  1:2844791/93‑1 (MQ=255)
                                                                          gcgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAAt               >  1:6158056/1‑93 (MQ=255)
                                                                            gcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCa             >  1:3209048/1‑93 (MQ=255)
                                                                             cAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCaa            >  2:682423/1‑93 (MQ=255)
                                                                              aaTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCaaa           <  2:3195431/93‑1 (MQ=255)
                                                                                tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACa         <  2:5126825/93‑1 (MQ=255)
                                                                                tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACa         <  2:2027493/93‑1 (MQ=255)
                                                                                tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACa         <  1:3721705/93‑1 (MQ=255)
                                                                                   tGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATa      >  2:4302557/1‑93 (MQ=255)
                                                                                   tGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATa      <  1:656346/93‑1 (MQ=255)
                                                                                   tGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATa      <  1:225858/93‑1 (MQ=255)
                                                                                    ggCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATAc     <  1:5684653/93‑1 (MQ=255)
                                                                                    ggCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATAc     >  1:2928364/1‑93 (MQ=255)
                                                                                    gccACCTTTCGTTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATAc     <  2:5372396/91‑1 (MQ=255)
                                                                                     gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACa    <  2:33997/93‑1 (MQ=255)
                                                                                      cACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAg   >  2:1481623/1‑93 (MQ=255)
                                                                                      cACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAg   <  1:364242/93‑1 (MQ=255)
                                                                                               |                                                                                      
GAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGC  >  NC_000913/3107377‑3107555

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