Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,282,928 T→A 7.2% I318N (ATC→AAC) ‡ 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,9280TA 7.2% 294.2 / 10.8 83I318N (ATC→AAC) ‡agaStagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base A (3/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

AAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGT  >  NC_000913/3282780‑3283076
                                                                                                                                                    |                                                                                                                                                    
aaaTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAt                                                                                                                                                      >  3:364282/1‑149 (MQ=255)
 aaTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtc                                                                                                                                                     <  1:911878/149‑1 (MQ=255)
    cGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcgtc                                                                                                                                                  <  4:364282/149‑1 (MQ=255)
     gCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcgtcg                                                                                                                                                 <  3:1080282/149‑1 (MQ=255)
              tGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTc                                                                                                                                        <  3:1406351/149‑1 (MQ=255)
                aCGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGc                                                                                                                                                <  4:1389722/139‑1 (MQ=255)
                aCGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGc                                                                                                                                                >  3:1389722/1‑139 (MQ=255)
                    aaCGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACata                                                                                                                                  >  1:701810/1‑149 (MQ=255)
                         tggtggtggTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatcatc                                                                                                                             >  3:806276/1‑149 (MQ=255)
                               tggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGc                                                                                                                                                >  4:438545/1‑124 (MQ=255)
                               tggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGc                                                                                                                                                <  3:438545/124‑1 (MQ=255)
                                    tttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                            >  4:823021/1‑122 (MQ=255)
                                    tttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                            <  3:823021/123‑2 (MQ=255)
                                          tcCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgtca                                                                                                            >  1:260025/1‑149 (MQ=255)
                                          tcCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgtca                                                                                                            >  3:534040/1‑149 (MQ=255)
                                             aGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGc                                                                                                                                             <  3:1360081/113‑1 (MQ=255)
                                             aGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGc                                                                                                                                             >  4:1360081/1‑113 (MQ=255)
                                               ccACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTa                                                                                                       >  2:1036648/1‑149 (MQ=255)
                                                 aCCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATc                                                                                                     >  4:1146712/1‑149 (MQ=255)
                                                 aCCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATc                                                                                                     >  1:615269/1‑149 (MQ=255)
                                                 aCCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATc                                                                                                     >  3:75709/1‑149 (MQ=255)
                                                   ccTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTc                                                                                                                                        >  4:836904/1‑112 (MQ=255)
                                                   ccTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTc                                                                                                                                        <  3:836904/112‑1 (MQ=255)
                                                     ttACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCt                                                                                                                                               <  1:1127988/103‑1 (MQ=255)
                                                     ttACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCCTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCt                                                                                                                                               >  2:1127988/1‑103 (MQ=255)
                                                       acacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGtt                                                                                               <  1:1036648/149‑1 (MQ=255)
                                                       acacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGtt                                                                                               >  1:448633/1‑149 (MQ=255)
                                                         acCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGa                                                                                             <  3:1146712/149‑1 (MQ=255)
                                                          cccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAg                                                                                            <  2:615269/149‑1 (MQ=255)
                                                           ccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCTGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatcatc                                                                                                                             <  2:578125/115‑1 (MQ=255)
                                                           ccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCTGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatcatc                                                                                                                             >  1:578125/1‑115 (MQ=255)
                                                           ccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAGTCGCCATCGCCGCGGAAAGCAGCGACATCGTCACTGCCGGTcc                                                                                                                                       >  4:1396036/1‑105 (MQ=255)
                                                           ccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTcc                                                                                                                                       <  3:1396036/105‑1 (MQ=255)
                                                           ccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatcatc                                                                                                                             <  2:78526/115‑1 (MQ=255)
                                                           ccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatcatc                                                                                                                             >  1:78526/1‑115 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGcc                                                                                                                                            >  4:395870/1‑98 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGcc                                                                                                                                            <  3:395870/98‑1 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGc                                                                                                                                             <  2:1357033/97‑1 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGc                                                                                                                                             <  2:1304115/97‑1 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGc                                                                                                                                             >  1:1357033/1‑97 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGc                                                                                                                                             >  1:1304115/1‑97 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAg                                                                                         <  2:448633/149‑1 (MQ=255)
                                                              tCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAgg                                                                                        >  4:619537/1‑149 (MQ=255)
                                                                 gTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcac                                                                                                                <  3:948940/122‑1 (MQ=255)
                                                                 gTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcac                                                                                                                >  4:948940/1‑122 (MQ=255)
                                                                  tATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAtt                                                                                    >  2:377585/1‑149 (MQ=255)
                                                                      atcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGc                                                                                >  2:461741/1‑149 (MQ=255)
                                                                       tcttgatctGCTGGCTGAACTTCGCCGTGACAACTAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGc                                                                                                                                             <  4:748892/87‑1 (MQ=255)
                                                                       tcttgatctGCTGGCTGAACTTCGCCGTGACAACTAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGc                                                                                                                                             >  3:748892/1‑87 (MQ=255)
                                                                       tcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTAt                                                                                                      <  4:222775/126‑1 (MQ=255)
                                                                       tcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTAt                                                                                                      >  3:222775/1‑126 (MQ=255)
                                                                         ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                             <  4:534040/149‑1 (MQ=255)
                                                                         ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                             <  4:75709/149‑1 (MQ=255)
                                                                              ctgctgGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGAt                                                                        <  1:377585/149‑1 (MQ=255)
                                                                              ctgctgGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGAt                                                                        <  1:461741/149‑1 (MQ=255)
                                                                               tgctgGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATg                                                                       >  2:1374302/1‑149 (MQ=255)
                                                                                   ggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGa                                                                                             >  3:442726/1‑123 (MQ=255)
                                                                                   ggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGa                                                                                             <  4:442726/123‑1 (MQ=255)
                                                                                    gcagAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTAGGc                                                                  >  3:970243/4‑149 (MQ=255)
                                                                                      tgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGc                                                                                <  1:14101/133‑1 (MQ=255)
                                                                                      tgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGc                                                                                >  2:14101/1‑133 (MQ=255)
                                                                                       gAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCa                                                               >  3:55313/1‑149 (MQ=255)
                                                                                             cGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGtt                                                         <  3:1084020/149‑1 (MQ=255)
                                                                                               ccGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTg                                                       <  4:55313/149‑1 (MQ=255)
                                                                                                   gACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAtttt                                                                                  <  2:48749/118‑1 (MQ=255)
                                                                                                   gACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAtttt                                                                                  >  1:48749/1‑118 (MQ=255)
                                                                                                      aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                                                                                              <  3:33507/87‑1 (MQ=255)
                                                                                                      aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                                                                                              >  4:33507/1‑87 (MQ=255)
                                                                                                      aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTAt                                                                                                      >  2:920096/1‑95 (MQ=255)
                                                                                                      aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTAt                                                                                                      <  1:920096/95‑1 (MQ=255)
                                                                                                      aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGAt                                                <  4:970243/149‑1 (MQ=255)
                                                                                                      aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGAt                                                >  1:1264828/1‑149 (MQ=255)
                                                                                                       aCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATg                                               <  4:806276/149‑1 (MQ=255)
                                                                                                       aCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATg                                               <  2:260025/149‑1 (MQ=255)
                                                                                                        ccAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGc                                              >  4:194146/1‑149 (MQ=255)
                                                                                                         cAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                             <  1:1465001/117‑1 (MQ=255)
                                                                                                         cAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                             >  2:1465001/1‑117 (MQ=255)
                                                                                                           ggCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGt                                           >  2:602890/1‑149 (MQ=255)
                                                                                                                      tAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATg                                <  1:1374302/149‑1 (MQ=255)
                                                                                                                       aaTCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATgg                               <  3:619537/149‑1 (MQ=255)
                                                                                                                                          agcagcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCaa                           >  3:1232527/1‑134 (MQ=255)
                                                                                                                                          agcagcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCaa                           <  4:1232527/134‑1 (MQ=255)
                                                                                                                                             agcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATc         >  4:691549/1‑149 (MQ=255)
                                                                                                                                              gcGACATCGTCGCTGCCGGTCCACATATCATCCTGCTACCGTCACGTCACTTTATCGACGTTGAGCAGGCAttt                                                                                   >  2:396014/1‑74 (MQ=255)
                                                                                                                                              gcGACATCGTCGCTGCCGGTCCACATATCATCCTGCTACCGTCACGTCACTTTATCGACGTTGAGCAGGCAttt                                                                                   <  1:396014/74‑1 (MQ=255)
                                                                                                                                                   aTCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCAt          <  2:832486/142‑1 (MQ=255)
                                                                                                                                                   aTCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCAt          >  1:832486/1‑142 (MQ=255)
                                                                                                                                                    tcgtcgCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCa                                                     >  2:1097555/1‑98 (MQ=255)
                                                                                                                                                    tcgtcgCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCa                                                     <  1:1097555/98‑1 (MQ=255)
                                                                                                                                                    tcgtcgCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGt  >  3:399145/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGT  >  NC_000913/3282780‑3283076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: