Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,390,401 C→A 6.5% intergenic (‑41/+56) 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,4010CA 6.5% 413.0 / 15.3 146intergenic (‑41/+56)psd/rsgAphosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base C (73/62);  new base A (5/5);  total (79/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 = 3.63e-01

Only 100 of 151 total aligned reads displayed.

CCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATT  >  NC_000913/4390252‑4390550
                                                                                                                                                     |                                                                                                                                                     
ccAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGcc                                                                                                                                                       <  2:494192/150‑1 (MQ=255)
ccAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGcc                                                                                                                                                       >  2:361612/1‑150 (MQ=255)
      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                                                                                                                                             >  2:415489/1‑150 (MQ=255)
          cAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                             >  1:145945/1‑150 (MQ=255)
             ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCATGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATc                                                                                                                                          <  2:755816/150‑1 (MQ=255)
             ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATc                                                                                                                                          >  1:1777566/1‑150 (MQ=255)
              cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATcc                                                                                                                                         <  1:855103/150‑1 (MQ=255)
              cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATcc                                                                                                                                         <  2:1614707/150‑1 (MQ=255)
              cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATcc                                                                                                                                         <  1:1130855/150‑1 (MQ=255)
              cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATca                                                                                                                                         <  2:1522976/150‑2 (MQ=255)
              cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGGCGTTCCTGATCCGGCCTACATGCCGGATcc                                                                                                                                         <  2:1251236/150‑1 (MQ=255)
               aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATcag                                                                                                                                        <  1:536851/150‑3 (MQ=255)
               aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATcag                                                                                                                                        <  2:118708/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)
                             cttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGAc                                                                                                                          <  2:1716445/150‑1 (MQ=255)
                                 cgcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGAtt                                                                                                                      <  2:1755739/150‑1 (MQ=255)
                                  gcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTc                                                                                                                     >  2:1295199/1‑150 (MQ=255)
                                  gcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTc                                                                                                                     >  1:390309/1‑150 (MQ=255)
                                   ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGCCGATTCt                                                                                                                    >  2:185458/1‑150 (MQ=255)
                                   ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCt                                                                                                                    >  1:1054371/1‑150 (MQ=255)
                                   ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGATTAATGTACTTTAGCTAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCt                                                                                                                    <  2:1715952/150‑1 (MQ=255)
                                    cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAgggg                                                                                                                            >  1:128472/1‑141 (MQ=255)
                                    cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAgggg                                                                                                                            <  2:128472/141‑1 (MQ=255)
                                    cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTa                                                                                                                   <  1:702190/150‑1 (MQ=255)
                                    cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCAGGATCCTGAAAAAAAGTGGACGATtata                                                                                                                   >  1:795941/1‑150 (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)
                                                         cccATAGTTTCGGCAGAATGTAATGTAGCGAAAGTTAAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGGCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtc                                                                                              <  2:309039/149‑1 (MQ=255)
                                                          ccATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGATTAATTgtca                                                                                             >  1:1509268/1‑150 (MQ=255)
                                                           cATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGAGCCGGAAAAAAAGGGGACGATTCTAACGCCGGGTAGCTTAATTgtcag                                                                                            >  2:1544047/1‑150 (MQ=255)
                                                             tAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtcagtc                                                                                          >  1:731250/1‑150 (MQ=255)
                                                               gTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCAt                                                                                        >  2:1873964/1‑150 (MQ=255)
                                                                tttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATc                                                                                       >  1:1094539/1‑150 (MQ=255)
                                                                 ttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATcc                                                                                      <  2:887660/150‑1 (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                                                                          <  1:792915/150‑1 (MQ=255)
                                                                             tACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCATCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAg                                                                          <  2:942697/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)
                                                                                       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)
                                                                                               aaTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTAcc                                                        >  2:1877290/1‑150 (MQ=255)
                                                                                                aTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGCCGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCt                                                       >  1:765917/1‑150 (MQ=255)
                                                                                                   aaTTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTgcg                                                    >  1:1078323/1‑150 (MQ=255)
                                                                                                      ttAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCa                                                 >  2:968298/1‑150 (MQ=255)
                                                                                                       tAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACATGCGCCAt                                                >  1:1694464/1‑150 (MQ=255)
                                                                                                       tAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGGCAGTCATCCGTAGGAGATAAGTTTTTACGCGTGTTTCCCTGCGCCag                                                >  2:144846/1‑149 (MQ=255)
                                                                                                         aCAACGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGc                                              <  2:243376/150‑1 (MQ=255)
                                                                                                          catggTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCt                                             <  2:948390/147‑1 (MQ=255)
                                                                                                            aGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCttt                                           >  1:1787037/1‑150 (MQ=255)
                                                                                                               tAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCa                                        >  2:563852/1‑150 (MQ=255)
                                                                                                                       aGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACgg                                <  1:851517/150‑1 (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)
                                                                                                                        ggCCATTGTATTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTATTTGTCAGTCATCCGTATCAGATAAGTTTTTACGCGGGTTTACCTGCGCGATGCTTTCCCGATTACggg                               >  2:869498/1‑149 (MQ=255)
                                                                                                                          ccATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGa                             <  2:1612341/150‑1 (MQ=255)
                                                                                                                             ttGTTTTGGCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAg                          <  1:1816792/150‑1 (MQ=255)
                                                                                                                               gTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGtt                        >  1:1352780/1‑150 (MQ=255)
                                                                                                                               gTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGtt                        <  2:159188/150‑1 (MQ=255)
                                                                                                                               gTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACATGCGCCATGCTTTCCAGAATACGGTGATAGtt                        >  2:204504/1‑150 (MQ=255)
                                                                                                                                        ttcaggATCCGGCATGTATGCAGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCAGCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCAATGCTTTCCAGAATACGGTGATAGTTTTAGAAACg               >  2:1146927/6‑150 (MQ=255)
                                                                                                                                           cTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt            >  1:1394794/1‑150 (MQ=255)
                                                                                                                                           cTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt            >  1:1691417/1‑150 (MQ=255)
                                                                                                                                           aggATCCGGCATGTAGGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt            >  1:395373/3‑150 (MQ=255)
                                                                                                                                            tGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGATACGGGtt           >  2:1135036/1‑150 (MQ=255)
                                                                                                                                            tGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGtt           >  1:828659/1‑150 (MQ=255)
                                                                                                                                            tGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGtt           >  1:1198982/1‑150 (MQ=255)
                                                                                                                                            ggATCCGGCATGTAGGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGtt           >  2:568618/2‑150 (MQ=255)
                                                                                                                                              aTCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTc         >  2:1232691/1‑150 (MQ=255)
                                                                                                                                                ccGGCGTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcg       >  2:790289/1‑150 (MQ=255)
                                                                                                                                                  ggCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgcg     >  2:440557/1‑150 (MQ=255)
                                                                                                                                                    ccTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTATTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAt   >  1:348741/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATT  >  NC_000913/4390252‑4390550

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: