Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,390,404 C→T 6.8% intergenic (‑44/+53) psd ← / ← rsgA phosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,390,4040CT 6.8% 465.5 / 14.5 148intergenic (‑44/+53)psd/rsgAphosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base C (75/62);  new base T (5/5);  total (81/67)
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.00e-01

Only 100 of 150 total aligned reads displayed.

TTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTC  >  NC_000913/4390258‑4390553
                                                                                                                                                  |                                                                                                                                                     
ttGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATg                                                                                                                                                    >  2:1510148/1‑150 (MQ=255)
 tGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGc                                                                                                                                                   >  1:1493282/1‑150 (MQ=255)
  gTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGcc                                                                                                                                                  <  1:504583/150‑1 (MQ=255)
   tCAGCCATCCTGCCCGCTTGCTTGCGCCCCCACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCg                                                                                                                                                 <  1:1147793/150‑1 (MQ=255)
    cAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                                <  2:655009/150‑1 (MQ=255)
    cAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                                >  1:145945/1‑150 (MQ=255)
     aGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGa                                                                                                                                               >  2:1856596/1‑150 (MQ=255)
       ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCATGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATc                                                                                                                                             <  2:755816/150‑1 (MQ=255)
       ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATc                                                                                                                                             >  1:1777566/1‑150 (MQ=255)
       ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATc                                                                                                                                             <  2:1895303/150‑1 (MQ=255)
       ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATc                                                                                                                                             <  2:1723867/150‑1 (MQ=255)
        cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATcc                                                                                                                                            <  1:1130855/150‑1 (MQ=255)
         aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCt                                                                                                                                           <  1:732991/150‑1 (MQ=255)
         aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATcag                                                                                                                                           <  2:118708/150‑3 (MQ=255)
         aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATcag                                                                                                                                           <  1:536851/150‑3 (MQ=255)
           ccTGCCCGCATGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGGTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATcagga                                                                                                                                         <  1:103515/150‑5 (MQ=255)
            cTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaa                                                                                                                                        >  2:454239/1‑150 (MQ=255)
              gCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTAAATGCCGGATCCTGaaaa                                                                                                                                      >  2:517631/1‑150 (MQ=255)
               cccGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaaaaa                                                                                                                                     >  1:802438/1‑150 (MQ=255)
                  gcttgcttGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAg                                                                                                                                  <  2:1194999/150‑1 (MQ=255)
                     tgcttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAgggg                                                                                                                               >  2:1323512/1‑150 (MQ=255)
                      gcttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGa                                                                                                                              >  1:747676/1‑150 (MQ=255)
                            gcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTc                                                                                                                        >  1:390309/1‑150 (MQ=255)
                             ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCt                                                                                                                       >  1:1054371/1‑150 (MQ=255)
                             ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGATTAATGTACTTTAGCTAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCt                                                                                                                       <  2:1715952/150‑1 (MQ=255)
                              cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAgggg                                                                                                                               <  2:128472/141‑1 (MQ=255)
                              cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAgggg                                                                                                                               >  1:128472/1‑141 (MQ=255)
                              cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTa                                                                                                                      <  1:702190/150‑1 (MQ=255)
                                  aCCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGAGTCTAacga                                                                                                                  >  1:1233518/1‑150 (MQ=255)
                                   cccGCTAGGCGAGTAATCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAacgac                                                                                                                 <  2:838997/150‑1 (MQ=255)
                                    ccGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAacgacg                                                                                                                <  1:1510148/150‑1 (MQ=255)
                                     cGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACgg                                                                                                               <  2:822101/150‑1 (MQ=255)
                                     cGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTGAAAATGAATTTAACAAGGGAGCCTCCAGGCCATGGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACgg                                                                                                               <  1:550667/150‑1 (MQ=255)
                                        cAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTa                                                                                                            >  1:159188/1‑150 (MQ=255)
                                         aGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAg                                                                                                           <  1:263212/150‑1 (MQ=255)
                                          ggCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGc                                                                                                          >  1:302334/1‑150 (MQ=255)
                                            cGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCtt                                                                                                        >  2:600248/1‑150 (MQ=255)
                                            cGAGTAAGACATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCtt                                                                                                        <  1:1127277/150‑1 (MQ=255)
                                             gAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTa                                                                                                       <  1:128413/150‑1 (MQ=255)
                                               gTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAAt                                                                                                     <  2:942407/150‑1 (MQ=255)
                                                    ccATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGATTAATTgtca                                                                                                >  1:1509268/1‑150 (MQ=255)
                                                     cATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtcag                                                                                               <  1:654009/150‑1 (MQ=255)
                                                     cATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGAGCCGGAAAAAAAGGGGACGATTCTAACGCCGGGTAGCTTAATTgtcag                                                                                               >  2:1544047/1‑150 (MQ=255)
                                                        aGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtccgtca                                                                                            <  2:1778007/150‑1 (MQ=255)
                                                         gTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCAt                                                                                           >  2:1873964/1‑150 (MQ=255)
                                                          tttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATc                                                                                          >  1:1094539/1‑150 (MQ=255)
                                                           ttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATcc                                                                                         <  2:887660/150‑1 (MQ=255)
                                                             ctgCAGAATGTACTGTATCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCTGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGt                                                                                       <  2:316578/148‑1 (MQ=255)
                                                              ggCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTa                                                                                      >  2:127823/1‑150 (MQ=255)
                                                              ggCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTa                                                                                      >  2:851517/1‑150 (MQ=255)
                                                               gCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTAt                                                                                     >  2:818959/1‑150 (MQ=255)
                                                                     tGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaaa                                                                               <  2:1807033/150‑1 (MQ=255)
                                                                     tGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaaa                                                                               >  2:1751927/1‑150 (MQ=255)
                                                                       tACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAg                                                                             >  2:433342/1‑150 (MQ=255)
                                                                       tACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAg                                                                             <  1:792915/150‑1 (MQ=255)
                                                                           gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtttt                                                                         >  2:289443/1‑150 (MQ=255)
                                                                             aGCGAAAGTTTAAATGAATTTAACAAGGTCGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACAGGCAGGATCCTGCAAAATAGGGGCCGATTCTAACGCCGTTTAGGTGAATTGTCAGTCAGCAGTATCAGAAAAGGTTTTa                                                                       >  1:116984/1‑150 (MQ=255)
                                                                             aGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTa                                                                       >  1:944064/1‑150 (MQ=255)
                                                                              gCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTAc                                                                      >  2:1597139/1‑150 (MQ=255)
                                                                                 aaaGTGTAATTGAATTTAAGAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGGCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTAcgcg                                                                   <  1:22314/150‑1 (MQ=255)
                                                                                  aaGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGt                                                                  >  2:645310/1‑150 (MQ=255)
                                                                                       tAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGTCGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTa                                                             <  1:1721170/150‑1 (MQ=255)
                                                                                         aaTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTAcc                                                           <  1:361612/150‑1 (MQ=255)
                                                                                             aaTTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTgcg                                                       >  1:1078323/1‑150 (MQ=255)
                                                                                                ttAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCa                                                    >  2:968298/1‑150 (MQ=255)
                                                                                                 tAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACATGCGCCAt                                                   >  1:1694464/1‑150 (MQ=255)
                                                                                                    catggTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCt                                                <  2:948390/147‑1 (MQ=255)
                                                                                                      aGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCttt                                              >  1:1787037/1‑150 (MQ=255)
                                                                                                        gTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTcc                                            <  2:145945/150‑1 (MQ=255)
                                                                                                         tAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCa                                           >  2:563852/1‑150 (MQ=255)
                                                                                                                 aGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACgg                                   <  1:851517/150‑1 (MQ=255)
                                                                                                                  ggCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTACCGACGGTTAGCTTTATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGt                                  >  2:1314389/1‑150 (MQ=255)
                                                                                                                  ggCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGt                                  <  2:630709/150‑1 (MQ=255)
                                                                                                                  ggCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGt                                  <  1:1597139/150‑1 (MQ=255)
                                                                                                                  ggCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTGAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATTCTTTCCAGAATACGGt                                  >  1:647807/1‑150 (MQ=255)
                                                                                                                  ggCCATTGTTTTGTCGTGCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGt                                  <  1:1757162/150‑1 (MQ=255)
                                                                                                                    ccATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGa                                <  2:1612341/150‑1 (MQ=255)
                                                                                                                    ccATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGa                                <  2:1533315/150‑1 (MQ=255)
                                                                                                                       ttGTTTTGGCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAg                             <  1:1816792/150‑1 (MQ=255)
                                                                                                                         gTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGtt                           >  1:1352780/1‑150 (MQ=255)
                                                                                                                         gTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACATGCGCCATGCTTTCCAGAATACGGTGATAGtt                           >  2:204504/1‑150 (MQ=255)
                                                                                                                             tGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCg                       <  2:1137723/150‑1 (MQ=255)
                                                                                                                                  ttcaggATCCGGCATGTATGCAGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCAGCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCAATGCTTTCCAGAATACGGTGATAGTTTTAGAAACg                  >  2:1146927/6‑150 (MQ=255)
                                                                                                                                    ccTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACggg                >  2:1658865/1‑150 (MQ=255)
                                                                                                                                     cTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt               >  1:1691417/1‑150 (MQ=255)
                                                                                                                                     cTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt               >  1:1394794/1‑150 (MQ=255)
                                                                                                                                     aggATCCGGCATGTAGGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt               >  1:395373/3‑150 (MQ=255)
                                                                                                                                      tGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGtt              >  1:1198982/1‑150 (MQ=255)
                                                                                                                                      ggATCCGGCATGTAGGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGtt              >  2:568618/2‑150 (MQ=255)
                                                                                                                                      ggATCCGGCATGTAGGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGtt              >  2:260684/2‑150 (MQ=255)
                                                                                                                                       gATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGttt             >  1:796210/1‑150 (MQ=255)
                                                                                                                                        aTCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTc            >  2:1232691/1‑150 (MQ=255)
                                                                                                                                         tCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTcc           <  1:769374/150‑1 (MQ=255)
                                                                                                                                          ccGGCGTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcg          >  2:790289/1‑150 (MQ=255)
                                                                                                                                           cGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgc         <  1:839709/150‑1 (MQ=255)
                                                                                                                                           cGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgc         >  1:867755/1‑150 (MQ=255)
                                                                                                                                            ggCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgcg        >  2:440557/1‑150 (MQ=255)
                                                                                                                                              ccTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTATTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAt      >  1:348741/1‑150 (MQ=255)
                                                                                                                                                  cATGCCGGATCCTGAAAAAAAGTGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTc  <  2:390309/150‑1 (MQ=255)
                                                                                                                                                  cATGCCGGATCCTGAAAAAAAGGGGACGATACTAACGACGGTTAGCTTAATTGTCAGACATCAGTAACAGAAAAGATTTTACGCGTTTTAAACTGCGCCATGCTATCCAGAATACGGTGATAGTATACGAAACGGCTTACCGCGATTATc  >  2:1243636/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTC  >  NC_000913/4390258‑4390553

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