Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,515 T→A 8.7% intergenic (+70/+17) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,5150TA 8.7% 1523.4 / 64.5 525intergenic (+70/+17)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (165/272);  new base A (41/12);  total (239/286)
Fisher's exact test for biased strand distribution p-value = 4.06e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 537 total aligned reads displayed.

GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCT  >  NC_000913/3959366‑3959664
                                                                                                                                                     |                                                                                                                                                     
gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                                       <  2:5142653/150‑1 (MQ=255)
gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                                       <  1:7790644/150‑1 (MQ=255)
   gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                                    <  1:1639362/150‑1 (MQ=255)
   gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                                    <  2:9290495/150‑1 (MQ=255)
      aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                 <  2:4183057/150‑1 (MQ=255)
       tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                                <  1:8631704/150‑1 (MQ=255)
        gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                                                                               <  2:7557815/150‑1 (MQ=255)
        gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAAAACTCAAGGCCATCTCCAGGAGAAgg                                                                                                                                               >  1:7138149/1‑150 (MQ=255)
           aTTGAGCAGGTCGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACACCCCCATGGGGGGGGAGTAAACACTCAAAGCCCTCAACGGGAGAAGGccg                                                                                                                                            >  2:3836947/1‑149 (MQ=255)
             tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                          <  1:2713628/150‑1 (MQ=255)
               aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAATCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCAAACCCACAGGGAGAGGGGGAAAACACTCAAGGCCATCTCCGGGCGAAGGCCTTGCt                                                                                                                                        >  2:94686/1‑150 (MQ=255)
               aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCACACCCCCAGGGGGAGGGGGAAAACACTCAACGCCGTCTCCAGGAGAAGGccattcg                                                                                                                                        >  1:5544263/1‑145 (MQ=255)
                 cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                      <  1:2168856/150‑1 (MQ=255)
                 cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                      >  2:4949005/1‑150 (MQ=255)
                  aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                                                                                                                                     <  1:1976958/150‑1 (MQ=255)
                   ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                                    <  1:8193504/150‑1 (MQ=255)
                   ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGCCTTgagtgtt                                                                                                                                    <  2:8889175/150‑7 (MQ=255)
                      aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTATCCTGGAGAAGGCCTTGCTATTAGtt                                                                                                                                 >  2:1545540/1‑150 (MQ=255)
                        gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                               <  1:2275726/150‑1 (MQ=255)
                        gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCGCACCCCCAGGGCGAGGGAGAAAACCCGCAAGGCCGTCGCCGGGCGAAGGCCGGGCTAGTAGGTGc                                                                                                                               >  1:1559922/1‑150 (MQ=255)
                          aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGTGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                             <  1:798952/150‑1 (MQ=255)
                          aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                             <  1:3821154/150‑1 (MQ=255)
                          aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                             <  1:5119209/150‑1 (MQ=255)
                           aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                            <  2:4301769/150‑1 (MQ=255)
                           aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGATAAGTGCGCGCAGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                            <  2:456620/150‑1 (MQ=255)
                            gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCCACCCACAACCACACGGGGAGGGGGCAAAAACAAAAGGCCGTATACAGGAGCAGGGCGTGCGATTAttttgggtt                                                                                                                           >  1:7640135/1‑144 (MQ=255)
                            gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACTCCCACAGGGGGTGGGGGAAAAAACCCAAGGCCATGGCCGGGGGGAGGCCGGGGTATTATTTGCGGtg                                                                                                                           >  2:3730978/1‑149 (MQ=255)
                                cggcgcGGCTATATGACAGCTATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCTCCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                                                                                                       <  2:8184761/148‑1 (MQ=255)
                                   cgcgGCTATATGACTGATATGAAACGTTTTGCTGTTGGGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGCGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                                    <  2:7218600/150‑1 (MQ=38)
                                    gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACACGGAGGGGGCGAAAAAACACAAGGCCATCAACAGGAGAAGGCCCTGGGATTAGTTTCGGGACAGCAcca                                                                                                                   >  1:10031777/1‑149 (MQ=255)
                                      ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTCATTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGCGAGGGAGAAAACACACAAGGCCGTCACCTGGAGCAGGCCTTGCTATAAGTTGCGGTACAGCACCTTa                                                                                                                 >  2:6510416/1‑150 (MQ=38)
                                       gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                                                                <  2:4781788/150‑1 (MQ=255)
                                        cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                                               <  2:2349136/150‑1 (MQ=255)
                                        cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                                               <  1:4468577/150‑1 (MQ=255)
                                        cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                                               <  1:8353954/150‑1 (MQ=255)
                                           tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                                            <  2:8405507/150‑1 (MQ=255)
                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                          <  2:3177421/150‑1 (MQ=255)
                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                          <  1:3515845/150‑1 (MQ=255)
                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCAAACCCCCACCCACAACCACAGGGCGAGGGTCTAAACACAAAAGGCCTTCTCCGGGAGAAAGCCATGCTAATATTTGCGGTACAGCACCATAatgattg                                                                                                          >  2:553993/1‑148 (MQ=255)
                                              gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCCCCCACACCCACAGGGAGAGGGAGAAAAAACACAAGGCCATCGCCGGGAGAAGGCCGTGCTAATAGTTGCGGTACAGCACCATAAATAtgtg                                                                                                         >  2:9785968/1‑150 (MQ=25)
                                                  gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                                                                                                     <  1:6871401/150‑1 (MQ=255)
                                                    tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCATCTACCACAGGGAGAGGGAGAAAAAACACAAGGCCATATACTTGCGCAGGCCGGGCTATTAGTTGCGGTAAAGCACCGTAATGATGTGatatac                                                                                                   >  2:8082092/1‑148 (MQ=25)
                                                      tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACTCCCACACGGAGAGGGAGAAAAAACTCAAGGCCATATACAGGAGAAGGCCGTGCTATTAGTTTCGGTACACCACCATAATTATTTGatataaga                                                                                                 >  2:3612702/1‑146 (MQ=255)
                                                        aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTACTGGAGAAGGCCTTGCTATTATTTGCGGTACAGCACCTTAATGTTGTGATATACGAAc                                                                                               >  2:6249363/1‑150 (MQ=255)
                                                        aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCAAACCCACAGGGAGAGGGAGAAAAAAAAAAAGGCCTTATACAGGGGGAGGCCGTGGTATTATTTGCGGTACAGCACCATAATTATGTGATATACGAAc                                                                                               >  1:9129663/1‑150 (MQ=21)
                                                           cGTTTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCTCCCCGACCCTCTCCCACAGGGAGAGTGAGATTACACTCTAGGCCTTCTCCTGGAGTAGGGCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                            <  2:7886298/150‑1 (MQ=255)
                                                           cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                            <  2:8633668/150‑1 (MQ=255)
                                                             tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                          <  1:7413728/150‑1 (MQ=255)
                                                                tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg                                                                                       <  1:3742207/150‑1 (MQ=255)
                                                                tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCCCCCCCACCCACACCCACAGGGCGTGGGGGAAAAAACACAAGGCCTTCATCGGGGGAAGGCCATGCTCTTAGTTGCGGGACAGCACCATAATGGTTTGATATCCGCACTGGGtgtg                                                                                       >  2:3042865/1‑150 (MQ=255)
                                                                  cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCa                                                                                     <  1:7272819/150‑1 (MQ=255)
                                                                   tGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCCAACCAAAACCCCACGGAGAGGGTGAAAAACCACAATGCCATATCCGTGAGAAAGCCATGCTATTATTTTCGGTACAGCCCCCTAATGAATTGAAAAACCAACTTGGTGTGCAg                                                                                    >  2:9335828/1‑150 (MQ=255)
                                                                     ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCg                                                                                  <  2:6747501/150‑1 (MQ=255)
                                                                       gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg                                                                                <  2:8941903/150‑1 (MQ=255)
                                                                       gCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCACAGGGAGAGGGAGAAAAAACAAAAGGCCATCACCTGGAGAATGCCTTGCTATTATTTTCGGTACATCACCATAATGATGTGATATCCGAACATGGTGTGCAGCggg                                                                                >  2:6008437/1‑150 (MQ=25)
                                                                         gggTGAAGTGCGCGCTGATTCCCCCAACCCGACCCACAACCACAGGGAGAGGGAGAAAAAACTCAAGGCCATCTCCTGGAGAAGGCCGTGCAATTACTTTCGGTACAGCACCTTAAAGATTTGTACAACGAACAGGGTGTGCGCCGGGcc                                                                              >  2:8735362/1‑150 (MQ=255)
                                                                             tAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                                          <  2:196630/150‑1 (MQ=255)
                                                                                  gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                                     <  2:66171/150‑1 (MQ=255)
                                                                                   cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                                                    <  2:9937603/150‑1 (MQ=255)
                                                                                      gcTGATGCCCTCAACCCCACCCTCACCCACACGGAGCGGGAGAAAAAACTCAAGGCCTTCACCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCATAAAGATGTGATATCCCAACTTGGTGTGCAGCGGGCCCGTCGGCACCacg                                                                 >  2:7447411/1‑148 (MQ=255)
                                                                                       cTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                                                                <  1:5807658/150‑1 (MQ=255)
                                                                                       cTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                                                                <  2:683794/150‑1 (MQ=255)
                                                                                       cTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                                                                <  1:5208792/150‑1 (MQ=255)
                                                                                            gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                                                           <  1:476412/150‑1 (MQ=255)
                                                                                               cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                                        <  1:3317727/150‑1 (MQ=255)
                                                                                                 cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACGGGACAaga                                                      <  2:7640135/150‑1 (MQ=255)
                                                                                                   ccccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                    <  2:111135/150‑1 (MQ=255)
                                                                                                     ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaa                                                  >  1:6272383/1‑150 (MQ=255)
                                                                                                      cGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                                 <  2:1438943/150‑1 (MQ=255)
                                                                                                       gACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                                <  1:6622036/150‑1 (MQ=255)
                                                                                                          ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGACAACCAc                                             >  2:4638200/1‑150 (MQ=255)
                                                                                                            tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                                           >  1:9382521/1‑150 (MQ=255)
                                                                                                            tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                                           >  1:7483526/1‑150 (MQ=255)
                                                                                                                  acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGa                                     >  1:88834/1‑150 (MQ=255)
                                                                                                                      ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc                                 >  1:1305093/1‑150 (MQ=255)
                                                                                                                       gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                                >  2:4010350/1‑150 (MQ=255)
                                                                                                                       gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                                >  1:5434119/1‑150 (MQ=255)
                                                                                                                        tgagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGTTACAGCACCTTAATGATGTGATATCCTAAGTGGGTGTGCAGGGGGCGGGTCGGATCCATTACCGGAAAAGATAAAACCACTTTATCGAACGCCg                               <  1:9228941/149‑1 (MQ=255)
                                                                                                                         gtgggtgaAAACACTCCCGGCCTTCTCCTGGAGATGGCCTTGCTATTTGTTGCGGTACAGCACCTTATTGTTGTTATATCCTAACTGGGTGTGCAGGGGGCCGGTCGGCTCCAGTACCGTACAAGATAAAACCACTTTATCGAACGCCgg                              <  1:9899210/148‑1 (MQ=255)
                                                                                                                         gtgggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGTTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGATAAAACCACTTTATCGAACGCCgg                              <  1:2621019/148‑1 (MQ=255)
                                                                                                                             gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAcc                          >  2:5265200/1‑150 (MQ=255)
                                                                                                                              agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCa                         <  2:2101202/150‑1 (MQ=255)
                                                                                                                                aaaaCAATCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTATTGATGTGATATCCGAACTGGGTGTGCTGCGGGCCGGTCGGCTCCAGTACCGGACAAGTGAAAACCACTTTATCGAACGCCGGAACCATc                       <  1:8521816/150‑1 (MQ=255)
                                                                                                                                 aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAACACCACTTTATCGAACGCCGGCACCAGCt                      >  1:6408600/1‑150 (MQ=255)
                                                                                                                                 aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCt                      <  2:545406/150‑1 (MQ=255)
                                                                                                                                 aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCt                      <  2:3379530/150‑1 (MQ=255)
                                                                                                                                  aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTg                     >  1:2382393/1‑150 (MQ=255)
                                                                                                                                   acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                    >  2:2728215/1‑150 (MQ=255)
                                                                                                                                       tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt                >  2:6999758/1‑150 (MQ=255)
                                                                                                                                        cAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCg                               <  2:2193512/134‑1 (MQ=255)
                                                                                                                                         aaGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTgg              >  2:7373344/1‑150 (MQ=255)
                                                                                                                                         aaGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTgg              <  1:1190992/150‑1 (MQ=255)
                                                                                                                                         aaGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAAGACCACTTTATCGAACGCCGGAACCATCTGACCCTgg              >  2:1812574/1‑150 (MQ=255)
                                                                                                                                          aGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGc             <  1:9785968/150‑1 (MQ=255)
                                                                                                                                          aGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGc             >  1:9515594/1‑150 (MQ=255)
                                                                                                                                             ccTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGa          <  1:6007834/150‑1 (MQ=255)
                                                                                                                                                  tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCa     <  1:7705906/150‑1 (MQ=255)
                                                                                                                                                  tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCa     <  1:3380544/150‑1 (MQ=255)
                                                                                                                                                  tcCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCa     >  2:5943826/1‑150 (MQ=255)
                                                                                                                                                     tGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCt  >  1:6780643/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCT  >  NC_000913/3959366‑3959664

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: