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.C98.9% 1097.7 / 2.7 355intergenic (‑296/‑101)b2965/yqgAornithine decarboxylase isozyme/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base C (152/199);  total (153/202)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 392 total aligned reads displayed.

GAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGC  >  NC_000913/3107377‑3107555
                                                                                              |                                                                                      
gAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                                         <  1:1427458/93‑1 (MQ=255)
 aaaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGACAATTCGCGCACTCATGGCACCtttc                                                                                        >  1:6413677/1‑92 (MQ=255)
 aaaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                                        <  1:2788073/93‑2 (MQ=255)
 aaaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                                        <  2:3358551/93‑2 (MQ=255)
  aaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcg                                                                                       <  2:5565887/93‑3 (MQ=255)
  aaCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcg                                                                                       <  2:1070566/93‑3 (MQ=255)
   aCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcgc                                                                                      <  1:3323015/93‑4 (MQ=255)
   aCGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcgc                                                                                      >  1:6719196/1‑90 (MQ=255)
      gCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaa                                                                                   >  1:4342828/1‑93 (MQ=255)
       cGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaaa                                                                                  <  2:1208962/93‑1 (MQ=255)
       cGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaaa                                                                                  <  2:1350125/93‑1 (MQ=255)
        gATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                                 >  2:4740668/1‑93 (MQ=255)
        gATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                                 >  2:2821913/1‑93 (MQ=255)
        gATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                                 <  1:3473582/93‑1 (MQ=255)
         aTAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGc                                                                                <  2:3691264/93‑1 (MQ=255)
         aTAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGc                                                                                <  1:4551085/93‑1 (MQ=255)
         aTAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGc                                                                                >  2:2404737/1‑93 (MQ=255)
          tAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCg                                                                               >  2:6513790/1‑93 (MQ=255)
          tAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCg                                                                               >  1:4696779/1‑93 (MQ=255)
            aGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTg                                                                             <  1:1645339/93‑1 (MQ=255)
              ggagagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCa                                                                           <  2:1743189/92‑1 (MQ=255)
                agagaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAg                                                                         >  1:6004445/1‑93 (MQ=255)
                    aaTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAg                                                                     >  2:6126701/1‑93 (MQ=255)
                     aTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                                    >  1:3849827/1‑93 (MQ=255)
                      tGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGa                                                                   >  2:3084524/1‑93 (MQ=255)
                      tGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGa                                                                   <  1:1105846/93‑1 (MQ=255)
                         aaGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACtt                                                                <  1:5313066/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               <  1:4976618/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               <  2:4638180/93‑1 (MQ=255)
                          aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                                               <  1:6515412/93‑1 (MQ=255)
                           ggAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTg                                                              >  2:2524941/1‑93 (MQ=255)
                            gagaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGc                                                             <  2:814412/93‑1 (MQ=255)
                              gaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGa                                                           >  1:2222917/1‑93 (MQ=255)
                                tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaac                                                         >  1:2711355/1‑93 (MQ=255)
                                 ggCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaaca                                                        >  1:2122993/1‑93 (MQ=255)
                                 ggCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaaca                                                        <  2:7481226/93‑1 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    >  1:6294323/1‑93 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    >  1:7598830/1‑93 (MQ=255)
                                     tGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                                    >  2:6747849/1‑93 (MQ=255)
                                       cATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAg                                                  <  1:3681049/93‑1 (MQ=255)
                                         tCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                                                <  1:2023877/93‑1 (MQ=255)
                                         tCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                                                >  2:5948830/1‑93 (MQ=255)
                                          cGGACAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCa                                               <  2:3930607/93‑1 (MQ=255)
                                           gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                                              <  1:3533177/93‑1 (MQ=255)
                                           gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                                              <  1:2584080/93‑1 (MQ=255)
                                           gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                                              <  1:7085252/93‑1 (MQ=255)
                                           gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                                              <  2:7178292/93‑1 (MQ=255)
                                           gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                                              <  2:7773897/93‑1 (MQ=255)
                                            ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATc                                             <  2:5363363/93‑1 (MQ=255)
                                            ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATc                                             <  1:58035/93‑1 (MQ=255)
                                            ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATc                                             <  1:4116777/93‑1 (MQ=255)
                                              gaaaCTCCGTACAGAGAGGCAAAAATACGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCAAGAACTTTGCGAACAAACAAGCCATcaa                                           <  1:7034508/92‑3 (MQ=255)
                                               aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                                          <  2:4206039/93‑1 (MQ=255)
                                               aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                                          <  1:515261/93‑1 (MQ=255)
                                               aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                                          <  2:2259862/93‑1 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         <  2:6022000/93‑1 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         >  1:6931589/1‑93 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         >  1:4897452/1‑93 (MQ=255)
                                                aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                                         <  1:4681365/93‑1 (MQ=255)
                                                  cTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                                       >  1:1176416/1‑93 (MQ=255)
                                                   tCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCa                                      >  2:4584309/1‑93 (MQ=255)
                                                   tCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCa                                      >  1:1589123/1‑93 (MQ=255)
                                                      gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCa                                   >  1:990898/1‑93 (MQ=255)
                                                      gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCa                                   >  1:6829524/1‑93 (MQ=255)
                                                       tACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaa                                  >  1:1141705/1‑93 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGCCTTTGCGAACAACCAAGCCATCGACGGCACGaaa                                 >  2:6782477/1‑93 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaaa                                 <  2:3484938/93‑1 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCaaa                                 <  1:4298097/93‑1 (MQ=255)
                                                        aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACCGCACCaaa                                 >  1:1324311/1‑93 (MQ=255)
                                                         cAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc                                >  2:1495264/1‑93 (MQ=255)
                                                         cAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc                                <  1:3024275/93‑1 (MQ=255)
                                                         cAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc                                <  2:627360/93‑1 (MQ=255)
                                                           gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg                              >  2:7627615/1‑93 (MQ=255)
                                                           gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg                              >  2:5688429/1‑93 (MQ=255)
                                                            agagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGt                             <  1:3382609/93‑1 (MQ=255)
                                                               ggagAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCat                          <  2:1992144/93‑1 (MQ=255)
                                                               ggagAAAATGCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCat                          <  2:1835227/93‑1 (MQ=255)
                                                                gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata                         <  1:5660534/93‑1 (MQ=255)
                                                                 agaAAATTCGCGCAATCATGGCCCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataa                        <  2:6661344/93‑1 (MQ=255)
                                                                  gaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataat                       <  2:4514439/93‑1 (MQ=255)
                                                                       ttCTCGCAATCATGGCTCTTTTCGTTAATGCGCGCAAGCCAGGACTTTTCGAACAAACAAGCCATTGACGGCACCAAACGGTCATAATATGaa                  >  2:7273296/1‑93 (MQ=255)
                                                                       ttCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaa                  <  1:6161528/93‑1 (MQ=255)
                                                                       ttCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGaa                  <  2:5793466/93‑1 (MQ=255)
                                                                          gcgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAAt               >  2:5297794/1‑93 (MQ=255)
                                                                             cAATCATGGCGCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCaa            <  2:1412770/93‑1 (MQ=255)
                                                                               aTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAAc          <  2:6461952/93‑1 (MQ=255)
                                                                                tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACa         <  2:3835279/93‑1 (MQ=255)
                                                                                tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACa         >  1:242553/1‑93 (MQ=255)
                                                                                 cATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACaa        >  2:4401319/1‑93 (MQ=255)
                                                                                 cATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACaa        <  1:3244883/93‑1 (MQ=255)
                                                                                  agggCACCTTCCGCTAAAGCGTGCAAGGCAGTACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAAt       <  2:4023102/91‑1 (MQ=255)
                                                                                  aTGGCCCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAAt       <  2:4974459/93‑1 (MQ=255)
                                                                                  aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAAt       >  1:7267401/1‑93 (MQ=255)
                                                                                  aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAAt       <  1:1835501/93‑1 (MQ=255)
                                                                                     gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACa    <  1:7769007/93‑1 (MQ=255)
                                                                                     gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACa    <  1:7471291/93‑1 (MQ=255)
                                                                                     gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACa    >  2:5168669/1‑93 (MQ=255)
                                                                                     gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACa    >  1:1902185/1‑93 (MQ=255)
                                                                                      cACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAg   <  2:1592942/93‑1 (MQ=255)
                                                                                       aCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGc  >  1:5920514/1‑93 (MQ=255)
                                                                                              |                                                                                      
GAAACGGCGATAAGCGAGAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGAAAATCAAACAATACAGC  >  NC_000913/3107377‑3107555

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