Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 788,385 G→T 6.2% T151T (ACC→ACA galM ← aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913788,3850GT 6.2% 267.4 / 5.0 97T151T (ACC→ACAgalMaldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base T (3/4);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 103 total aligned reads displayed.

GGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCT  >  NC_000913/788236‑788533
                                                                                                                                                     |                                                                                                                                                    
ggCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTAtcgtcg                                                                                                                                                      >  2:1179712/1‑150 (MQ=255)
 gCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCgg                                                                                                                                                     <  2:700285/150‑1 (MQ=255)
   aGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTc                                                                                                                                                   <  1:1589427/150‑1 (MQ=255)
    gATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCa                                                                                                                                                  >  1:1248068/1‑150 (MQ=255)
     atatTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                 >  1:1571843/1‑150 (MQ=255)
      tatTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCTACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                >  1:1543461/1‑150 (MQ=255)
      tatTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATCGAAATAGACGTGATTAGTCATATTCACCGAGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                >  1:1177225/1‑150 (MQ=255)
        ttCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACg                                                                                                                                              >  1:686858/1‑150 (MQ=255)
         tCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGAGTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGa                                                                                                                                             <  2:1608729/150‑1 (MQ=255)
               tccAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGGTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGc                                                                                                                                       <  1:700212/149‑1 (MQ=255)
               gCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTCTCGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTga                                                                                                                                       >  2:705158/1‑149 (MQ=255)
               gCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGGGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGataaccg                                                                                                                                       <  1:897729/150‑5 (MQ=255)
               gCCAGAATCTGAAACTGGTCATAGCTCACGTAAGACTGCTCTCCGTCAATATTGAAATAGACGTGATTACGCATATTCACCGGGCAAGGTTTATCAACTGTGGCTCGATAATTAATGGAGATACGTTAATCGTCTGACCGACtataaccg                                                                                                                                       <  2:179674/150‑5 (MQ=255)
                        tGCACCTTGTGATTGCGCACGTCCGACTGCTCGCGGGCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGGTTAGCAACTGTGGCGCGATAAGTCATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcgcg                                                                                                                              <  2:311236/150‑1 (MQ=255)
                        tGCAACTTTTGATTGCGAACGTCAGACTGCTCGCAGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGGTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcgcg                                                                                                                              <  1:1308840/150‑1 (MQ=255)
                         gCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCgcgc                                                                                                                             <  1:35308/150‑1 (MQ=255)
                         gCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTACCGTCGGTCAGACGATATTGCACCGTCgcgc                                                                                                                             <  1:687983/150‑1 (MQ=255)
                          cAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATATGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGcc                                                                                                                            >  2:844072/1‑150 (MQ=255)
                          cAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGcc                                                                                                                            >  1:1275664/1‑150 (MQ=255)
                            aCTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCga                                                                                                                          >  2:491544/1‑150 (MQ=255)
                                gtgATTGCGCACGTCAGCCTGCGCGCCGTCAAGATTGAAATAGACGTGATTAGTCAAATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGTCAGACGATATTGCACCGTCGCGCCGAGAtt                                                                                                                      <  2:264597/150‑1 (MQ=255)
                                   aTTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTAcc                                                                                                                   >  1:1474224/1‑150 (MQ=255)
                                       cgcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCCCCGTCGCGCCGAGATTACCCggg                                                                                                               >  1:322617/1‑150 (MQ=255)
                                         cACGTCAGACTGCTCGCCGTCTAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGaa                                                                                                             >  1:342522/1‑150 (MQ=255)
                                          aCGTCAGACTGCTCGCAGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGTGAAg                                                                                                            >  1:588125/1‑150 (MQ=255)
                                            gTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGcc                                                                                                          >  1:1317007/1‑150 (MQ=255)
                                              cAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCt                                                                                                        <  1:1511065/150‑1 (MQ=255)
                                              cAGACTGCTCGCCGTCAAGATTGAAAAAGACGTGATTACTCATATTCACCGGGCAAGGTATATCAACGGTGGCGCGATAAGTAAGGGACATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCt                                                                                                        <  2:1330152/150‑1 (MQ=255)
                                               aGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                       <  2:32480/150‑1 (MQ=255)
                                                gACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGa                                                                                                      >  1:1686432/1‑150 (MQ=255)
                                                  cTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATc                                                                                                    >  2:246583/1‑150 (MQ=255)
                                                   tGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                   >  2:117587/1‑150 (MQ=255)
                                                    gCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                  >  1:593899/1‑150 (MQ=255)
                                                    gCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                  >  2:1196437/1‑150 (MQ=255)
                                                    gCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                  >  2:480882/1‑150 (MQ=255)
                                                      tCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACca                                                                                                >  1:375474/1‑150 (MQ=255)
                                                      tCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACca                                                                                                <  2:1544190/150‑1 (MQ=255)
                                                      tCGCCGTCAAGATTGAAATAGACGTGATTAGACATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATTGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCAGATCACca                                                                                                >  1:609107/1‑150 (MQ=255)
                                                        gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTCGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCAGGGAAGCCCTGATCACcatc                                                                                              >  1:909959/1‑150 (MQ=255)
                                                         ccGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatca                                                                                             <  1:819428/150‑1 (MQ=255)
                                                         ccGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatca                                                                                             <  1:805337/150‑1 (MQ=255)
                                                           gTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCatcatc                                                                                           <  2:1190188/150‑1 (MQ=255)
                                                           gTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCatcatc                                                                                           >  1:600725/1‑150 (MQ=255)
                                                            tCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCt                                                                                          <  1:1481376/150‑1 (MQ=255)
                                                             cAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTg                                                                                         <  1:246583/150‑1 (MQ=255)
                                                              aaGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGa                                                                                        <  2:148866/150‑1 (MQ=255)
                                                               aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTAATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGaa                                                                                       <  2:1473512/150‑1 (MQ=255)
                                                               aGATTGAAATAGACGTGATTAGTCATATTAACCGCGCAAGGTTTATCAACTGTGTCGCGATAAGTAATGGAGATACGTTTATCGTCGGTCAGACGATATTGCACCTTCGCGCCGAGATTACACGGGAAGCCCTGATCACCATCATCTGaa                                                                                       >  1:1468623/1‑150 (MQ=255)
                                                                gATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAAc                                                                                      >  1:1565164/1‑150 (MQ=255)
                                                                  ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGCTTACCCGGGAAGCCCTGATCACCATCATCTGCACTc                                                                                    >  2:345858/1‑150 (MQ=255)
                                                                    gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGATGCCCTGATCCCCATCATCTGAACTCAg                                                                                  >  2:314121/1‑150 (MQ=255)
                                                                    gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCATGATCACCATCATCTGAACTCAg                                                                                  <  2:1474224/150‑1 (MQ=255)
                                                                     aaaTAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCTCCGAGATTACCAGGGAAGCCCTGATCACCATCAACTGAACTCAgg                                                                                 >  1:173878/1‑150 (MQ=255)
                                                                      aaTAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAggg                                                                                <  2:593899/150‑1 (MQ=255)
                                                                        tAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCa                                                                              >  2:452338/1‑150 (MQ=255)
                                                                         aGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCaa                                                                             >  1:1710636/1‑150 (MQ=255)
                                                                           aCGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                           <  1:291243/150‑1 (MQ=255)
                                                                            cGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACa                                                                          <  2:1705940/150‑1 (MQ=255)
                                                                              tGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGc                                                                        >  2:426320/1‑150 (MQ=255)
                                                                               taTTAGTCATATTCAACGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCa                                                                       <  1:1062993/149‑1 (MQ=255)
                                                                               gATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCa                                                                       <  1:844072/150‑1 (MQ=255)
                                                                                 ttAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGAAACGGTTATCGTCGGTCAGACTATATTGCACCGTCTCGCCGAGATTACCCGGGAAGCCCTGATCACCATGATCTGAACTCAGGGCAAACAGCAcc                                                                     <  1:1233830/150‑1 (MQ=255)
                                                                                  tAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGTATATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCt                                                                    <  2:1068900/150‑1 (MQ=255)
                                                                                    gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGa                                                                  >  1:1337163/1‑150 (MQ=255)
                                                                                    gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGa                                                                  >  2:1527109/1‑150 (MQ=255)
                                                                                    gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGAGCACCATCATCTGAACTCAGGGCAACCAGCAcagga                                                                  >  2:659781/1‑146 (MQ=255)
                                                                                      cATATTCACCGGGCATGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACTCTTATCGTCGGTCAGACGATATTTCACCGTCGCGCCGAGATTAACCGGGAAGCCCAGATCACCATCAACTGAACTCAGGGCAAACCGCACCTGACg                                                                >  1:107979/1‑150 (MQ=255)
                                                                                      cATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACg                                                                >  2:1276637/1‑150 (MQ=255)
                                                                                         atTCACCGGGCAAGGGTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATc                                                             <  1:216273/150‑1 (MQ=255)
                                                                                          ttCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCg                                                            <  1:1261249/150‑1 (MQ=255)
                                                                                           tCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATTTTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGt                                                           >  1:221960/1‑150 (MQ=255)
                                                                                           tCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGt                                                           <  2:1317007/150‑1 (MQ=255)
                                                                                             aCCGGGCAAGGTTTATCAACTGTGGTGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTc                                                         <  2:1192080/150‑1 (MQ=255)
                                                                                              ccGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCt                                                        >  1:1630071/1‑150 (MQ=255)
                                                                                                gggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTgg                                                      <  1:567928/150‑1 (MQ=255)
                                                                                                gggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTgg                                                      <  1:1641062/150‑1 (MQ=255)
                                                                                                   cAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTc                                                   >  1:1181700/1‑150 (MQ=255)
                                                                                                     aGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTcac                                                 <  1:752118/150‑1 (MQ=255)
                                                                                                              aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCa                                        >  2:210542/1‑150 (MQ=255)
                                                                                                              aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCa                                        <  1:480882/150‑1 (MQ=255)
                                                                                                               aCTGTGGCGCGATAAGTAATGGAGATACGTTTATCGTCGGTCCGACGATATTGCACCTTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAg                                       <  2:1203928/150‑1 (MQ=255)
                                                                                                                cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGc                                      <  2:1177225/150‑1 (MQ=255)
                                                                                                                  gtgGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGa                                    >  2:1092605/1‑150 (MQ=255)
                                                                                                                   tgGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGAc                                   <  2:1543461/150‑1 (MQ=255)
                                                                                                                    ggCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTTCACCGTCGAGCCGAGATTACCCGGGAAGCCCTGATGACCATCATCTGAACTCATGGCAAACAGCACCTGACGATCGTTCTAGTTCACAATCTGCCAGCGACg                                  >  1:1100735/1‑150 (MQ=255)
                                                                                                                          aTAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTc                            <  1:491544/150‑1 (MQ=255)
                                                                                                                           aaaGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCg                           <  2:712043/149‑1 (MQ=255)
                                                                                                                            aaGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGa                          <  2:1686432/150‑1 (MQ=255)
                                                                                                                               tAATGGAGATACGGTTATCGTCGGACAGACGATATTGCACCGTATCGCCGAGATTAACCGGGAAGCACTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAAcc                       <  1:282094/150‑1 (MQ=255)
                                                                                                                                 aTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCt                     >  1:628178/1‑150 (MQ=255)
                                                                                                                                     agaTACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCg                 <  1:1634928/150‑1 (MQ=255)
                                                                                                                                       ggtgcaatATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGc               >  1:50592/8‑150 (MQ=255)
                                                                                                                                       aTACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGc               >  2:1470698/1‑150 (MQ=255)
                                                                                                                                        tACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGcc              <  1:1196437/150‑1 (MQ=255)
                                                                                                                                          gcaatATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCCCCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGcccc            >  2:632823/5‑149 (MQ=255)
                                                                                                                                            aatATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGCACTCAGGGCAAACGGCACCTGACGAGCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGCTCCCTTGCGGCccgcc          >  2:240126/3‑150 (MQ=255)
                                                                                                                                                tcgtcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGc      <  1:1630653/150‑1 (MQ=255)
                                                                                                                                                tcgtcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGc      <  1:1527109/150‑1 (MQ=255)
                                                                                                                                                   gcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGAGCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGc   <  1:459588/149‑1 (MQ=255)
                                                                                                                                                    cgGGCAGACGATATTGCCCCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCt  <  1:899417/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCT  >  NC_000913/788236‑788533

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