Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,282,927 A→G 5.2% I318V (ATC→GTC) ‡ agaS → tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,282,9270AG 5.2% 351.5 / 6.8 115I318V (ATC→GTC) ‡agaStagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base A (66/43);  new base G (2/4);  total (68/47)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 117 total aligned reads displayed.

CAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCC  >  NC_000913/3282778‑3283073
                                                                                                                                                     |                                                                                                                                                  
cAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACa                                                                                                                                                    >  2:375811/1‑150 (MQ=255)
   aaTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcg                                                                                                                                                 >  1:1313094/1‑150 (MQ=255)
    aTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcgt                                                                                                                                                <  2:1042013/150‑1 (MQ=255)
       gCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGc                                                                                                                                             >  2:938181/1‑150 (MQ=255)
          ggTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGc                                                                                                                                          <  2:1171594/150‑1 (MQ=255)
          ggTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGc                                                                                                                                          <  2:1165548/150‑1 (MQ=255)
           gTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                         <  1:1761776/150‑2 (MQ=255)
           gTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                         <  1:14549/150‑2 (MQ=255)
            tCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTTCTGGCTGACCTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCg                                                                                                                                        >  1:1044944/1‑150 (MQ=255)
            tCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCg                                                                                                                                        >  2:1738162/1‑150 (MQ=255)
             cGATGACGAAACGCTGGTGGTGGTCTTTGTCACCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCgg                                                                                                                                       <  1:1098724/150‑1 (MQ=255)
             cGATGACGAAACGCTGGTGGTGGTATTTGTCTCCATCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGATGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAACGCAGCGACATCGTCGCTGCCgg                                                                                                                                       >  2:954019/1‑150 (MQ=255)
               aTGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTc                                                                                                                                     >  1:742/1‑150 (MQ=255)
                 gACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCca                                                                                                                                   <  1:692263/150‑1 (MQ=255)
                   cGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCcaca                                                                                                                                 <  1:1800584/150‑1 (MQ=255)
                    gAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACat                                                                                                                                <  1:493073/150‑1 (MQ=255)
                     aaaCGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACata                                                                                                                               >  1:1301242/1‑150 (MQ=255)
                       aCGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCCGGCAATGCGTGTAATCGCCATCGCCGCGGAACGCAGCGACATCGTCGCTGCCGGTCCACATatc                                                                                                                             <  2:1408394/150‑1 (MQ=255)
                           tggtggtggTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATcc                                                                                                                         >  2:172228/1‑150 (MQ=255)
                            ggtgggggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGCTCTGCTGGCTGACCTTCGCCGTGACAACCAGGCAATGCGTTTAATCGCCATCGCCGCGGAAAGCACCGACATCGTCGCTGCCGGTCCACATATCATCCt                                                                                                                        <  2:262228/150‑1 (MQ=255)
                                 tggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACACCCAGGCAATGCGAGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGACGGTCCACATATCATCCTGCCAc                                                                                                                   >  2:1709981/1‑150 (MQ=255)
                                    tATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCCCcgt                                                                                                                >  1:1110688/1‑150 (MQ=255)
                                     aTTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtc                                                                                                               <  2:1232343/150‑1 (MQ=255)
                                       ttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcac                                                                                                             >  2:223129/1‑150 (MQ=255)
                                        tGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacg                                                                                                            >  2:1441400/1‑150 (MQ=255)
                                         gTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                                                                                           <  2:1085830/150‑1 (MQ=255)
                                         gTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                                                                                           >  2:118234/1‑150 (MQ=255)
                                          tctcCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgtc                                                                                                          <  2:1071541/150‑1 (MQ=255)
                                           ctcCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCCGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgtca                                                                                                         >  2:710835/1‑150 (MQ=255)
                                             ccAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACt                                                                                                       >  2:504451/1‑150 (MQ=255)
                                               aGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACttt                                                                                                     >  1:837647/1‑150 (MQ=255)
                                                  cACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCGGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATc                                                                                                  <  2:1816582/150‑1 (MQ=255)
                                                   aCCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCg                                                                                                 >  2:834860/1‑150 (MQ=255)
                                                       ttACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGt                                                                                             <  2:1313094/150‑1 (MQ=255)
                                                        tACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGtt                                                                                            >  2:1744277/1‑150 (MQ=255)
                                                         acacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTg                                                                                           <  2:471450/150‑1 (MQ=255)
                                                          cacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGa                                                                                          >  1:901444/1‑150 (MQ=255)
                                                            cccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGc                                                                                        >  1:1794823/1‑150 (MQ=255)
                                                            cccGTCAGTATGATCTTGATCTGCTGGCAGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCTTCCTGCCACCGTCACGTCACTTTATCGACGTTGAGc                                                                                        >  1:153940/1‑150 (MQ=255)
                                                                    tATGATCTTGATCTGCTGTCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAttt                                                                                <  2:689270/150‑1 (MQ=255)
                                                                    tATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAttt                                                                                <  1:434396/150‑1 (MQ=255)
                                                                       gatcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCCCCGTCACGTCCCTTTATCGACGTTGAGCAGGCATTTTGc                                                                             >  1:1344585/1‑150 (MQ=255)
                                                                       gatcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGCGCAGGCATTTTGc                                                                             >  2:84500/1‑150 (MQ=255)
                                                                         tcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAACGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCtt                                                                           >  2:343000/1‑150 (MQ=255)
                                                                         tcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCtt                                                                           <  1:613155/150‑1 (MQ=255)
                                                                          cttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGTGACATCGTCGCTGCCGGTCCAAATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                          >  2:827063/1‑150 (MQ=255)
                                                                          cttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                          >  2:1470844/1‑150 (MQ=255)
                                                                          cttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                          >  2:376799/1‑150 (MQ=255)
                                                                          cttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATATTGCTTc                                                                          >  2:1405084/1‑150 (MQ=255)
                                                                           ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTcc                                                                         <  1:759158/150‑1 (MQ=255)
                                                                            agatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGATGCCGGTCCACATATCATCCTGCCACCGTCAAGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCt                                                                        >  2:1732394/2‑150 (MQ=255)
                                                                             gatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTg                                                                       <  2:938679/150‑1 (MQ=255)
                                                                              atctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGa                                                                      <  1:1540885/150‑1 (MQ=255)
                                                                               tctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGAt                                                                     <  1:994302/150‑1 (MQ=255)
                                                                                ctgctgGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATg                                                                    <  2:1150129/150‑1 (MQ=255)
                                                                                ctgctgGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATg                                                                    >  2:149399/1‑150 (MQ=255)
                                                                                  gctggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTa                                                                  <  1:472757/150‑1 (MQ=255)
                                                                                   ctggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTACACATATCATCCTGCCCCCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTAc                                                                 >  1:373855/1‑150 (MQ=255)
                                                                                      gctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGcc                                                              <  1:554547/150‑1 (MQ=255)
                                                                                       ctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGccc                                                             >  2:1782981/1‑150 (MQ=255)
                                                                                       ctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGccc                                                             >  1:1583822/1‑150 (MQ=255)
                                                                                           aCTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGAc                                                         <  1:709624/150‑1 (MQ=255)
                                                                                            cTTCGCCGTGACAACCCGGCAATGCGGGTAATCGCCATCGCCGCGGAACGCCGCGACATCGTCGCGGCCGGTCCACATATCATCCTGCCCCCGTCCCGTCACTGTATCGACGGTGAGCAGGCAGTTTGCTTCCTGATGTACGACCAGACg                                                        >  1:267341/1‑150 (MQ=255)
                                                                                             ttCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGt                                                       >  2:1034661/1‑150 (MQ=255)
                                                                                               cGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGACAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGttt                                                     <  2:282875/150‑1 (MQ=255)
                                                                                                gCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTGTATCGACGTTGAGCAGGGATTTTGCGTACTGATGTACCCCCAGCCGTTTg                                                    >  2:907056/1‑150 (MQ=255)
                                                                                                gCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTGTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTg                                                    >  1:360959/1‑150 (MQ=255)
                                                                                                 ccGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGc                                                   >  2:805535/1‑150 (MQ=255)
                                                                                                   gTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCAc                                                 <  1:1232819/150‑1 (MQ=255)
                                                                                                   gTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCAc                                                 >  1:1701043/1‑150 (MQ=255)
                                                                                                    tGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACt                                                >  1:1343018/1‑150 (MQ=255)
                                                                                                    tGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACt                                                <  1:1425885/150‑1 (MQ=255)
                                                                                                     gACAACCTGGCAATGCGTGGAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTg                                               <  2:1602324/150‑1 (MQ=255)
                                                                                                       cAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGAt                                             >  1:1106781/1‑150 (MQ=255)
                                                                                                          ccAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCa                                          >  1:22442/1‑150 (MQ=255)
                                                                                                           cAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAg                                         >  1:1162469/1‑150 (MQ=255)
                                                                                                              gCAATGCGTGTAAACGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCTGTCCACATATCTTCGTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCg                                      <  1:118234/150‑1 (MQ=255)
                                                                                                               cAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGc                                     <  1:375811/150‑1 (MQ=255)
                                                                                                                aaTGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCt                                    >  2:369115/1‑150 (MQ=255)
                                                                                                                aaTGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCt                                    <  1:402634/150‑1 (MQ=255)
                                                                                                                  tGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGc                                  <  2:1162469/150‑1 (MQ=255)
                                                                                                                    cGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGcac                                >  1:1103034/1‑150 (MQ=255)
                                                                                                                         aaTCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATggg                           >  1:1559461/1‑150 (MQ=255)
                                                                                                                          atcgccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGc                          >  1:642828/1‑150 (MQ=255)
                                                                                                                           tcgccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGATGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTTCACTGATGCAGTCGCTGCACATGGGCa                         >  2:1207233/1‑150 (MQ=255)
                                                                                                                             gccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAAt                       >  1:1432099/1‑150 (MQ=255)
                                                                                                                                  cgccgcGGAAAGCAGCGACATCGTCGCTGCCTGTCCACATATCATCCAGCCCCCGTCACCTCACATTATCGACGTTGAGCAGGCAATTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACAAGGGAAATACGcc                  >  1:182162/1‑150 (MQ=255)
                                                                                                                                   gccgcGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCg                 <  2:500137/150‑1 (MQ=255)
                                                                                                                                       cgGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATa             <  2:901444/150‑1 (MQ=255)
                                                                                                                                        ggAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATAc            <  1:1694519/150‑1 (MQ=255)
                                                                                                                                        ggAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATAc            <  2:1745104/150‑1 (MQ=255)
                                                                                                                                        ggAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGATTGCACTGATGCAGACGCTGCACAAGGGCAATACGCCGGATAc            >  1:1110468/1‑150 (MQ=255)
                                                                                                                                         gAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATAcc           >  1:1228085/1‑150 (MQ=255)
                                                                                                                                           cggcagcGACGATGTCGCTGCTGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATAcccc         >  1:199506/3‑150 (MQ=255)
                                                                                                                                            ggcagcGACGATGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCa        >  2:1521872/2‑150 (MQ=255)
                                                                                                                                              cagcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATc      >  1:1458941/1‑150 (MQ=255)
                                                                                                                                              cagcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATc      >  2:15292/1‑150 (MQ=255)
                                                                                                                                               agcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCCCTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGCTGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCa     >  2:1482515/1‑150 (MQ=255)
                                                                                                                                                gcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAg    >  2:1033748/1‑150 (MQ=255)
                                                                                                                                                  gACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGcc  >  1:661983/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCC  >  NC_000913/3282778‑3283073

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: