Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,813,2720CT30.9% 47.5 / 15.3 42S181L (TCG→TTG) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase and phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base C (4/25);  new base T (13/0);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 9.33e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGTGAATGCTTCT  >  NC_000913/3813039‑3813522
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                              
gattgccctgcgtccgacgaTTTGTCGGTCTTCCACTGCGCCTGTACCCGGCCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTTCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGGCAACCATTGTGGATATGGCGGTAGCGCATGTTTCGCCCGTCAACGACCTGa                                                                                                                                                                                                                                              <  2:73154/237‑1 (MQ=255)
 aaagccctgggaccgactagttGTCTGGCTACCACCGCCCGTGTAGCCGTGCCCCACGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa                                                                                                                                                                                                                                             <  2:156638/231‑1 (MQ=255)
           tATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGGATTTTTTGCCCGTccaacacccgcaaaatctgaac                                                                                                                                                                                                                                    >  1:149796/1‑230 (MQ=255)
               gacgaTTTGTCGGGATACCCCGGCGCCTGTAGCCGTGCCCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTa                                                                                                                                                                                                                              <  1:378550/251‑1 (MQ=255)
                acgaTTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGCGCATTTTTTTCCCGCCAACCCCCTTGAACATCTGAACATTAt                                                                                                                                                                                                                             >  1:392195/1‑251 (MQ=255)
                      gtgtgtgtcaacgccggggccggTAGCGGTGCTCCCCGCCATGAGCCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGTGCCAGACAGTGGCAGTCAGGCTTGTGTTGATATCGGTCCTGGGCTAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTATCGCATTTTTCGCCCGTCAAGGACCTGAAACATCTGAGCAGTATGAtt                                                                                                                                                                                                                         <  2:199695/228‑1 (MQ=255)
                                ccccccgccccgtatcctgTGCCCCCCGCTGTGAACCAGCATAGGTCCCGTGCCGTTGCCCCGAATCATAATTAGGAGGTGCTTGTTCCCCGTGGTTTCCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATACCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTTGATATGGCGGTAGCGCATTTTCCCCCCGTCACCGACCTGAACCATCTGACCATTATGATTACCGCCGGccc                                                                                                                                                                                                              <  2:392195/233‑1 (MQ=255)
                                  acCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGACCATTTTTTTCCCCTCAACCACCTGAAAAATCTGTACATTTTGATTACCCcccgcccgc                                                                                                                                                                                                            >  1:317415/1‑250 (MQ=255)
                                       cgcCTGTAGCGGGGCCCCCCGCCATGAACCAGCAGATGTCCCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGt                                                                                                                                                                                                       <  1:106061/251‑1 (MQ=255)
                                         ccTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCGCAACGACCTGAAACAACTGAACATTATGATTACCGCCGGCCCGAcgcgggg                                                                                                                                                                                                     >  1:88491/1‑248 (MQ=255)
                                                       cgccGCCATGAACCAGCAGAGGTACCGTGCCGCTGCCACCCAGCATAATTTAGAGGTGCTTGCTCCCCGGGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTACCCATTTTGGGTATGGCGGTAGCGCATTTTTCGCCCGTCAACGCCCTGAAACATCTGAACATGATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg                                                                                                                                                                                       <  1:284579/249‑1 (MQ=255)
                                                       ccccGCCATGAACCAGCAGATGTCGCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCGGGGGGCCAGACAGTGGCAGTCAGGCTTGGGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGCCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg                                                                                                                                                                                       <  2:415367/251‑1 (MQ=255)
                                                           gCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGCGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTACTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCCACACCCTGTAACATCTGTAAACTATTATTAACGCCCGCCcgccgccggaacagcccgctccggggcg                                                                                                                                                                                   >  1:464876/1‑229 (MQ=255)
                                                              aTGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGGGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTta                                                                                                                                                                                <  1:4912/251‑1 (MQ=255)
                                                                              aCCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGAGGACTTTTTTGCCGGCAACGCACGGAAACAATTTGACATTATGATTACCGCCGGCCCCCCCCGTGAACCGCTCGAACCCGCGCGTTATATCTCTaaaccc                                                                                                                                                                   >  1:17860/1‑244 (MQ=255)
                                                                                cgggCCGCTGCCACGCTGCCTTTTTTCGAGGTGCTTGCTTCCCGTTGTTGGCTCTGCTGGGGGCTAGACAGAGTCCGTTGTGCTTTGGGGGATATCGGTCCTGGGCGATGGCTGGATCCGTTAACCAGTGTGTTAATGGCGTTAGCGCATTTTTCGCCCGTCACCGTCCTGAAACATCTTAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCTGTGCGTTATATCGCTAATCACAGcgc                                                                                                                                                              <  2:336642/248‑3 (MQ=255)
                                                                                         gCCACGCAGCATAATTTAGAGGTTCTTGCTCCCCGTGGTTTGCTCTTCTGGGGGCCAGACAGTGGCAGTTAGGCTTGTTGTGATATCGGTCCGGGGCGAATGCTGGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCaaa                                                                                                                                                       <  2:70119/249‑2 (MQ=255)
                                                                                            ccGCAGCATAATTGAGGGGTGCTTGCTTCCCGTGGTTGGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCGGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGCCCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATggg                                                                                                                                                 <  2:358352/250‑1 (MQ=255)
                                                                                                 gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTGGCGCCTTTTTTGCCCGTCCACCACCTTAAACATCCGAACATTATGATTACCCCCCGCCCCACCCGCGAACCGCTCCATCCGGGGCGCTATATTTCTAATTACAgacccgggcaaaaggggtttt                                                                                                                                             >  1:80018/1‑231 (MQ=255)
                                                                                                 gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTg                                                                                                                                            <  1:131936/251‑1 (MQ=255)
                                                                                                          atgggTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACGGGGGCAGTCAGGCTTGTGGTGATAGCGGGCCGGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCGGAACATGAGGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccg                                                                                                                                   <  1:192771/248‑1 (MQ=255)
                                                                                                          agagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCCTTTTTTGCCCCTCAAACACCCGAAACATCTAAAAATTATGCTTACCCCCGGCCCCACGCGGGAAACGCCCCACCCGGTGCGGTATATCTCTAAACACACCTccccgaaaaaggggtttgttatcgccc                                                                                                                                    >  1:203782/1‑226 (MQ=255)
                                                                                                           gagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGGAAACACTGAACATTATTATTACCGCCGGCCCGACGCGGGGACCGCTCGATCCGGTGGGTTTTTTCTCTAATTACACCTCCCGGAAGAATGGGTTTGCTATtgcccc                                                                                                                                  >  1:246466/1‑249 (MQ=255)
                                                                                                            agGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCCACGACCTGAAAAATCTGAACACTATTTTTACCGCCGGGCCGCCGCCGGAAACGCTCGCTCCGCCGGTTTGTATCTTTAATCACAGCCCCGGCCAGAGGGGTTTTgagttccgcgcc                                                                                                                                 >  2:377737/1‑240 (MQ=255)
                                                                                                                        cccGTGGTTGGCGCATCGGGGGGCCAGAAAGTGGCAGTCAGGCTTGTGGGGATGTGGGTCTTGGGCGAATCCTCTGTGCGTGACCCATTGGGGTTATGGCGGTAGCGCATTTTCCGCCCGTAAACGACCGGAAACTTCGGAACATAAGGATTCCCGCCGGCCCGAGGCGTTAACCGCTTGATGCGGTGCGTTATATCTCTATTCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCcgt                                                                                                                     <  1:451206/251‑1 (MQ=255)
                                                                                                                                               ccAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGGCCTGGGCGAATGCTCGATCCGTTAACCATTGGGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCGGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGa                                                                                                                                                      <  2:437020/196‑1 (MQ=255)
                                                                                                                                                     aGTGGCAGTCAGGCTTGTGGTGATATCGGTCCCGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGATCGGCGTTTTTCCCCGTCCACGACCCGCACCATCTTAAAAATATTGTTTCCCCCGGCCCCCCGCGCGAACCGCCCGATCCGGGGCGGTCTATCTCTAATACCAGCCGCGGCGCGACGGGGTTTGATTTCGCCGGCCCCGCCGCCCGTCCTGGcgcgggcgccccgctggtgtccgg                                                                                        >  2:101780/1‑231 (MQ=255)
                                                                                                                                                       gggcaggcaGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATTGGGCGGTAGCGCATTTTCCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTc                                                                                      <  1:201121/250‑1 (MQ=255)
                                                                                                                                                           aggcAGGCTTGTGGGGAGAGGGGTCCGGGGGGAATGCTCGATCCGTTAACCATTGTGGATAGGGCGGGAGCGCTTTTTTCGCCCGTAAACGACCTGAAACATCTGAACATTATTATACCCGCCGGCCCGACGCGGGAACCGCTCGATCCGGGGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGGGGCGCGAACGTCACGCGGGTATCAGGTCCGGt                                                                                  <  2:393514/248‑1 (MQ=255)
                                                                                                                                                              cAGGCTTGTGGTGATTTCGGTCCTGGGCGAATGCTCGATCCTTTAACCATGGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTc                                                                               <  1:85397/251‑1 (MQ=255)
                                                                                                                                                                     ggggtgATATCGGGCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTGCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTAcc                                                                         <  1:452556/248‑1 (MQ=255)
                                                                                                                                                                        gtgATATGGGGCCTGGGCGAATGCTCGACCCTTTAACCATTGTGGATATGGCGGTAGCGCATTTTCCGCCCGTCAACGCCCTGAAACATCTAAACATTATGATTCCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACg                                                                     <  2:26955/251‑1 (MQ=255)
                                                                                                                                                                               gggTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATAGGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGt                                                             <  1:304279/250‑1 (MQ=255)
                                                                                                                                                                                    cTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACCACCTGAAACATCTGGACATTATGGTTACCGCCGGCCCCACGCGTGAACCGCCCGATCCGGGGCGTTATATTTCTAATCACCACCCCGGCCAGATGGGGTTTTTTATCGCCGCCCCCCCCGCCCCTCGTGTGGCGTACGTCCCGCTGGTGTCAGGGCCGGGTTCCCTACCCCCGCCCCCGgttgtt                                                        >  2:76826/1‑251 (MQ=255)
                                                                                                                                                                                        gCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTCGCGCATTTTTTGCCCGACAACGCCCTTGCACATATGACAATTATGATTATCGCCGGCCCCACCCGTGAACCCCCCCATCCCCGGCGGTTCATATCCAATCACCCCACCGGCCAGATGGGATTTGCTCTCGCCCCCGCCCCCGGCCCCCGCGGCGCGGACGTCCCGCTGGGGTCAGGGCCGGTTTCACTCCCGCCGCCCCCCTTTGtttacc                                                    >  2:24180/1‑247 (MQ=255)
                                                                                                                                                                                           aaggCTCGTTCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTAAACCATCTGACCATTATGATTCCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACgtg                                                 <  1:101780/248‑1 (MQ=255)
                                                                                                                                                                                              gctggttCCGTTAACCATTGGGGTTTTGGGGGTAGCGCTTTTTCCGCCCGTCAACGACCTGAAACATCTGAACATTATGATAACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTGATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCTACCGTGTGTTAAACGTGtt                                               <  2:84978/244‑1 (MQ=255)
                                                                                                                                                                                                  gATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCCGCCCCACGCGTGAACCCCCCGATCCGGGGCGTTTTATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCTCCCCCGCCCCTGCCCCGCGTGGCGCGGACGTCACGCTGGTATCCGGTCCGGGTTCACTACCGACGCCCCCGCTTTTTTAACGTGTtgaggt                                          >  1:57814/1‑248 (MQ=255)
                                                                                                                                                                                                  gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACAACCTGAAACATCTGAACATTATGATTACCCCCGGCCCGACGCGTGAACCGCCCGATCCCGTGCGGTATATCTCTAATCACAGCCCCGGCAAGATGGGGTTTTCTATCTCCGCCGCCGCTGCCCGTCGGGGCGCGGACGTCCCGCTGGTATCAGGTCCGGTTTCACTTCCGACGCCACCCGTTTTTTAACGTGTtgatgt                                          >  1:197919/1‑251 (MQ=255)
                                                                                                                                                                                                                     atatGGGGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAGCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATTGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCGGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATg                       <  1:230886/251‑1 (MQ=255)
                                                                                                                                                                                                                       atggcggccgggcATTTTTCTCCCGTCAAAGACCTTAAAAATCTGCACATTATGATTACCGCCGGCCCGCCGCGAGACCAGCCCGATCCGGTGCGTTATATCTCTAATCACAGCTCCCGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGACACGCCGGTATCAGGTCCGGTTTCACTACCGACTCCCCCCATTGTTAAACGTGTTGATGTGCTTACCCCGCTGGAAATgg                      >  2:190544/12‑250 (MQ=255)
                                                                                                                                                                                                                                  gcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGTg           >  2:140820/1‑250 (MQ=255)
                                                                                                                                                                                                                                           ccGCCCGTCAACGCCCTGAAACATCTGAACATTATGATTCCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGTGAATGCtt    <  1:76826/249‑1 (MQ=255)
                                                                                                                                                                                                                                            cGCCCGTCAACGACCTGAACCTTCTGAACATTATGATTCCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAGCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGTGAATGCTTCt  <  1:84244/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                              
GAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGTGAATGCTTCT  >  NC_000913/3813039‑3813522

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: