Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,892,552 C→T 7.7% T113T (ACC→ACT) ‡ queD → 6‑carboxy‑5,6,7,8‑tetrahydropterin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,892,5520CT 7.7% 192.1 / 4.9 65T113T (ACC→ACT) ‡queD6‑carboxy‑5,6,7,8‑tetrahydropterin synthase
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base T (3/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

CCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCA  >  NC_000913/2892403‑2892701
                                                                                                                                                     |                                                                                                                                                     
ccaaccTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCAcc                                                                                                                                                       >  1:96357/1‑150 (MQ=255)
           tcgcCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCGGTGCAGGtttctt                                                                                                                                            <  2:121611/149‑5 (MQ=255)
            cgcCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGGTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGtatc                                                                                                                                           <  1:408916/150‑1 (MQ=255)
             gcCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCGGTGCAGGtttctttt                                                                                                                                          <  1:480117/150‑7 (MQ=255)
              ccTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTTAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGtatcta                                                                                                                                         >  2:1258062/1‑150 (MQ=255)
                     cACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCgcggc                                                                                                                                  <  1:1015936/150‑1 (MQ=255)
                      aCCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCgcggcg                                                                                                                                 >  1:1374772/1‑150 (MQ=255)
                       ccATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGa                                                                                                                                >  2:998429/1‑150 (MQ=255)
                       ccATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGa                                                                                                                                >  2:144373/1‑150 (MQ=255)
                         aTTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAAt                                                                                                                              >  2:1355779/1‑150 (MQ=255)
                              ctcAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATaa                                                                                                                         <  2:677844/150‑1 (MQ=255)
                                 aaTGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAagag                                                                                                                      <  1:1333623/150‑1 (MQ=255)
                                  aTGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAgt                                                                                                                     >  1:1334161/1‑150 (MQ=255)
                                   tGATATTCCAGGCCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAgtg                                                                                                                    <  2:261573/150‑1 (MQ=255)
                                    ggttttCCAGACCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACGCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGGAGGTTGTATCTATCGCGGCGATTGATAAGAgtgt                                                                                                                   >  1:162903/5‑150 (MQ=255)
                                     atatTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAgtgtg                                                                                                                  <  1:998429/150‑1 (MQ=255)
                                        ttCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTcg                                                                                                               >  1:622199/1‑150 (MQ=255)
                                        ttCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTcg                                                                                                               <  2:1279896/150‑1 (MQ=255)
                                         tCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTCTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTcgg                                                                                                              >  1:527020/1‑150 (MQ=255)
                                            aGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTcggcgg                                                                                                           <  2:172902/150‑1 (MQ=255)
                                              gTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTc                                                                                                         <  1:1361631/150‑1 (MQ=255)
                                               tCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCa                                                                                                        <  1:617291/150‑1 (MQ=255)
                                                  ggAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAAtt                                                                                                     >  2:340532/1‑150 (MQ=255)
                                                    aaaaCCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAACCCCGTTGTGCCGCTGTTAAGTGCGGGGATGGTAAAAGAAACCTGCACCGCAGGGTGTATCTATCGCGGCGAATGATAAGAGTGCGCCGGCGGTCAATTTc                                                                                                   <  2:892878/150‑1 (MQ=255)
                                                    aaaaCCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTc                                                                                                   <  1:636395/150‑1 (MQ=255)
                                                      aaCCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTccc                                                                                                 <  1:546183/150‑1 (MQ=255)
                                                             ccAGCGAGGGGTTAGCAAAATGGATTTGGGAGCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAACCCTGCACCGCAGGTTGTATCTCTCGCGGCGAATGATAAGAGTGAGGCGGCGGTCAATTTCCCTTAAGTa                                                                                          <  2:156519/150‑1 (MQ=255)
                                                                  gAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGta                                                                                                                                             >  2:75653/1‑94 (MQ=255)
                                                                  gAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGta                                                                                                                                             <  1:75653/94‑1 (MQ=255)
                                                                          aGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGt                                                                                                          <  1:1211031/121‑1 (MQ=255)
                                                                          aGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGt                                                                                                          >  2:1211031/1‑121 (MQ=255)
                                                                           gCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACGGGAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGGGGTGAATTTCCCTTAAGTAACGCTATGTTAggg                                                                            <  2:529557/150‑1 (MQ=255)
                                                                           gCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAATAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAggg                                                                            <  2:941858/150‑1 (MQ=255)
                                                                           gCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAggg                                                                            >  2:150734/1‑150 (MQ=255)
                                                                                  ggATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGgtgttgtg                                                                     >  2:86907/1‑150 (MQ=255)
                                                                                  ggATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGgtgttgtg                                                                     <  2:796829/150‑1 (MQ=255)
                                                                                       tGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTg                                                                >  2:494351/1‑150 (MQ=255)
                                                                                        gggATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTACCGCTATGTTAGGGTGTTGTGTTCTgg                                                               >  1:184691/1‑150 (MQ=255)
                                                                                         ggATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGa                                                              <  2:96357/150‑1 (MQ=255)
                                                                                           aTCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGata                                                            <  2:1374772/150‑1 (MQ=255)
                                                                                                gTTAAACCCTTTGTGCCGCGGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTgg                                                       <  2:484458/150‑1 (MQ=255)
                                                                                                gTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTgg                                                       >  1:451393/1‑150 (MQ=255)
                                                                                                 ttAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTggg                                                      <  2:1012196/150‑1 (MQ=255)
                                                                                                 ttAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTggg                                                      <  2:411578/150‑1 (MQ=255)
                                                                                                   aaaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGc                                                    <  1:805805/150‑1 (MQ=255)
                                                                                                    aaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCa                                                   <  1:725408/150‑1 (MQ=255)
                                                                                                    aaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCa                                                   <  1:340532/150‑1 (MQ=255)
                                                                                                     aCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCAt                                                  >  2:855797/1‑150 (MQ=255)
                                                                                                              gCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGaa                                         >  2:666302/1‑150 (MQ=255)
                                                                                                                  cTGTTATGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGGGGCGAATGATAAGAGTGTGGCGGGGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAgacg                                     <  2:764394/150‑1 (MQ=255)
                                                                                                                  cTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAgacg                                     <  2:622199/150‑1 (MQ=255)
                                                                                                                           gCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACa                            >  2:411242/1‑150 (MQ=255)
                                                                                                                              gTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACatta                         <  1:1258062/150‑1 (MQ=255)
                                                                                                                                 aTGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACattatta                      <  1:855797/150‑1 (MQ=255)
                                                                                                                                   ggTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACattattatt                    >  1:898652/1‑150 (MQ=255)
                                                                                                                                     tAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGGGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACattattattat                  <  1:1149940/150‑1 (MQ=255)
                                                                                                                                      aaaaGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACattattattatt                 >  2:545314/1‑150 (MQ=255)
                                                                                                                                        atacaaCCTGCGGTGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGCCTGCGACATTATTATTATTgg               >  2:725831/5‑150 (MQ=255)
                                                                                                                                         tacaaCCTGCGGTGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGt              >  1:763840/4‑150 (MQ=255)
                                                                                                                                         tacaaCCTGCGGTGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGt              >  1:758719/4‑150 (MQ=255)
                                                                                                                                          gAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGCCTGCGACATTATTATTATTGGTg             >  1:1070368/1‑150 (MQ=255)
                                                                                                                                               cTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGt        <  2:451393/150‑1 (MQ=255)
                                                                                                                                                tGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTa       >  1:842189/1‑150 (MQ=255)
                                                                                                                                                tGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTa       >  1:369619/1‑150 (MQ=255)
                                                                                                                                                 gCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTAt      >  1:862288/1‑150 (MQ=255)
                                                                                                                                                   aCCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTg    >  1:1340506/1‑150 (MQ=255)
                                                                                                                                                    ccGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGc   <  1:1355779/150‑1 (MQ=255)
                                                                                                                                                    ccGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGc   >  2:1149080/1‑150 (MQ=255)
                                                                                                                                                     cGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCa  <  1:739696/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCA  >  NC_000913/2892403‑2892701

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