Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,892,550 A→G 7.7% T113A (ACC→GCC) ‡ 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,5500AG 7.7% 185.4 / 5.6 65T113A (ACC→GCC) ‡queD6‑carboxy‑5,6,7,8‑tetrahydropterin synthase
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base G (3/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

AACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTG  >  NC_000913/2892401‑2892699
                                                                                                                                                     |                                                                                                                                                     
aaccaaccTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCa                                                                                                                                                       >  2:1333623/1‑150 (MQ=255)
  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                                                                                                             <  2:1279896/150‑1 (MQ=255)
                                          ttCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTcg                                                                                                             >  1:622199/1‑150 (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:411578/150‑1 (MQ=255)
                                                                                                   ttAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTggg                                                    <  2:1012196/150‑1 (MQ=255)
                                                                                                     aaaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGc                                                  <  1:805805/150‑1 (MQ=255)
                                                                                                      aaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCa                                                 <  1:340532/150‑1 (MQ=255)
                                                                                                      aaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCa                                                 <  1:725408/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:758719/4‑150 (MQ=255)
                                                                                                                                           tacaaCCTGCGGTGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGt            >  1:763840/4‑150 (MQ=255)
                                                                                                                                            gAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGCCTGCGACATTATTATTATTGGTg           >  1:1070368/1‑150 (MQ=255)
                                                                                                                                                 cTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGt      <  2:451393/150‑1 (MQ=255)
                                                                                                                                                  tGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTa     >  1:369619/1‑150 (MQ=255)
                                                                                                                                                  tGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTa     >  1:842189/1‑150 (MQ=255)
                                                                                                                                                   gCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTAt    >  1:862288/1‑150 (MQ=255)
                                                                                                                                                     aCCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTg  >  1:1340506/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTG  >  NC_000913/2892401‑2892699

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: