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.5% 689.4 / NA 217intergenic (‑296/‑101)b2965/yqgAornithine decarboxylase isozyme/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (92/123);  minor base G (0/1);  total (92/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

GAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGC  >  NC_000913/3107377‑3107555
                                                                                              |                                                                                      
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         <  1:3707903/93‑1 (MQ=255)
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         >  2:2654699/1‑93 (MQ=255)
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         >  1:3460657/1‑93 (MQ=255)
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         >  1:4418801/1‑93 (MQ=255)
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         >  1:2154355/1‑93 (MQ=255)
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         <  2:1451064/93‑1 (MQ=255)
 aaaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                                        <  2:3171867/93‑2 (MQ=255)
 aaaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                                        <  1:1737574/93‑2 (MQ=255)
   aCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcgc                                                                                      <  1:872219/93‑4 (MQ=255)
   aCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcgc                                                                                      >  2:654487/1‑90 (MQ=255)
        gATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                                 >  2:2745679/1‑93 (MQ=255)
         aTAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGc                                                                                <  2:2598823/93‑1 (MQ=255)
          tAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTTAAGCg                                                                               >  1:4177582/1‑93 (MQ=255)
          tAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGAAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATAATGGCACCTTTCGCTAAAGCg                                                                               >  2:561419/1‑93 (MQ=255)
            aGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTg                                                                             <  1:1488753/93‑1 (MQ=255)
               gagagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                                          <  1:2436040/93‑1 (MQ=255)
               gagagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATAGCACCTTTCGCTAAAGCGTGCaa                                                                          >  2:3171692/1‑93 (MQ=255)
                 gagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGc                                                                        <  1:3154132/93‑1 (MQ=255)
                 gagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCATACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGc                                                                        >  2:828901/1‑93 (MQ=255)
                  agaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGcc                                                                       <  1:3916882/93‑1 (MQ=255)
                  agaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGcc                                                                       <  1:1352015/93‑1 (MQ=255)
                  agaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGcc                                                                       >  1:4580148/1‑93 (MQ=255)
                  agaAGGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGcc                                                                       >  2:2093769/1‑93 (MQ=255)
                  aaaaTTTAAGGAGATGGCGTGCATCGTGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCATTCATGGCACCTTTCGCTAAAGCGTGCAAGcc                                                                       <  2:2569522/91‑1 (MQ=255)
                   gaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                                      >  1:1680570/1‑93 (MQ=255)
                   gaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                                      <  2:409068/93‑1 (MQ=255)
                     aTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                                    <  2:2027029/93‑1 (MQ=255)
                     aTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                                    >  1:2725614/1‑93 (MQ=255)
                     aTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGCGCAAGCCAgg                                                                    >  1:816609/1‑93 (MQ=255)
                       gTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAc                                                                  <  2:947240/93‑1 (MQ=255)
                        tAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACt                                                                 >  1:2811607/1‑93 (MQ=255)
                         aaGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACtt                                                                >  2:945329/1‑93 (MQ=255)
                         aaGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACtt                                                                >  2:10078/1‑93 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               >  1:2129776/1‑93 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               <  2:34002/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               <  2:1202524/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               <  1:3187695/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               >  2:3846174/1‑93 (MQ=255)
                            gagaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGc                                                             <  2:2909892/93‑1 (MQ=255)
                              gaTGGCGTGCATCGGGCAAACTCTGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGa                                                           >  1:4647074/1‑93 (MQ=255)
                               aTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaa                                                          >  1:4284203/1‑93 (MQ=255)
                                   cGTGCATCGGGCAACCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaa                                                      <  2:2428752/93‑1 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    >  2:3710230/1‑93 (MQ=255)
                                       cATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAg                                                  <  1:2945229/93‑1 (MQ=255)
                                       cATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAg                                                  >  2:4523313/1‑93 (MQ=255)
                                        aTCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGc                                                 <  1:1164567/93‑1 (MQ=255)
                                         tCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                                                <  1:2032429/93‑1 (MQ=255)
                                         tCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                                                >  2:1612800/1‑93 (MQ=255)
                                           gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                                              >  2:3112889/1‑93 (MQ=255)
                                            ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATc                                             <  2:817447/93‑1 (MQ=255)
                                            ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATc                                             <  2:4289394/93‑1 (MQ=255)
                                             gCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCg                                            >  2:1295407/1‑93 (MQ=255)
                                             gCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCg                                            >  2:1105647/1‑93 (MQ=255)
                                             gCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCg                                            >  1:4467094/1‑93 (MQ=255)
                                              cAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGa                                           <  2:596029/93‑1 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         >  2:280103/1‑93 (MQ=255)
                                                  cTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                                       >  1:1346592/1‑93 (MQ=255)
                                                   tCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCa                                      >  1:606199/1‑93 (MQ=255)
                                                     cGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAcc                                    <  1:1051637/93‑1 (MQ=255)
                                                     cGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAcc                                    >  2:3510818/1‑93 (MQ=255)
                                                     cGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAcc                                    >  1:1640753/1‑93 (MQ=255)
                                                      gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCa                                   <  1:4148245/93‑1 (MQ=255)
                                                      gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCa                                   <  2:2518758/93‑1 (MQ=255)
                                                       tACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaa                                  <  1:1993386/93‑1 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaaa                                 <  2:4665129/93‑1 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaaa                                 >  1:4220377/1‑93 (MQ=255)
                                                          agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACg                               <  2:1416851/93‑1 (MQ=255)
                                                          agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACg                               >  2:1033060/1‑93 (MQ=255)
                                                           gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg                              >  2:490364/1‑93 (MQ=255)
                                                           gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg                              >  2:3239731/1‑93 (MQ=255)
                                                            agagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAGGCCATCGACGGCACCAAACggc                             >  1:953276/1‑92 (MQ=255)
                                                              aggagAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCa                           >  1:3879088/1‑93 (MQ=255)
                                                                gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata                         <  2:3652843/93‑1 (MQ=255)
                                                                gagaAAATTCGCGCAATATTGGCACGTTGCCCTAAAGCGTGCAAGCCAGGACTTTGCGAACAACCAAGCCATCGACGCCTCCAAACGGTCata                         <  2:287578/93‑1 (MQ=255)
                                                                  gaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCatact                       >  1:140332/1‑91 (MQ=255)
                                                                  gaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataat                       <  1:1889066/93‑1 (MQ=255)
                                                                    aaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa                     >  2:1946803/1‑93 (MQ=255)
                                                                    aaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa                     >  1:4384986/1‑93 (MQ=255)
                                                                     aaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAg                    >  2:3332701/1‑93 (MQ=255)
                                                                     aaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAg                    <  1:1281574/93‑1 (MQ=255)
                                                                      aTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGa                   <  2:1337673/93‑1 (MQ=255)
                                                                       ttCGCGCAATCATGGCACCTTTGGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaa                  <  2:3951697/93‑1 (MQ=255)
                                                                       ttCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaa                  <  2:1614221/93‑1 (MQ=255)
                                                                        tCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaaa                 <  1:1141902/93‑1 (MQ=255)
                                                                         cgcgCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaaaa                <  1:3871457/93‑1 (MQ=255)
                                                                           cgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATc              <  2:666796/93‑1 (MQ=255)
                                                                             cAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCaa            <  2:2150054/93‑1 (MQ=255)
                                                                              aaTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCaaa           >  1:3379544/1‑93 (MQ=255)
                                                                              aaTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCaaa           >  2:1846448/1‑93 (MQ=255)
                                                                               aTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAAc          >  2:494747/1‑93 (MQ=255)
                                                                               aTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAAc          <  1:1456138/93‑1 (MQ=255)
                                                                                  aTGGCACCTTTCGCTAAAGCGTGCAGGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAAt       <  2:4169917/93‑1 (MQ=255)
                                                                                  aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAAt       >  1:3749826/1‑93 (MQ=255)
                                                                                  aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAAt       <  1:3016844/93‑1 (MQ=255)
                                                                                   tGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCACAATAAGAAAATCAAACAATa      <  2:4572679/93‑1 (MQ=255)
                                                                                   tGGCACCTTTCGCTAAAGCATGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATa      <  2:1787627/93‑1 (MQ=255)
                                                                                    ggCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATAc     <  2:4423992/93‑1 (MQ=255)
                                                                                    ggCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATAc     <  2:2198439/93‑1 (MQ=255)
                                                                                     gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACa    <  2:577201/93‑1 (MQ=255)
                                                                                       aCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGc  >  2:4716979/1‑93 (MQ=255)
                                                                                              |                                                                                      
GAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGC  >  NC_000913/3107377‑3107555

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