Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,891,842 C→T 6.2% G465D (GGT→GAT) ‡ yoaA ← putative 5' to 3' DNA helicase implicated in DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,891,8420CT 6.2% 300.5 / 5.3 96G465D (GGT→GAT) ‡yoaAputative 5' to 3' DNA helicase implicated in DNA repair
Reads supporting (aligned to +/- strand):  ref base C (50/40);  new base T (4/2);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGC  >  NC_000913/1891695‑1891987
                                                                                                                                                   |                                                                                                                                                 
agCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAACATAAAACAACGACCg                                                                                                                                                 >  2:109972/1‑150 (MQ=255)
   gcATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttg                                                                                                                                              <  1:1263000/150‑1 (MQ=255)
       tACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAg                                                                                                                                          >  2:949921/1‑150 (MQ=255)
        aCCGGCGCTGACAAATTGCTGCAACAGTAGCCCTTTGCTGGTTTACCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGc                                                                                                                                         >  2:181690/1‑150 (MQ=255)
         ccGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgttttat                                                                                                                                        <  2:904541/150‑5 (MQ=255)
         ccGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCt                                                                                                                                        >  1:503942/1‑150 (MQ=255)
           ggCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTc                                                                                                                                      <  1:1620763/150‑1 (MQ=255)
           ggCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTc                                                                                                                                      >  2:526047/1‑150 (MQ=255)
               cTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatg                                                                                                                                  >  1:258085/1‑150 (MQ=255)
                tGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatga                                                                                                                                 >  2:1367961/1‑150 (MQ=255)
                  aCAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATc                                                                                                                               >  1:423894/1‑150 (MQ=255)
                  aCAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATc                                                                                                                               >  2:498684/1‑150 (MQ=255)
                   cAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGGAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCg                                                                                                                              <  1:487497/150‑1 (MQ=255)
                   cAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCg                                                                                                                              >  2:993147/1‑150 (MQ=255)
                         gctgcAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgca                                                                                                                        <  2:1423861/150‑1 (MQ=255)
                          ctgcAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgcag                                                                                                                       >  2:321721/1‑150 (MQ=255)
                          ctgcAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgcag                                                                                                                       <  1:37959/150‑1 (MQ=255)
                            gcAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgcagca                                                                                                                     >  1:429276/1‑150 (MQ=255)
                                cAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGc                                                                                                                 >  1:312905/1‑150 (MQ=255)
                                      cccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTaa                                                                                                           >  1:519914/1‑150 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAAc                                                                                                          >  2:1376877/1‑150 (MQ=255)
                                           tGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGc                                                                                                      <  1:526047/150‑1 (MQ=255)
                                            gCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCg                                                                                                     <  2:519914/150‑1 (MQ=255)
                                            gCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCg                                                                                                     >  1:781021/1‑150 (MQ=255)
                                            gCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAACATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCg                                                                                                     <  1:109972/150‑1 (MQ=255)
                                             cTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGa                                                                                                    >  2:498152/1‑150 (MQ=255)
                                               ggTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAACCAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGc                                                                                                  >  2:1285132/1‑150 (MQ=255)
                                               ggTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGc                                                                                                  >  2:975265/1‑150 (MQ=255)
                                                gTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCa                                                                                                 >  2:1390135/1‑150 (MQ=255)
                                                gTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCa                                                                                                 <  2:1412985/150‑1 (MQ=255)
                                                 tttCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAg                                                                                                <  2:764974/150‑1 (MQ=255)
                                                   tCCCCCTGCAACAATACAGGAAGCGTCTTGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGaa                                                                                              <  1:1490812/150‑1 (MQ=255)
                                                    cccccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAc                                                                                             <  1:903419/150‑1 (MQ=255)
                                                     ccccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAcc                                                                                            >  2:1380917/1‑150 (MQ=255)
                                                     ccccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAcc                                                                                            >  1:704024/1‑150 (MQ=255)
                                                         tccaacaaTACAGGAAGCGTCATCGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATCGCGTGAGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggc                                                                                        >  2:417239/3‑150 (MQ=255)
                                                         tGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggc                                                                                        <  2:1273384/150‑1 (MQ=255)
                                                          gCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggct                                                                                       >  1:1255530/1‑150 (MQ=255)
                                                            aacaaTACAGGAAGCGGCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGTTCGGTCGCAGCATTGCCGCTAAGTGGCGAGCAGAACctggctgg                                                                                     <  2:1242263/150‑1 (MQ=255)
                                                              caaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggtt                                                                                   <  1:805255/150‑1 (MQ=255)
                                                              caaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggtt                                                                                   <  1:321721/150‑1 (MQ=255)
                                                               aaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttg                                                                                  <  2:900442/150‑1 (MQ=255)
                                                                 tACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttggt                                                                                <  1:1600933/150‑1 (MQ=255)
                                                                  aCAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTggttggtt                                                                               <  2:503942/150‑1 (MQ=255)
                                                                          cGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAg                                                                       <  1:498152/150‑1 (MQ=255)
                                                                          cGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAg                                                                       <  1:181690/150‑1 (MQ=255)
                                                                           gTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGa                                                                      >  2:309293/1‑150 (MQ=255)
                                                                             cATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGAtt                                                                    >  2:723115/1‑150 (MQ=255)
                                                                               tGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCGGAACCTGGCTGGTTGGTTTGCGGCAGATTgc                                                                  >  1:251668/1‑150 (MQ=255)
                                                                                  tAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgca                                                               <  1:1367961/150‑1 (MQ=255)
                                                                                   aGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgcag                                                              <  1:1376877/150‑1 (MQ=255)
                                                                                    gcgcGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgcagc                                                             <  2:34759/150‑1 (MQ=255)
                                                                                        ggAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcacac                                                         <  1:596098/150‑1 (MQ=255)
                                                                                         gAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcacaca                                                        <  2:704024/150‑1 (MQ=255)
                                                                                         gAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcacaca                                                        <  1:778917/150‑1 (MQ=255)
                                                                                             tGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGt                                                    >  1:462430/1‑150 (MQ=255)
                                                                                              gCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcacac                                                         <  1:208699/144‑1 (MQ=255)
                                                                                              gCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcacac                                                         >  2:208699/1‑144 (MQ=255)
                                                                                                tCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAAc                                                 >  1:867412/1‑150 (MQ=255)
                                                                                                tCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAAc                                                 <  2:1255530/150‑1 (MQ=255)
                                                                                                 cGGCCAGATCGCGCATCATGACGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACg                                                >  2:1192443/1‑150 (MQ=255)
                                                                                                    ccAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCt                                             <  2:429276/150‑1 (MQ=255)
                                                                                                       gATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTggc                                          >  2:696294/1‑150 (MQ=255)
                                                                                                          cgcgCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTggcggc                                       <  1:975265/150‑1 (MQ=255)
                                                                                                          cgcgCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTggcggc                                       <  1:498684/150‑1 (MQ=255)
                                                                                                            cgcATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTg                                     <  1:1261697/150‑1 (MQ=255)
                                                                                                              catcatGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTa                                   <  1:741285/150‑1 (MQ=255)
                                                                                                              catcatGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTa                                   >  1:1275146/1‑150 (MQ=255)
                                                                                                                tcatGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGC‑‑AACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCa                               <  1:1535605/150‑1 (MQ=255)
                                                                                                                  atGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCa                               >  1:1609767/1‑150 (MQ=255)
                                                                                                                   tGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaa                              >  1:1413821/1‑150 (MQ=255)
                                                                                                                        tGCGAGGTACAAAGCATAAAAGAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGATCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtggg                         <  1:1267070/150‑1 (MQ=255)
                                                                                                                         gCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGTTCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggc                        >  2:353730/1‑150 (MQ=255)
                                                                                                                          cGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggct                       >  2:81836/1‑150 (MQ=255)
                                                                                                                              gTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggctgggc                   >  2:1034215/1‑150 (MQ=255)
                                                                                                                                aCAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCaa                 >  2:20919/1‑150 (MQ=255)
                                                                                                                                 cAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAAc                <  2:345819/150‑1 (MQ=255)
                                                                                                                                 cAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAAc                <  1:1186632/150‑1 (MQ=255)
                                                                                                                                  aaaGCATAAAACAACGACAGTTGTTAGCTTCGATGATCGGTCGCAGCATTGACGCTAACTGGCGAGCAGAACATGGCTCGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACa               >  2:714564/1‑150 (MQ=255)
                                                                                                                                   aaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAg              >  1:405860/1‑150 (MQ=255)
                                                                                                                                   aaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAg              <  2:787859/150‑1 (MQ=255)
                                                                                                                                    aGCATAAAACACCGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGAGGCAGATTGCGCAGCACACAGAGTAACGCATGGCGGCTGTAATCAAATGGGCTGGGCAACAGc             >  1:1553135/1‑150 (MQ=255)
                                                                                                                                     gaagctaacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCa            >  2:855010/7‑150 (MQ=255)
                                                                                                                                     gCATAAAACAACGACCGTTGTTAGCTTCCATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCa            >  2:527400/1‑150 (MQ=255)
                                                                                                                                      cATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCaa           >  1:1087164/1‑150 (MQ=255)
                                                                                                                                      cATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCaa           >  1:1488789/1‑150 (MQ=255)
                                                                                                                                       agctaacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAAc          >  2:1227585/5‑150 (MQ=255)
                                                                                                                                        tAAAACAACGACCGTTGTTAGCTTCGATGCTCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGGTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACg         >  2:481431/1‑150 (MQ=255)
                                                                                                                                        gctaacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACg         >  1:464166/4‑150 (MQ=255)
                                                                                                                                          taacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcaca                                                          <  1:1080409/98‑1 (MQ=255)
                                                                                                                                          taacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcaca                                                          >  2:1080409/2‑99 (MQ=255)
                                                                                                                                          aaaCAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGAc       >  2:929169/1‑150 (MQ=255)
                                                                                                                                             caaggtccGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCGGGGTGGTTTGCGTCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCATCGACTCg    <  2:926398/144‑1 (MQ=255)
                                                                                                                                             caacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCg    >  2:67022/1‑150 (MQ=255)
                                                                                                                                              aacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCgg   <  1:327767/150‑1 (MQ=255)
                                                                                                                                              aacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCgg   >  2:643255/1‑150 (MQ=255)
                                                                                                                                               acGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGc  >  1:486087/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGC  >  NC_000913/1891695‑1891987

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: