Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,380:1 +CG 100% intergenic (+586/+56) 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.C95.2% 169.2 / 4.1 62intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (35/22);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0009134,296,3802.G95.2% 173.3 / 3.4 62intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base G (35/23);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTT  >  NC_000913/4296266‑4296520
                                                                                                                   ||                                                                                                                                               
tATGCAGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTAAGCCAGACACATCAATAGCCTTAATTAGACTCATTTTGTAATTCGGATAAAACGTTTAATC‑‑C‑TCATCAGAAATCAATGCATCATGctaacctt                                                                                                             >  2:201363/1‑143 (MQ=0)
                            gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGA‑TTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                                >  2:1435948/1‑150 (MQ=11)
                             cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                                >  2:1508343/1‑150 (MQ=9)
                              gACGCTTACCGCTTCTTATCTGGCCTACGCCAGACAGCGCACTACCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                               <  1:765317/150‑1 (MQ=11)
                                cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                             >  2:383030/1‑150 (MQ=255)
                                 gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                            >  2:222350/1‑150 (MQ=255)
                                    tAACGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGACTGATGCGACGCTTGCCGCGTATTATCAGGCCTATCTTAACCGTAGGTTa                                                                         >  1:1601725/1‑150 (MQ=11)
                                     aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                        >  1:887476/1‑150 (MQ=255)
                                     aaCGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGAAGGTCGGATAAGGCGTTTACGC‑‑C‑GCAACCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                        >  2:2062655/1‑150 (MQ=11)
                                        gcgTCTTATCATGCCTACGCCAGTCAGCGCAATAGCCTGATTTAGATTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                     <  2:58622/150‑1 (MQ=14)
                                        gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGTTTACGC‑‑C‑GCTGCCGACATCAATGCCTGATGCGACGATGGCCGCGTCTTATAAGGCGTAGCTTACCCGTTGGTTAAttt                                                                     >  1:468407/1‑150 (MQ=11)
                                         cgTCTTATCAGGCCTACTCCTTACATCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTTATAAGTCGTTTACGC‑‑C‑GAATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                                    <  2:1483808/150‑1 (MQ=14)
                                            cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                 <  1:596246/150‑1 (MQ=255)
                                             ttttCAGGCCTAATCCAGACAGAGCAATAGATTTATTTAGCGTGATTTTGTAGGTTGTATTAGGTGTTTAAGA‑‑C‑TCATCCGAAATCAATGCCTTTTGCGACTCTTGTCGCGTCTTATCATGCCTATCTTAACCGTTGGTTAATTTTTCat                                                                <  1:201363/147‑2 (MQ=9)
                                             ttATCTGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGCATACGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                <  1:1480108/150‑1 (MQ=14)
                                               aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                              >  1:1482898/1‑150 (MQ=255)
                                               aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                              >  2:609315/1‑150 (MQ=255)
                                                 cagctcTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                            <  1:404495/145‑1 (MQ=14)
                                                 cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                            >  1:512027/1‑150 (MQ=255)
                                                 cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                            >  2:1603743/1‑150 (MQ=255)
                                                  aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                           >  2:2096338/1‑150 (MQ=255)
                                                   ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                          >  1:1039923/1‑150 (MQ=255)
                                                      cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                       <  1:1207263/150‑1 (MQ=255)
                                                      cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                       <  2:1780103/150‑1 (MQ=255)
                                                           ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                  <  1:668944/150‑1 (MQ=255)
                                                            cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                 >  2:546769/1‑150 (MQ=255)
                                                             aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                                <  2:648505/150‑1 (MQ=255)
                                                                   cgcAATAGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCATGATGCGACGCTTGCCGCGTATTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                          >  1:107673/1‑150 (MQ=18)
                                                                    gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGTCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                         >  1:1468011/1‑150 (MQ=255)
                                                                     cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                        <  1:1744795/150‑1 (MQ=255)
                                                                     cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                        <  1:546548/150‑1 (MQ=255)
                                                                       aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCt                                      >  1:1778134/1‑150 (MQ=255)
                                                                        tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                     >  1:2076643/1‑150 (MQ=255)
                                                                        tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                     >  1:1888277/1‑150 (MQ=255)
                                                                         aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                    <  1:726700/150‑1 (MQ=255)
                                                                          gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                   >  1:733866/1‑150 (MQ=255)
                                                                          gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                   <  2:1871733/150‑1 (MQ=255)
                                                                          gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                   >  1:1370646/1‑150 (MQ=255)
                                                                           ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CG‑CAGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAACGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                  >  2:1297999/1‑150 (MQ=21)
                                                                           ccTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCAACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                  >  2:1504961/1‑150 (MQ=21)
                                                                           ccTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                  <  1:1037036/150‑1 (MQ=255)
                                                                              gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                               >  2:123128/1‑150 (MQ=255)
                                                                               aTTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CG‑CAGCATCCGACATCAATGCATGATGCGACGCTTGCCCCGACTTATCAGGCCTATCTTAACAGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAcc                              >  2:1079571/1‑149 (MQ=21)
                                                                                tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                             >  2:1498404/1‑150 (MQ=255)
                                                                                    gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                         >  1:849938/1‑150 (MQ=255)
                                                                                      gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                       >  1:2069695/1‑150 (MQ=255)
                                                                                      gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                       <  1:322566/150‑1 (MQ=255)
                                                                                       tGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                      <  1:433680/150‑1 (MQ=255)
                                                                                        gATGTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGAGTCTTATCAGGCCTATCTTAACAGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAAACCCTTCCATACAGCTCAGGtt                     <  1:121952/150‑1 (MQ=21)
                                                                                          ttttGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa                   >  1:240418/1‑150 (MQ=255)
                                                                                            ttGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                 <  1:61427/150‑1 (MQ=255)
                                                                                            ttGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                 <  1:1816293/150‑1 (MQ=255)
                                                                                                aGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGAAGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACTGCTCAGGTTCAACGAGt             >  1:1799618/1‑150 (MQ=25)
                                                                                                 ggTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg            >  2:155766/1‑150 (MQ=255)
                                                                                                 ggTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg            >  1:741712/1‑150 (MQ=255)
                                                                                                 ggTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg            >  1:1929648/1‑150 (MQ=255)
                                                                                                 ggTCGGATAAGGCGTTTA‑‑CGAC‑GCATCCGACATCAATGCCTGAAGCGACGCTTGCCGCGTCTTATCAGGACTATCTTAAAAGGTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCAGACAGATCAGGTTCAACCAGTg            >  1:1538032/1‑150 (MQ=21)
                                                                                                  gTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc           <  2:887476/150‑1 (MQ=255)
                                                                                                    cGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg         >  2:1021482/1‑150 (MQ=255)
                                                                                                        tAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc     <  1:103259/150‑1 (MQ=255)
                                                                                                         aaGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCACCCAGTGCAGCGCCt    >  2:529519/1‑150 (MQ=255)
                                                                                                           ggCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCttt  >  2:295975/1‑150 (MQ=255)
                                                                                                                   ||                                                                                                                                               
TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTT  >  NC_000913/4296266‑4296520

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