Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,267,489 A→G 6.0% E57G (GAG→GGG)  ychA → transglutaminase‑like/TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,267,4890AG 6.0% 203.9 / 2.8 68E57G (GAG→GGG) ychAtransglutaminase‑like/TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base G (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

CGAATTTAATAAAGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCGGTTTTATTA  >  NC_000913/1267340‑1267637
                                                                                                                                                     |                                                                                                                                                    
cGAATTTAATAAAGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTagaaga                                                                                                                                                      <  1:160907/150‑1 (MQ=255)
 gAATTTAATAAAGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAagag                                                                                                                                                     >  1:180419/1‑150 (MQ=255)
 gAATTTAATAAAGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAagag                                                                                                                                                     <  1:1072598/150‑1 (MQ=255)
     ttAATAAAGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGt                                                                                                                                                 >  1:268282/1‑150 (MQ=255)
           aaGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGaaa                                                                                                                                           >  2:832094/1‑150 (MQ=255)
           aaGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGaaa                                                                                                                                           <  1:1062173/150‑1 (MQ=255)
             gcgcCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGaaaaa                                                                                                                                         >  1:1281101/1‑150 (MQ=255)
              cgcCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAAc                                                                                                                                        <  2:1403024/150‑1 (MQ=255)
               gcCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACt                                                                                                                                       <  1:38250/150‑1 (MQ=255)
                ccATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTg                                                                                                                                      <  2:1344640/150‑1 (MQ=255)
                 cATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGa                                                                                                                                     >  1:755117/1‑150 (MQ=255)
                 cATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGa                                                                                                                                     >  1:760800/1‑150 (MQ=255)
                   ttGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATc                                                                                                                                   <  2:400783/150‑1 (MQ=255)
                         gAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTg                                                                                                                             >  2:615997/1‑150 (MQ=255)
                         gAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTg                                                                                                                             >  2:1403223/1‑150 (MQ=255)
                          aaGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGt                                                                                                                            <  1:615997/150‑1 (MQ=255)
                              cATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTa                                                                                                                        >  1:1414864/1‑150 (MQ=255)
                                tGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGaaaa                                                                                                                                          <  1:74231/130‑1 (MQ=255)
                                tGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGaaaa                                                                                                                                          >  2:74231/1‑130 (MQ=255)
                                     cTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGAc                                                                                                                 >  2:702046/1‑150 (MQ=255)
                                       ggCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGc                                                                                                                    >  2:323908/1‑145 (MQ=255)
                                       ggCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGc                                                                                                                    <  1:323908/145‑1 (MQ=255)
                                       ggCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTg                                                                                                               >  2:515138/1‑150 (MQ=255)
                                           tGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGGGCAGTTGGAAAAACt                                                                                                                                       >  1:1383113/1‑122 (MQ=255)
                                           tGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGGGCAGTTGGAAAAACt                                                                                                                                       >  1:1377657/1‑122 (MQ=255)
                                           tGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGGGCAGTTGGAAAAACt                                                                                                                                       <  2:1377657/122‑1 (MQ=255)
                                           tGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGGGCAGTTGGAAAAACt                                                                                                                                       <  2:1383113/122‑1 (MQ=255)
                                            gCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGAt                                                                                                          >  1:1033409/1‑150 (MQ=255)
                                              gAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGAttt                                                                                                        >  1:322480/1‑150 (MQ=255)
                                               aaGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTa                                                                                                       >  1:226036/1‑150 (MQ=255)
                                                 gCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGCGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTaaa                                                                                                     >  2:1243993/1‑150 (MQ=255)
                                                  cAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAg                                                                                                    <  1:422293/150‑1 (MQ=255)
                                                   aaTCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGc                                                                                                   >  1:592894/1‑150 (MQ=255)
                                                       cgccgcGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCa                                                                                               >  2:465302/1‑150 (MQ=255)
                                                         ccgcGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCAcg                                                                                             <  1:928/150‑1 (MQ=255)
                                                          cgcgATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCAcgc                                                                                            >  2:416747/1‑150 (MQ=255)
                                                          cgcgATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCAcgc                                                                                            <  1:1180475/150‑1 (MQ=255)
                                                           gcgATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCAcgcg                                                                                           >  1:60285/1‑150 (MQ=255)
                                                             gATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGgt                                                                                         <  2:1197797/150‑1 (MQ=255)
                                                              tttttCGCTCGCAAGATGTTTACGACGTACTGTAGGGTCTCGTTAGTCTGGCGAAGGAAGATATCAGCCAGCTTCTGCCTTTAGAAGAGGAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGTCTGGGGATTTAAAGCCTCTCGCGgtg                                                                                        <  2:727407/149‑1 (MQ=255)
                                                                        gCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCt                                                                              >  1:1044514/1‑150 (MQ=255)
                                                                        gCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCt                                                                              <  2:180419/150‑1 (MQ=255)
                                                                               gTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGAt                                                                       <  1:1403223/150‑1 (MQ=255)
                                                                                tttACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATg                                                                      >  2:1494246/1‑150 (MQ=255)
                                                                                   acgccgAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCAt                                                                   >  1:169221/1‑150 (MQ=255)
                                                                                   acgacgAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGGTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCAt                                                                   <  2:1281101/150‑1 (MQ=255)
                                                                                             ggAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGa                                                         <  1:442413/150‑1 (MQ=255)
                                                                                             ggAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGa                                                         >  1:212290/1‑150 (MQ=255)
                                                                                               aGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGAcc                                                       <  1:1255422/150‑1 (MQ=255)
                                                                                                  gTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAgg                                                    <  2:268282/150‑1 (MQ=255)
                                                                                                  gTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAgg                                                    >  1:173841/1‑150 (MQ=255)
                                                                                                   tctcGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGgt                                                   <  1:898105/150‑1 (MQ=255)
                                                                                                   tctcGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGgt                                                   >  2:346106/1‑150 (MQ=255)
                                                                                                    ctcGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGgtg                                                  <  1:572265/150‑1 (MQ=255)
                                                                                                     tcGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGgtgt                                                 <  1:519724/150‑1 (MQ=255)
                                                                                                      cGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGtt                                                >  1:117716/1‑150 (MQ=255)
                                                                                                         tAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTaaa                                             <  2:1477644/150‑1 (MQ=255)
                                                                                                          aGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAg                                            <  1:603154/150‑1 (MQ=255)
                                                                                                          aGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAg                                            <  1:1230703/150‑1 (MQ=255)
                                                                                                          aGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAg                                            <  2:418336/150‑1 (MQ=255)
                                                                                                          aGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGACTCACGCTGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAg                                            <  1:598055/150‑1 (MQ=255)
                                                                                                             cTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAAt                                         >  2:367592/1‑150 (MQ=255)
                                                                                                              tGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATc                                        >  1:246447/1‑150 (MQ=255)
                                                                                                                  gaaggaagAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACa                                    <  2:760800/150‑1 (MQ=255)
                                                                                                                  gaaggaagAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACa                                    <  2:755117/150‑1 (MQ=255)
                                                                                                                      gaagaaATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGc                                <  1:426320/150‑1 (MQ=255)
                                                                                                                               agccagcTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTa                       <  1:702046/150‑1 (MQ=255)
                                                                                                                                     cTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTa                 >  1:764330/1‑150 (MQ=255)
                                                                                                                                           ccTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCg           >  1:1027749/1‑150 (MQ=255)
                                                                                                                                            cTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCgg          >  2:847636/1‑150 (MQ=255)
                                                                                                                                             tttAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCGGt         <  1:367592/150‑1 (MQ=255)
                                                                                                                                                 gaagaGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCGGTTtta     >  1:762692/1‑150 (MQ=255)
                                                                                                                                                  aagaGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCGGTTttat    <  1:1425426/150‑1 (MQ=255)
                                                                                                                                                   agagCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCCGTTttatt   >  1:678790/1‑150 (MQ=255)
                                                                                                                                                    gagCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCGGTTttatta  <  2:212290/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAATTTAATAAAGCGCCATTGTGCGAAGGCATGATCCTGGCTTGCGAAGCAATCCGCCGCGATTTTCCCTCGCAAGATGTTTACGACGAACTGGAGCGTCTCGTTAGTCTGGCGAAGGAAGAAATCAGCCAGCTTCTGCCTTTAGAAGAGCAGTTGGAAAAACTGATCGCGCTGTTTTACGGCGACTGGGGATTTAAAGCCTCACGCGGTGTTTATCGTCTTTCCGATGCATTATGGCTGGACCAGGTGTTAAAGAATCGACAGGGCAGTGCGGTATCATTAGGTGCGGTTTTATTA  >  NC_000913/1267340‑1267637

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: