Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,1910AC30.2% 41.2 / 20.2 53intergenic (+397/+245)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (22/15);  new base C (10/6);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTA  >  NC_000913/4295892‑4296475
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 
gATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACGTAAACGGCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTaa                                                                                                                                                                                                                                                                                                   <  2:905665/300‑1 (MQ=25)
 ttctttattGATTCGGTAGGTGTGATAAGTCTTTTATGCCGCTCCCGAATTCACCGCCTATTGCGACGCTTAACCCGTGTTATCAGGCTGGTGGCACACTGCGCAATAGCTTGATTTAGCTTGATTTTGTAGTTCGGATATGGCGTTTACGCCGCGTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTACGCCAGACAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAAc                                                                                                                                                                                                                                                                                                  <  2:1088222/292‑1 (MQ=11)
   ttAGGGTGATGTTGTAGGTCGGAGAAGGCATTTATGCGGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAAcgc                                                                                                                                                                                                                                                                                                <  2:598996/300‑1 (MQ=255)
       gttGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCATCCTACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTc                                                                                                                                                                                                                                                                                             <  2:576266/297‑1 (MQ=25)
       cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCCGATGCG‑ACTCTTACCGCGTCtt                                                                                                                                                                                                                                                                                           >  2:1089133/1‑301 (MQ=25)
              ttGGAGGGCGGATAAGGCTGTTATTCGGCGTCCGGCACCACCGACTGATGACAAGCTTACAGCGTCTTGTCAGGCCTACGCCAGACAGTGCAATATACTGGTTTAGCTGGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAgg                                                                                                                                                                                                                                                                                     <  2:77220/300‑1 (MQ=11)
                   gggcGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATAAG‑ACGCGTTAAGCGTCGCATCAGGcgttga                                                                                                                                                                                                                                                                               <  2:241012/298‑6 (MQ=255)
                    gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCCGACCGCGCAATAGCCTGATTTACCGTGATTTCGTAGTTCGGAAAATGCCTTCATGCCGCAGCCCCAATCAATCCTTAATGCG‑CCTCATAAGGCGTTTTCTCAGGCTTAAGc                                                                                                                                                                                                                                                                              >  1:610067/1‑301 (MQ=11)
                      gggATAATGCGTTTATGCCGCATCCGACATCAACGCTCGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGGGCAATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCa                                                                                                                                                                                                                                                                            <  2:533701/300‑1 (MQ=255)
                        gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTCAGCCGCATCCGCCATCAACACCGGTTGCG‑ACGCTTCACGAGTCTTATCACGCATCCGCTAgc                                                                                                                                                                                                                                                                          >  1:875985/1‑300 (MQ=11)
                             ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTATGCCGCATCCGCCATCAACGCCTGATGCG‑ACGCTTACCTCGTCTTATCAGGCCTACGCCACACCgcc                                                                                                                                                                                                                                                                     >  2:14569/1‑300 (MQ=255)
                                 gtgaggccgcatccgcactgcactgctggagCGAAGTTTATCGCGTCTTATCGTGGCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                                                                 <  2:267959/271‑1 (MQ=11)
                                                    cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCCCATCCCACATAAACGCCTGCTGCG‑ACGCTTGAGCCCTCCTCTCCAGCCTCTCCCCGACAGCGCAAAATCCTGGTTCTGCCTGAtt                                                                                                                                                                                                                                              >  1:127112/1‑301 (MQ=11)
                                                    cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTATGCCGCATCCGACATCAACTCCTGATGCG‑CCGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCCATAGCCTGATTAAGCGTGAtt                                                                                                                                                                                                                                              >  2:377778/1‑301 (MQ=255)
                                                      aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtttt                                                                                                                                                                                                                                            >  2:264640/1‑301 (MQ=34)
                                                         ccTGATAAGACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAccg                                                                                                                                                                                                                                                                     >  2:812303/1‑273 (MQ=255)
                                                                       tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGCACTATCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTATGCCGCATCCGACCTCAACGCCTGACGCG‑TCGCTTAACGCGTCTTTTCAAGCCTCCTCCCAACCGAGGAACAGCCTGCGTTAGCGTAAtttggccagtccgataggggg                                                                                                                                                                                                                           >  1:175158/1‑283 (MQ=11)
                                                                            cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTac                                                                                                                                                                                                                      >  2:435474/1‑300 (MQ=17)
                                                                                  aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATAAGGCCTACGCCAGACAGCCCACTAGCCTGCGTTACCGTGATTTTGTAGGTAGGATACGGCGTTTAG‑Gcctcc                                                                                                                                                                                                               >  2:670398/1‑298 (MQ=11)
                                                                                  aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGACTTAGTGGGTCGGATAAGGCGTTTAc‑gccgcc                                                                                                                                                                                                               >  1:424624/1‑300 (MQ=16)
                                                                                      ggcctaacccataaggcgaAATAGGTTGATTTAGCATGATTTTGTAGGGCGGAGAAGGCGTTTACGCCGCAACCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCg                                                                                                                                                                                                           <  2:737460/286‑1 (MQ=11)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGGG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCGACATc                                                                                                                                                                                                      >  2:1037494/1‑301 (MQ=25)
                                                                                           aCGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGGATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCGACa                                                                                                                                                                                                        <  2:268673/299‑1 (MQ=25)
                                                                                              ccgaacagcctaaactccTGTTTTCGCGTGATTTTTGAGTTAGATTAAGGGTTGTCCGCCGCGTCCGACCTCACCGCCTGATGGGACGCTTGCCGCGTCTGTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCGACATCaat                                                                                                                                                                                                   <  1:1107360/285‑2 (MQ=11)
                                                                                              ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTc                                                                                                                                                                                                                                                                                             >  2:783511/1‑212 (MQ=21)
                                                                                                 gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGCTGCG‑ACGCTTACCGCGTCTTATCCGTCCTAAGCCAAGCCGCGGTATAGCTTGATTTAGCGGGATTCTGTGGGTCGGTCAAGGC‑TGTATCGCCGCATCGGAccgcaagcgc                                                                                                                                                                                                >  1:73959/1‑292 (MQ=11)
                                                                                                 gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTGACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGTGTTTCTGATGCACGCGACATAAAGGCATGATGGC‑ACGTCTAACGCGTCTTATAACGACTACGACAGACCGCGCTATAGCCTGATTTAGAGTGATTTTGTAGGTCGTAGAGGGGGTTTCT‑GCAGGATTCGACATCacagcg                                                                                                                                                                                                >  2:649226/1‑296 (MQ=11)
                                                                                                     gcgcAATAGCTTGATTTAGCGTGATTTTGTAGGTCGCAGAAGGCTGTTACGCCGCATCCGACATCAACGCCTGATGGGAGGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCGACATCAATGCCTGAt                                                                                                                                                                                            <  2:926152/301‑1 (MQ=11)
                                                                                                     gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGCCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGTTAAGGCGTGTCC‑GCCGCATCCGCAATCAATGCCGGAt                                                                                                                                                                                            >  1:439927/1‑301 (MQ=11)
                                                                                                     gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTAACAGGCCTACGCCACACAGCCGCATAGCCTGTTTTAGCGTTCATGTGTCGATCGGATCAGGCGGTTAC‑GCCGCGTCCGACactcatgtctg                                                                                                                                                                                              >  1:362391/1‑289 (MQ=11)
                                                                                                       gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTAC‑GCCGCATCCGACATCACTTCCTGAttc                                                                                                                                                                                          >  2:454395/1‑299 (MQ=19)
                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCa                                                                                                                                                                                                                                                                   <  2:359069/223‑1 (MQ=33)
                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGACAGACAGCGCa                                                                                                                                                                                                                                                                   >  1:358989/1‑223 (MQ=18)
                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCC‑ACGCTCAACGCGTCCTATTAGGCCTACGCTAGACAGCGCCAATGCCCGTTTTAGCGCCCCTTCATAGGTCGGATCAGGAGTTTTC‑GCCGCTTCCGCTATCactgcatgtggccacgcc                                                                                                                                                                                    >  1:416713/1‑284 (MQ=11)
                                                                                                               cggATTTCGCGTGTTTTGCTAGGTCTGATAAGGCTTTTACGCCGCATCCGCAATCAACGCCTGATTCGGCGCTTGACGCGTCTTATCAGTCCTACGCCAGACAGCCCAATAGCCTGATTTAGCGTGATTTTCTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCGACATCAATGCCTGATGCGACGCttg                                                                                                                                                                                  <  2:105771/299‑2 (MQ=11)
                                                                                                                         gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGTGCAATAGCCTCATCTAGAGTGATTTTGTCGGTCGGACAAGGCGTTTAC‑GCCCGATCCCGCATCAATGCCCGATGCGACACttgccgctgcttc                                                                                                                                                                        >  1:11000/1‑290 (MQ=11)
                                                                                                                            aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTGAC‑GCCGCATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                      >  2:1118227/1‑300 (MQ=12)
                                                                                                                                tGTCGGTCGGGAGAGGCGTCTACGCCAGAACCGACAGCACCGCCGGATGCGATGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATGTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACTCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCGACATCAATGCCTGATAAGACGCGGCAAGCGTCGCATCAGGc                                                                                                                                                                  <  2:903389/300‑1 (MQ=11)
                                                                                                                                   aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAC‑GCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                                                              >  2:928604/1‑300 (MQ=14)
                                                                                                                                    ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑AAGCTTAACGCGTCTTATCAGGACGACGCCAGCCAGAGCAATAGACGGTCACAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCA‑GCGACATCCGACACTAATGACTGATGCGCCGAttgcagcgagctcggacgatcgtg                                                                                                                                                             >  1:92109/1‑279 (MQ=11)
                                                                                                                                                                                                                                           ggATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCACCAAAATAAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTTATTTAGCGTAATTTTGTAGGTAGGATAAGGCTTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                    <  2:460316/300‑1 (MQ=11)
                                                                                                                                                                                                                                                            tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                   >  2:700207/1‑301 (MQ=2)
                                                                                                                                                                                                                                                              aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTACCAGGCCTACGCCAGACAGCGCAGTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACCTCAATGCCTGATGCGACGCTTGCCGCGTCTTAGCAGGCCTATCTTAACCGtt                                 >  2:89296/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    ttACGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCCGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTCTCCGGCCTATCTTAACCGTTGGTTa                            >  1:531611/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                         ccGCTTCCTACACCACCCCCTG‑TGCGTACGCTTTCCGCGTCTGACCAGGCCGACGACAGAAAGCGTAATATCCTGATTTAGCGGGATTTTGTAGGTAGGATACGGCGTTTAT‑GCCGCATCCGACATCACCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCAACGGTTaagatag                      <  1:377694/301‑6 (MQ=11)
                                                                                                                                                                                                                                                                         ccGCATCCGATATTAACGCCGTATGCG‑ACGCTTGCCGCGTCTTATCCGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                      <  2:165208/301‑1 (MQ=11)
                                                                                                                                                                                                                                                                          cGCATCCGACCGCACTGCCGGAGGCG‑TCGCTTGCGTCGTCTTATCAGGCCTACACCAGCAATCGCAATAGCCTGATTTGACGTGATGTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                     <  1:700044/301‑1 (MQ=11)
                                                                                                                                                                                                                                                                           cccccccaCATCAACGCTTGATGCG‑AAGCTTGCCGTGTCTTATCAGGCCTACGCCAGACAGCGGAATAGCCTGATTTTGCTTGATTTTGTAGGTCGGATAAGGCGTTTAG‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                    <  2:531733/294‑1 (MQ=11)
                                                                                                                                                                                                                                                                                cGACATCAACGCCTGATGCG‑ACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCATGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACCTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACctc               >  2:927688/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                         cGCCTGAGTAG‑ACGCTTGCCGCGTCTTTTCAGGACTACGCCAGCAAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGACTACGCCAGACAGCGCAATAGCATGATTTAGCATGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc      <  1:927455/301‑1 (MQ=14)
                                                                                                                                                                                                                                                                                          gCCTGATGCG‑ACGCTTGCCGCGTCTTATAAGGCCTACGCCAGACAGCGCAATAGCCTGATCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATTCGACATTCACGCCTGATGTGACGCTTAACGCGTCTTATCAGGCTTACGCTAGACATCGCAATAGCCTGGTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTG‑‑TATGCCGGATCCGACAGAAATGTCCGATGCGTTGCTGGCAGCGTCTTATACGGCCTATCTTACTCGCTGGTTCATTTTTCGAACTTTTcaatcac      >  2:766453/1‑295 (MQ=11)
                                                                                                                                                                                                                                                                                          gCCTGATGCG‑ACGCGTGCCGCGTCTTATCAGGCCTACGCCATACAGCGCAATATCCTAATTTAGCGTGGTTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCATACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg     <  2:416802/301‑1 (MQ=14)
                                                                                                                                                                                                                                                                                             tGATGCG‑ACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATCGCCTGATTTAGCGTGATTTTGTAGGTCGGACAATGCGTTGCCG‑‑CCGTAACCGACATCAATGCCTGACGCGACGATTGCCGTGTCTTATCCGGCCTCTTTCACCCGTCGGTTCATTTTTCCAAACCTTctacccaggtt  >  1:684101/1‑291 (MQ=11)
                                                                                                                                                                                                                                                                                             tGATGCG‑ACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATTCTGCCTTTCTTAACCGTTGGTTAATTTTTTAAACTATTcaaccccagtt  >  2:826312/1‑292 (MQ=11)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 
GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCG‑ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT‑GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTA  >  NC_000913/4295892‑4296475

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑