Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,256,0690GC36.6% 46.2 / 14.3 42P133R (CCG→CGG) plsBglycerol‑3‑phosphate 1‑O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/23);  new base C (15/0);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.29e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG‑TGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCA  >  NC_000913/4255836‑4256306
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
gCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCA‑CTTCGCCTTTTTTTCGCCCCCGCGCGCGACCAAAAATc                                                                                                                                                                                                                                 >  1:106224/1‑251 (MQ=255)
 cAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTAT‑CCTCGCCTTTTTCACCCCCCgccccccgccccacaatcc                                                                                                                                                                                                                                >  2:476813/1‑233 (MQ=255)
      gTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATCCGCGCCACGCTCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGTCGATAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGac                                                                                                                                                                                                                           <  1:320207/251‑1 (MQ=255)
      gTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGAGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCAACTT‑ACCTTTTTCTCACCCCCGAGCGCGACCAAACATcccccac                                                                                                                                                                                                                           >  1:518908/1‑245 (MQ=255)
          gCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTCTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑AATCGCCCTTTTTACGGCCCggcgggcggacccacaaccacccacacc                                                                                                                                                                                                                       >  2:398879/1‑231 (MQ=255)
                ggCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCa                                                                                                                                                                                                                <  1:123262/251‑1 (MQ=255)
                  cAGCCAGACGTTGACGGGCAAAGGGCATACGCTCCGCGCGCTCCTGTTTCGGAGCGATATTTTTACCCGGGCCTGGTTCATCCTCCATCCGGCGCAGCGAAACGGACGGCGAGAAACGCACAAACC‑TTTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCTCCGACACTGGC‑ACCATc                                                                                                                                                                                                              <  2:77346/251‑1 (MQ=255)
                       aGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCA‑CTCTGCCCTTTTTACCCCCCCGCGCGCGGCCCAACCACACCCACACTGGccccctactccc                                                                                                                                                                                                          >  1:371721/1‑240 (MQ=255)
                         aCGTTGTCGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCCTTTTTACGCCCCCGCGCGCGGACAAAAATCAACGacacccgtaccctctgccca                                                                                                                                                                                                         >  1:144366/1‑234 (MQ=255)
                                ggggaaaagcgtttctggccccggcgcgcCAGGTTGTGAGCGAGTGTTTTACCGGGGCCGTTTTCATCCGCCATACGGCGCAGCGAAACTGACGGCAAGAAACGCAAAAACC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCTTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTGTCCCGCCCCTGCGCGCGACCAAACATCACCGACACTGGC‑ACCAGCTGCACAGCCAGAtg                                                                                                                                                                                                <  2:489779/228‑2 (MQ=255)
                                          tcaatcgcgCAACGGGCGCCAGTTTCTGAGCGATAGTTTTACCCGTGCCGTTTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCTCAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACTCAGCGGCGGGTTCT‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTTGGTTGCt                                                                                                                                                                                       <  1:188943/245‑1 (MQ=255)
                                                          cgcCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTA‑CTTCGCCCTTTTCACGCCCCCGCGCGCGACCCAAAACTACCGACACCGGCCACC‑TCTTCACACCCAGATTTGGGTTGCTACAGGGCaaaacaa                                                                                                                                                                         >  2:507707/1‑244 (MQ=255)
                                                             cAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCCCGCCCCCGCGCGCGACCAAACATCACCCACACTGGC‑CCCATCCGCACCTCCCGATTTGGGTTGCTACGGTGCAAATC‑GCGATAGt                                                                                                                                                                  >  1:170580/1‑251 (MQ=255)
                                                             cAGTTTCGGAGCGATAGTTTTACCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATcc                                                                                                                                                                                                     <  1:441984/217‑1 (MQ=255)
                                                               gTTTCTGACCGATAGTTTTAACCGTGCCGTGTTCTCCCGCCAGACGGCGCAGCTAAACTGACGGCGAGAAACGCCCAAAAC‑GGTCGCGACCGAGCCACAGTACAGCGAAAAATTGCTGTACGCCGTTAAGCAGACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCg                                                                                                                                                                <  1:156314/251‑1 (MQ=255)
                                                                  tCTGAGCGATAGTTTTATCCGTGCCGTTTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACGGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTg                                                                                                                                                              <  1:43547/250‑1 (MQ=255)
                                                                     gagtaatgtttttATCCGTGCGGTTTCCCGCCGCCAGACGGGGAAGCAAAACTGACGGCGAAAAACG‑ATAAAACCTGTCGCGCCCGGGCCACAGTACAGCGAAAAAGTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGCCCAAACAGCACCGACACTGGC‑ACCATCTGCACAGCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCTTGGAAc                                                                                                                                                          <  2:236097/242‑1 (MQ=255)
                                                                            aGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCCGCGCGCGACCAAAAATCACCGACACTGGC‑ACCCTCTGGCCATCCAGATTTGGGGTGCTACGGTGCAAAAC‑GGGATAGTCGTGGGACAGCTTa                                                                                                                                                    >  1:235004/1‑250 (MQ=255)
                                                                              tttttcccgggccgtgTTCACCCGCCATACGGCGCAGCGAAACTGCCGGCGAGAAACGCACAAAAC‑TGTCGCTCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAAt                                                                                                                                                  <  2:405915/240‑1 (MQ=255)
                                                                                gaaccgagcccgtgTTCAGCCGCAATACGGCGCAGCAAAACTGACGGCAAGAAATGCAAAAAGC‑CGTCGCGTCCGAGCCACGGTACAGCGAAAAATGTCTGTACGCCGTTAAGCCTACGCGGCGGTGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGGGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAg                                                                                                                                                <  2:158267/241‑1 (MQ=255)
                                                                                   tCCGGGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGAct                                                                                                                                             <  1:432812/250‑1 (MQ=255)
                                                                                                 acccGCCATACGGCGCAGCGAAACTGACGGCGAAAAACGCACAAACC‑TGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTCCGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGGCGTGGAACAGCTTAATAGACTCTTCTTTCGGCg‑tgt                                                                                                                             <  2:288020/249‑1 (MQ=255)
                                                                                                    cGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTA‑CTTCTCCCTTTTTACGCCCCCGCGCGCGGCCAAACAACACCGACACCGGC‑ACCCTCTGCACAACCCGGTTTGGGGTGCTACGGTGCAAAAA‑GAGATAGTCGTGGAAAAGCTTAATAGACTCTTTTTTTCGCG‑TGTaa                                                                                                                           >  1:421931/1‑250 (MQ=255)
                                                                                                                   agAAACTGACGGCGAAAAACGCACAAAAC‑TGCCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGTCGCGCGACCAAACATCTCCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG‑TGTAATAGGTGAACACAcgcg                                                                                                           <  2:144366/250‑1 (MQ=255)
                                                                                                                          gACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACAACAACGCCACTGGC‑ACCATCTTCACCTCCCGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGACTCTTTTTTTGGCG‑GGTAATAGGGGAACACACGCGGCccgcc                                                                                                    >  1:474257/1‑251 (MQ=255)
                                                                                                                              gCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTA‑CTTCTCCTTTTTTACGCCCCCGCGCGCGGCCAAACATCACCGACACTGGC‑CCCCTCTGCACATCCAGATTTTGGTTGCTTCGGGGCCAAAC‑GAGGTAGTCGTGGGACAGCTTAATAGAATCTTTTTTCTGCG‑GGGAATAGGTGGACACACGCGGCCCGCCGtgg                                                                                                >  1:2524/1‑250 (MQ=255)
                                                                                                                                gagaAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCCGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTTCACAACCAGATTTGGGGTGCTACGGTGCAAAAA‑GAGGAAGTCGGGGGACAGGTTAAAAGAACCTTCTTTTTGCG‑GGTAATAAGGGAACACAAGCGGCCCCCCGtgaa                                                                                               >  1:292277/1‑250 (MQ=255)
                                                                                                                                      cGCAAAAAAC‑GGTCGCGCCCGAGCCACAGTAAAGCGAAAAATTTCTTTCGGCCGTTAAGCATACGCAGCGGCGGTTTCA‑CTGCGCCTTTTCCCCGCCCCGGCGCGCGACAAAACTCCCCCGACACTGGC‑ACCATCTGCACATCCAGATTGGGGTGGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG‑TGTAATTGGTGAACACACGCGGCCCGCCGTGAATGAacac                                                                                        <  1:268263/251‑1 (MQ=255)
                                                                                                                                           aaaaccttccgcgACCGACCCACAGCACAGCAAAAAATTTTTGTCCCCCGTCAAGCATCGGCGGCGGGGGTTTCA‑CGTCCCCTTTTCCCCGCCCCGCCGGCCCACAAAACACCCCCCACACGGGC‑CCCATCTGCACATCCATATTTGGGTTGTTACGGTGAAAATG‑GAGAAAGTCGTGCAACAGCTTAATAGACTCTTCTTCGGGCG‑GGTAATGGGTGAACCCACGCGGCCCGCCGTGAATGAACACATAg                                                                                    <  2:422933/242‑1 (MQ=255)
                                                                                                                                               ccttgcgggccccacccACAGTACAGCAAAAAATTCCTGTCCGCCGTTAAGCATCCGCGGCGGCGGTTTCA‑CGCCCCCTTTTCCCCGCCCCGGCGCGCAACAAAACACCCCCTACACTGGC‑ACCATCTGCACACCCAAATTTGGGTTGGTACGTTGCAAATC‑GTGATAGGCTTGGAACAGCTTAATAGACTCTT‑TTTCCTCGGTGTAATAGGGGAACACACGCGGCCCGCCGTGAATGAACACATAGcgcg                                                                                <  1:78648/236‑1 (MQ=255)
                                                                                                                                                       gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTTCCGCCCCCGCGCGCGACCa                                                                                                                                                                                                                                       >  1:3065/1‑95 (MQ=255)
                                                                                                                                                       gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCa                                                                                                                                                                                                                                       <  2:3065/95‑1 (MQ=255)
                                                                                                                                                          cGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAg                                                                                                                                                        <  2:471812/169‑1 (MQ=255)
                                                                                                                                                          cGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCCGCGCGCGACCAAACATCACCGACACTGGC‑CCCATCTGCACATCCAGATTTGGGTTGCTACGGGGCAAATC‑GAGATAGTCGTGGGACAg                                                                                                                                                        >  1:471812/1‑169 (MQ=255)
                                                                                                                                                            aGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCAt                                                                                                                                                                                                               <  1:215698/113‑1 (MQ=255)
                                                                                                                                                                            aaaaTTCCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTA‑CTTCGCCTTTTTCCCGCCCCGGCGCGCGGCCAAACATCACCGCCACTGGC‑ACCCTCTCCACATCCAGATTTGGGGTGCTACGGGGCAAATC‑GAGATAGTCGTGGGACAGCTTAATAGACTCTTCTTTTCGCG‑TGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGGAGCGTGCCGTCGATTTCCAGCGGc                                                   >  1:145389/1‑251 (MQ=255)
                                                                                                                                                                                  tCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCCCCTTTTTTCCGCCCCCGCGCGCGACCCAACAACACCCACAATGGC‑ACCACCTGCCCATCCAGATTTGGGTTGCTACGGTGCACAAC‑GAGATAATCGTGGGACAGGTTAATAGACCCTTCTTTTGGCG‑TGTAATAAGGGAACACACGCGGCCCCCCCTGGATGAACAAATAGCGCGGCAGTAGGGGGCCGGCGCTTTCCCGCGGCCCTAAc                                             >  1:454870/1‑251 (MQ=255)
                                                                                                                                                                                     gTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCCCGCCCCCGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTTGGGTGCTACGGTGCAAATC‑GAGATAGTCGGGGAACAGCCTAATAGACTCTTCTTTTGGCG‑TGTAATAGGGGAAAACCCGCGGCCCGCCGGGAATGAACACATAACGCGGCAGGAGCGGGCCCGCGATTTCCAGCGGCTCTAAAggg                                          >  1:238933/1‑251 (MQ=255)
                                                                                                                                                                                           cGTTAAGCATACGCTGCGGTGGTTTCT‑CTTCGCTTTTTTCCCGCCCCGGCGCGCCACCAAACACCACCGACACTGGC‑CCCATCTGCCCATCCAGATTTGGTTTGCTACGGTGCAAATC‑GGGATAGTCGTGGATCAGCTTAATAGACTCTTTTTTGGGCG‑TGTAAGAGGTGAACACACGCGTCCCTCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGc                                    <  1:473728/251‑1 (MQ=255)
                                                                                                                                                                                                          gcggcgGGTTCC‑CTCCGCCTTTTCCCCGCCCCGGCGCGCGACCAAACACCCCCGACCCTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTAAATAGACTCTTCTTTCGGCG‑TGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGa                     <  1:398879/251‑1 (MQ=255)
                                                                                                                                                                                                               gggTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG‑TGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACCCTg                >  1:444020/1‑251 (MQ=255)
                                                                                                                                                                                                                ggTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGCTTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG‑TGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTgg               >  2:217152/1‑251 (MQ=255)
                                                                                                                                                                                                                           cccttttcccccccccGGCGCGCGACCAAACAACCCCGACACTGGC‑CCCATCTCCCCATCCAGTTTTGGTTTGCTACGGT‑CAAAACAGAGATAGTCGTGAAACAGCTTAATAACCTCTTCTTCCGGCG‑TGTAATAGGTGACCCCCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACg     <  2:235004/239‑1 (MQ=255)
                                                                                                                                                                                                                               tttccACGCCCCGGCGCGCGTCCAACCACCCCCGACCCTGGC‑ACCATATGCACATCCAGATTTGGGTTGCTAGGGTGCAATGC‑GAGATAGTCTTGGACCACTTTAATAGACCCTTCTTGCCGCG‑TGTATTGGGTAACCGCACGCGCTCCGCCGTGAATTAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCa  <  1:73065/246‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
GCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAAC‑TGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCA‑CTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGC‑ACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATC‑GAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG‑TGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCA  >  NC_000913/4255836‑4256306

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