Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C98.8% 272.2 / ‑5.1 84intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (39/43);  minor base T (1/0);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
*NC_0009134,296,3802.G98.8% 279.1 / NA 84intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base G (40/42);  minor base C (0/1);  total (41/43)
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

GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913/4296294‑4296519
                                                                                         ||                                                                                                                                                  
gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                <  1:40688/150‑1 (MQ=9)
 cGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                               <  1:267242/150‑1 (MQ=255)
   aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCTCGTCTTATCAGGCCTATCTTAACCGtt                                                                             <  1:1021761/150‑1 (MQ=11)
    cGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                            <  3:36411/150‑1 (MQ=11)
    cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                            <  3:3210/150‑1 (MQ=255)
    cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                            <  3:3264/150‑1 (MQ=255)
     gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAAAGGGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                           >  1:25179/1‑150 (MQ=11)
       ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                         >  4:4250/1‑150 (MQ=255)
         aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                       >  3:183409/1‑150 (MQ=255)
          aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTCAt                                                                      >  2:1058871/1‑150 (MQ=255)
           cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                     <  1:1335443/150‑1 (MQ=255)
           cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                     <  1:136276/150‑1 (MQ=255)
            gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGACGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGACGCTA‑‑‑‑‑‑‑CGCTTGCCGTGTCTTATCAGGCCTATATTAACCGTTGGTTAAttt                                                                    >  4:319130/1‑150 (MQ=11)
            gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                    >  4:398856/1‑150 (MQ=255)
            gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                    >  2:789672/1‑150 (MQ=255)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                  >  3:325861/1‑150 (MQ=255)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                  >  2:40429/1‑150 (MQ=255)
               tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                 >  1:1114997/1‑150 (MQ=255)
                cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                >  4:558052/1‑150 (MQ=255)
                 ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                               >  1:45213/1‑150 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                              <  3:476841/150‑1 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                              <  1:579490/150‑1 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                              <  4:514636/150‑1 (MQ=255)
                   aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                             >  2:365130/1‑150 (MQ=255)
                   aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGCGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                             >  2:1288705/1‑150 (MQ=255)
                    tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                            <  3:14311/150‑1 (MQ=255)
                    tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                            >  2:1314976/1‑150 (MQ=255)
                     cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                           <  1:1314976/150‑1 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ATGCCGCATCCGACATCAACGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTCGGTTAATTCTTCAATctct                                                          >  2:161630/1‑150 (MQ=11)
                       ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  2:712123/150‑1 (MQ=255)
                       ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGACTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  3:319130/150‑1 (MQ=255)
                       ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGACGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  1:367538/150‑1 (MQ=255)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                        >  3:21709/1‑150 (MQ=255)
                        gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                        <  4:69760/150‑1 (MQ=255)
                          cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                      >  1:1112858/1‑150 (MQ=255)
                           tACTCCAGACAGCGCAACGGCCTCATCTAGCGTGATCTTGTAGGTCGGCTCACTCTCTC‑‑ACCCCGCCTCCCACACCAACGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                     <  2:980752/150‑1 (MQ=12)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATTCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                  >  2:978381/1‑150 (MQ=255)
                               ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccag                                                 <  1:736712/150‑3 (MQ=255)
                               ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 <  1:955356/150‑1 (MQ=255)
                                cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                >  2:1012899/1‑150 (MQ=255)
                                cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                >  3:106677/1‑150 (MQ=255)
                                  gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                              <  1:1049274/150‑1 (MQ=255)
                                  gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCCTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                              >  2:616936/1‑150 (MQ=255)
                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  4:136862/150‑1 (MQ=255)
                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  2:1338731/150‑1 (MQ=255)
                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         >  2:972877/1‑150 (MQ=255)
                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  2:184448/150‑1 (MQ=255)
                                        gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                        <  1:334374/150‑1 (MQ=255)
                                          aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                             >  4:109488/1‑150 (MQ=17)
                                          aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAAGCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                      >  4:548121/1‑150 (MQ=255)
                                          aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                      >  1:33827/1‑150 (MQ=255)
                                                 tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  1:338599/150‑1 (MQ=255)
                                                 tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  4:21709/150‑1 (MQ=255)
                                                  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                              <  2:47568/150‑1 (MQ=255)
                                                  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                              >  2:1362811/1‑150 (MQ=255)
                                                     ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATATTAAACGTTGGCTCATTTATCATATTATTCTCAACCGGAATCTAACCCTTCCATACAGc                           >  3:562102/1‑150 (MQ=14)
                                                        gCGTGACTTTGTAGGTCGCATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGACTATCTGTACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACCGCTCa                        <  4:199819/150‑1 (MQ=21)
                                                         cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       >  1:781128/1‑150 (MQ=255)
                                                         cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       <  4:439964/150‑1 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                   >  2:1021140/1‑136 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                   <  1:1021140/136‑1 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                     >  1:1091186/1‑150 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                     >  3:141153/1‑150 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                     <  1:161630/150‑1 (MQ=255)
                                                             aTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc                   >  2:144703/1‑150 (MQ=255)
                                                              ttttGTAGGTCGGATAAGGCGTTTCCG‑‑CCGCATCCGACATCGATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCACGGTGCCAGACCCTTCCATACAGCTCAGGTCCa                  >  2:616055/1‑150 (MQ=18)
                                                                ttGTAGGTCGGATAAGGCGTCT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGC‑‑‑‑‑‑‑CGCTTGCCCCGTCCTATCAGGCACATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                <  1:816114/150‑1 (MQ=255)
                                                                   tAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg             <  3:504438/150‑1 (MQ=255)
                                                                   tAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg             <  1:932059/150‑1 (MQ=255)
                                                                    aGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCCGTATCCAACCGTTCCATACCGCACAGGTTCAACCaat            >  1:736739/1‑148 (MQ=21)
                                                                    aGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt            <  1:307294/150‑1 (MQ=255)
                                                                    aGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt            <  1:40429/150‑1 (MQ=255)
                                                                      gTCGGATAAGGCGTTT‑‑ACGCCGCATNCGACATCAATGCNTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCATAGCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc          >  3:499639/1‑150 (MQ=37)
                                                                      gTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc          >  1:988454/1‑150 (MQ=255)
                                                                      gTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc          >  3:432716/1‑150 (MQ=255)
                                                                      gTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc          >  1:168748/1‑150 (MQ=255)
                                                                      gTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc          >  1:568438/1‑150 (MQ=255)
                                                                       tCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACTCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa         >  2:576316/1‑150 (MQ=255)
                                                                       tCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa         <  1:789672/150‑1 (MQ=255)
                                                                         ggATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc       <  3:185340/150‑1 (MQ=255)
                                                                         ggATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc       <  2:984232/150‑1 (MQ=255)
                                                                           aTAAGGCGTTT‑‑ACGCCGCCTCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcgc     <  2:916081/150‑1 (MQ=255)
                                                                           aTAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcgc     <  1:1310539/150‑1 (MQ=255)
                                                                              aGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt  <  1:1061993/150‑1 (MQ=255)
                                                                                         ||                                                                                                                                                  
GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGA‑‑‑‑‑‑‑CGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913/4296294‑4296519

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