Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 788,389 A→C 5.6% L150R (CTG→CGG)  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,3890AC 5.6% 290.1 / 5.2 100L150R (CTG→CGG) galMaldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (45/49);  new base C (3/3);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 101 total aligned reads displayed.

GATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGG  >  NC_000913/788240‑788535
                                                                                                                                                     |                                                                                                                                                  
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)
       gTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTGGCCGTCAAGATTGAAATAGGCGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTTATGGAGTTACGGTTATCGTCGGTCAGCCGata                                                                                                                                             >  1:313472/1‑150 (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                                                                                                    >  2:1196437/1‑150 (MQ=255)
                                                gCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                    >  1:593899/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:1641062/150‑1 (MQ=255)
                                                                                            gggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTgg                                                        <  1:567928/150‑1 (MQ=255)
                                                                                               cAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTc                                                     >  1:1181700/1‑150 (MQ=255)
                                                                                                 aGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTcac                                                   <  1:752118/150‑1 (MQ=255)
                                                                                                          aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCa                                          <  1:480882/150‑1 (MQ=255)
                                                                                                          aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCa                                          >  2:210542/1‑150 (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)
                                                                                                                      aTAAGTCATGGAGATACGTTTATCGGAGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACTATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTc                              <  2:1246858/150‑1 (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)
                                                                                                                                                  gTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTgg  <  2:742300/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGG  >  NC_000913/788240‑788535

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: