Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,517,571 C→T 5.3% G306R (GGG→AGG)  aroB ← 3‑dehydroquinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,517,5710CT 5.3% 254.1 / 3.1 76G306R (GGG→AGG) aroB3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base T (2/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

TTAAGAACAAGCTCGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTAACCAATTGCC  >  NC_000913/3517426‑3517709
                                                                                                                                                 |                                                                                                                                          
ttAAGAACAAGCTCGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCAt                                                                                                                                          <  4:58403/148‑1 (MQ=255)
  aaGAACAAGCTCGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCAtt                                                                                                                                         >  2:266676/1‑147 (MQ=255)
     aaCAAGCTCGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGAcc                                                                                                                                     >  4:241491/1‑148 (MQ=255)
           cTCGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAAc                                                                                                                               <  8:169285/148‑1 (MQ=255)
             cGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGccc                                                                                                                                            <  8:79154/133‑1 (MQ=255)
             cGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGccc                                                                                                                                            >  7:79154/1‑133 (MQ=255)
               tGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAg                                                                                                                           >  1:224797/1‑148 (MQ=255)
                   aaaCGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTa                                                                                                                                 >  1:212924/1‑138 (MQ=255)
                   aaaCGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTa                                                                                                                                 <  2:212924/138‑1 (MQ=255)
                   aaaCGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAAc                                                                                                                               >  4:195385/1‑140 (MQ=255)
                   aaaCGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAAc                                                                                                                               <  3:195385/140‑1 (MQ=255)
                   aaaCGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCg                                                                                                                       >  2:119779/1‑148 (MQ=255)
                    aaCGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGc                                                                                                                      <  2:248512/148‑1 (MQ=255)
                      cgccgcTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGa                                                                                                                                       >  4:53592/1‑129 (MQ=255)
                      cgccgcTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGa                                                                                                                                       <  3:53592/129‑1 (MQ=255)
                          gcTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGc                                                                                                                <  1:246958/148‑1 (MQ=255)
                               gAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGt                                                                                                           <  1:119779/148‑1 (MQ=255)
                                  cTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCTCGGCCCATTGTCCGGTAACCCAGCCCGCTTGAGCAGGGTtat                                                                                                        <  7:143125/148‑1 (MQ=255)
                                    tCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTtata                                                                                                       >  4:260286/1‑147 (MQ=255)
                                      aCTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTataata                                                                                                    >  4:17208/1‑148 (MQ=255)
                                       ctctTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAg                                                                                                              >  1:253583/1‑137 (MQ=255)
                                       ctctTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAg                                                                                                              <  2:253583/137‑1 (MQ=255)
                                               aaTTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGtt                                                                                           <  6:19303/148‑1 (MQ=255)
                                                aTTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGttt                                                                                          <  8:95234/148‑1 (MQ=255)
                                                   gCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCg                                                                                        <  1:266676/147‑1 (MQ=255)
                                                    ccAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTGGAGCAGGGTTATACTACGCTGCGGTTCGGc                                                                                      <  7:88729/148‑1 (MQ=255)
                                                       aCGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGa                                                                                   <  3:17208/148‑1 (MQ=255)
                                                         ggAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAAc                                                                                 >  7:71263/1‑148 (MQ=255)
                                                            agaaTTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGAGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTa                                                                               >  3:225847/1‑147 (MQ=255)
                                                                    gcgcATCTGTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGAGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCTGCCCGCTTGAGCAGGGTTATAATACGCTGCTTTTCGGAAGAACTAAACTGCCCg                                                                      <  6:203447/148‑1 (MQ=255)
                                                                      gcATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCgag                                                                    <  3:260286/148‑1 (MQ=255)
                                                                       cATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCgaga                                                                   <  3:257432/148‑1 (MQ=255)
                                                                        aTCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGAc                                                                  >  8:23222/1‑148 (MQ=255)
                                                                          ctctcCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGt                                                                >  7:293050/1‑148 (MQ=255)
                                                                          ctcacCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGt                                                                >  5:6647/1‑148 (MQ=255)
                                                                           tctcCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGCCAGAACTAAACTGCCCGACACGtt                                                               >  3:165820/1‑148 (MQ=255)
                                                                                cGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGAc                                                          >  7:255355/1‑148 (MQ=255)
                                                                                    aGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTcc                                                      >  6:181664/1‑148 (MQ=255)
                                                                                       aCTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCcgcg                                                   >  8:151819/1‑148 (MQ=255)
                                                                                       aCTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCcgcg                                                   <  4:276291/148‑1 (MQ=255)
                                                                                          ttCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGagccg                                                >  5:176549/1‑148 (MQ=255)
                                                                                          ttCTAGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACAATTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGATCCTACGTCCGcgccg                                                >  6:192216/1‑148 (MQ=255)
                                                                                           tCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGcgccg                                                <  5:181664/147‑1 (MQ=255)
                                                                                             ttGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCTATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGcgccg                                                <  6:218455/145‑1 (MQ=255)
                                                                                             ttGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCTATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGcgccg                                                <  8:162084/145‑1 (MQ=255)
                                                                                             ttGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCTATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGcgccg                                                >  5:218455/1‑145 (MQ=255)
                                                                                             ttGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCTATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGcgccg                                                >  7:162084/1‑145 (MQ=255)
                                                                                               gTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCgccgcc                                              <  8:293050/145‑1 (MQ=255)
                                                                                                tCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCg                                                                      <  2:233709/120‑1 (MQ=255)
                                                                                                tCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCg                                                                      >  1:233709/1‑120 (MQ=255)
                                                                                                  aCGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCTGCGCCGCCATCAcc                                        <  6:6647/148‑1 (MQ=255)
                                                                                                  aCGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCAcc                                        <  2:224797/148‑1 (MQ=255)
                                                                                                  aCGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCa                                       >  5:209301/1‑148 (MQ=255)
                                                                                                   cGCAGCATATGCGGTAAATACGCCTGCGCGGCCATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCa                                       >  1:42759/1‑148 (MQ=255)
                                                                                                   cGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCAcc                                        >  3:64972/1‑147 (MQ=255)
                                                                                                    gcagcaTATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCa                                       >  5:293682/1‑147 (MQ=255)
                                                                                                     cagcaTATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCGTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATa                                     <  3:241491/148‑1 (MQ=255)
                                                                                                           atGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCgag                                                                    >  1:79074/1‑111 (MQ=255)
                                                                                                           atGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCgag                                                                    <  2:79074/111‑1 (MQ=255)
                                                                                                           atGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCg                                 >  8:147844/1‑146 (MQ=255)
                                                                                                           atGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCa                               >  4:41701/1‑148 (MQ=255)
                                                                                                             ggggTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGCCGTTCCGACGTCCGCGCCGCCATCACCATAc                                    >  1:266516/3‑141 (MQ=255)
                                                                                                             ggggTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATAc                                    <  2:266516/139‑1 (MQ=255)
                                                                                                               ggTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGa                           >  1:88868/1‑148 (MQ=255)
                                                                                                               ggTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGa                           >  7:160479/1‑148 (MQ=255)
                                                                                                                     tACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTataat                                                                                                     <  7:198345/68‑1 (MQ=255)
                                                                                                                     tACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTataat                                                                                                     >  8:198345/1‑68 (MQ=255)
                                                                                                                       cGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCAGCGCCGCCATCACCATACCCGCAGCGACCGCTTc                    <  4:165820/147‑1 (MQ=255)
                                                                                                                          cTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGAGCCGCCATCACCATACCCGCAGCGACCGCTTCACCa                <  6:176549/148‑1 (MQ=255)
                                                                                                                           tGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCTCTTGAGCAGGGTTATAATACGCTGCCTTTCGTCAGAACTAAACTGCACGAGACGTTCCGAAGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCAt               >  3:208495/1‑148 (MQ=255)
                                                                                                                            gcgcGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATg              >  7:299110/1‑148 (MQ=255)
                                                                                                                             cgcgGACATTTCGCGCGGCCCATTGAACGGTATCCCTGCCCGCGTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCTCAGCGACCTCTTCACCAt               <  8:233687/146‑1 (MQ=255)
                                                                                                                              gcgGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTa            <  4:64972/148‑1 (MQ=255)
                                                                                                                                        cgcgcgGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTAACCAATTGc   <  7:23222/147‑1 (MQ=255)
                                                                                                                                        cgcgcgGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTAACCAATTGcc  <  6:209301/147‑1 (MQ=255)
                                                                                                                                         gcgcgGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTAACCAATTGcc  <  3:41701/147‑1 (MQ=255)
                                                                                                                                         gcgcgGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTAACCAATTGcc  <  8:299110/147‑1 (MQ=255)
                                                                                                                                         gcgcgGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTAACCAATTGcc  <  7:147844/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
TTAAGAACAAGCTCGTGCGAAACGCCGCTGCGAACTTCACTCTTACCAATTGCCAACGGAAGAATTAAGCGCATCTCTCCCGCAAGGACTTTCTTGTCACGCAGCATATGCGGTAAATACGCCTGCGCGGACATTTCGCGCGGCCCATTGACCGGTAACCCAGCCCGCTTGAGCAGGGTTATAATACGCTGCGTTTCGGCAGAACTAAACTGCCCGAGACGTTCCGACGTCCGCGCCGCCATCACCATACCCGCAGCGACCGCTTCACCATGTAACCAATTGCC  >  NC_000913/3517426‑3517709

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