Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,515 T→A 5.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 5.7% 348.4 / 7.4 129intergenic (+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 (56/60);  new base A (3/5);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 129 total aligned reads displayed.

CAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACC  >  NC_000913/3959367‑3959663
                                                                                                                                                    |                                                                                                                                                    
 aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCTCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCACTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTACTCCTga                                                                                                                                                    >  1:236632/1‑149 (MQ=255)
     aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTATACCACAGGGAGAGGGAGAAAACACTAAATGCCATATCATGGAGaa                                                                                                                                                >  1:488533/1‑150 (MQ=38)
      tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                               <  1:65818/150‑1 (MQ=255)
        cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                                             >  2:1319178/1‑150 (MQ=255)
          aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                           <  1:738971/150‑1 (MQ=255)
           ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGCTAGGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCGGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCCAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                                                          <  1:1261307/150‑1 (MQ=38)
           ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCtt                                                                                                                                          <  1:1603348/150‑1 (MQ=255)
           ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGCCtt                                                                                                                                          <  1:494231/150‑1 (MQ=255)
            tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTg                                                                                                                                         <  1:227894/150‑1 (MQ=255)
             gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                                                        >  2:734738/1‑150 (MQ=255)
             gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                                                        >  2:1336073/1‑150 (MQ=255)
              aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                                                       <  1:1064053/150‑1 (MQ=255)
               gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                                      <  2:1528661/150‑1 (MQ=255)
                cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                     <  1:1324638/150‑1 (MQ=255)
                 aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGTGAAAAAACTCAAGGCCTTCTCCTGGATAAGGCCGTGCTAtt                                                                                                                                    >  1:1704469/1‑150 (MQ=255)
                 aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAAACCTCAAGGCCTTCTCCGGGAGAAGGCCTTGCtttt                                                                                                                                    >  2:1166241/1‑150 (MQ=255)
                      ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGTGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                               <  2:1294064/150‑1 (MQ=255)
                      ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                               <  1:1076019/150‑1 (MQ=255)
                       gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                              <  1:1549691/150‑1 (MQ=255)
                       gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTAAACTGGAGAAGGCCTGGGTATTAGTTGc                                                                                                                              >  2:1835820/1‑150 (MQ=38)
                        taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                             <  1:1050847/150‑1 (MQ=255)
                        taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCAACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                             >  1:48893/1‑150 (MQ=255)
                         aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                            <  1:1639100/150‑1 (MQ=255)
                         aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                            <  1:1493439/150‑1 (MQ=255)
                          aGAAACGGCGCGGCTAAATGACAGATCTGAACCGGATTGCGGGGGCGGGTTAAGGGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                           <  2:137684/150‑1 (MQ=37)
                               cTGCGCGGCTATATGACAGATATGATACGTATTGCTGTTTCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGGCGTGCTATTAGTTGCGGTACAGc                                                                                                                      >  2:106580/1‑150 (MQ=255)
                                    cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACGGCACCtt                                                                                                                 <  2:1181152/150‑1 (MQ=255)
                                        tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatg                                                                                                             <  1:675640/150‑1 (MQ=255)
                                         atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatga                                                                                                            <  2:1825405/150‑1 (MQ=255)
                                           atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                                          >  1:1315748/1‑150 (MQ=255)
                                             gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                                        <  2:671329/150‑1 (MQ=255)
                                              aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                                                       <  1:392010/150‑1 (MQ=255)
                                              aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                                                       <  2:1551349/150‑1 (MQ=255)
                                              aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCGGGAGAAGGCCTTGCTATTAGTTGCGGGACAGCACCGTAATGAtgtgc                                                                                                       >  2:420985/1‑149 (MQ=255)
                                                aGATATGAAACGTTTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCACGACCATCTCCCACAGGGAGAGGGAGAAATCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata                                                                                                     <  1:301546/150‑1 (MQ=255)
                                                    atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAACGACGCACAACGGGGTAGGGGTCAAACACAAAAGGCCTTCTCCGGGTGAAAGGCTTTGTTTTATTTTCGGTTCAGGAACGTATATATGTGATATcac                                                                                                 >  2:1412975/1‑148 (MQ=255)
                                                      gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGTGAAAACACTCAAGGCCTTCTCCTGGAGCAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATACGaa                                                                                               >  2:140350/1‑150 (MQ=255)
                                                      gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTATCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGACAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATACGaa                                                                                               >  1:889441/1‑150 (MQ=255)
                                                             aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                                                                                        >  1:810103/1‑150 (MQ=255)
                                                             aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                                                                                        <  1:1154228/150‑1 (MQ=255)
                                                                gcggTTGCGGGTAAAGTGCGCGCGGATGCCCTCACCCCGAGCCTCTGCGACAGGGTGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATCATGTGATATCCGATCTGGGTGTGc                                                                                     <  2:1375140/147‑1 (MQ=255)
                                                                gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGc                                                                                     >  1:1312653/1‑150 (MQ=255)
                                                                gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGc                                                                                     <  1:734738/150‑1 (MQ=255)
                                                                gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGc                                                                                     <  2:35508/150‑1 (MQ=255)
                                                                gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTAAAGCACATTTATTATTTGATATACGAACTGGgtgttc                                                                                     >  2:1000301/1‑148 (MQ=37)
                                                                   gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGc                                                                                  <  1:1857301/150‑1 (MQ=255)
                                                                     tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                                                                                >  2:132469/1‑150 (MQ=255)
                                                                         ggTTAAGTGCGCGCTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                                                            >  2:755317/1‑150 (MQ=255)
                                                                          gTTAAGTGCTCGCTGATGCCCTCACCCCGACCCTCACCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg                                                                           <  2:1819834/150‑1 (MQ=255)
                                                                          gTTAAGTGCGCGCTGATGCCCTCACCCCGAGCCTCTCCCACTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTAGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg                                                                           <  2:523579/150‑1 (MQ=255)
                                                                               gggcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGc                                                                      <  2:1174052/148‑1 (MQ=255)
                                                                                agcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCACCGGGCCGGTCgcca                                                                     >  2:50458/2‑147 (MQ=255)
                                                                                  cgcgcgGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTACTGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                                                   <  1:19634/150‑1 (MQ=255)
                                                                                  cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCCTTCTCCTGGAGAAGGCCTTGCTTTTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                                                   <  1:143153/150‑1 (MQ=255)
                                                                                  cgcgcTGATGCCCGCACCCCGACACTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCGGGTGACGGCCTGGCTATTAGTTGCGGTACACCACCTTAAAGATGTGAAATCCGAACTTGGTGTGCAGCGGGCCGGTCGGCtgc                                                                   >  2:1065460/1‑148 (MQ=37)
                                                                                   gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGATGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                                  <  1:169378/150‑1 (MQ=255)
                                                                                   gcgcTGAGGCCCTCACCCCGACCCTCTCCCACAGCGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                                  <  1:377735/150‑1 (MQ=255)
                                                                                        gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                                                             >  2:1045135/1‑150 (MQ=255)
                                                                                          tccccTCACCGCGACCCTCTCCGACATGGAGAGGGAGAAAACTCTCAAGGCATTCTCCTGGAGAAGGCCTTGCTATAAGTTGCGGTACTGAACCTTAATGATGTGATATCCGTACTGTGTGTGCAGCGGGCCGGTCGGCTCCAGTAGCgg                                                           <  1:1311494/148‑1 (MQ=255)
                                                                                          tGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCgg                                                           <  1:417780/150‑1 (MQ=255)
                                                                                           gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGCAGGCCTTGCTATTAGTTGCGGTACAGCACCGTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAACggc                                                          >  2:1843162/1‑149 (MQ=255)
                                                                                           gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                                                          >  2:1254366/1‑150 (MQ=255)
                                                                                            cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGGCTTGCTATTAGTTGCGGTACAGCCCCTTTATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGGCGGCGCCAGTACCGGAc                                                         >  2:179267/1‑150 (MQ=39)
                                                                                              cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                                       >  1:919594/1‑150 (MQ=255)
                                                                                                cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                     >  1:370162/1‑150 (MQ=255)
                                                                                                cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCGTAATGATGTGATATCCGAACTGGGTGTGCGGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                     >  1:241992/1‑150 (MQ=255)
                                                                                                   cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                  >  1:402517/1‑150 (MQ=255)
                                                                                                     tctcccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                                >  1:611934/4‑150 (MQ=255)
                                                                                                      gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                               >  1:1628574/1‑150 (MQ=255)
                                                                                                      gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                               >  1:28085/1‑150 (MQ=255)
                                                                                                      caCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                               <  1:1713969/149‑1 (MQ=255)
                                                                                                         ccTCTCCCACAGGGAGAGGGAGAATACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTTCAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAGCCAc                                            >  2:406781/1‑150 (MQ=255)
                                                                                                            cacccACAGGGAGAGGGAGTAAACACTCACGGCATTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGCACAGCACCTAAATGATGGGAAATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACttt                                         <  1:1451868/148‑1 (MQ=255)
                                                                                                               ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAACACCACTTTATc                                      >  1:1178131/1‑150 (MQ=255)
                                                                                                                cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGTAGGCCTTGCTTTTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg                                     <  2:598637/150‑1 (MQ=255)
                                                                                                                cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg                                     <  2:126331/150‑1 (MQ=255)
                                                                                                                   agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAc                                  <  1:132469/150‑1 (MQ=255)
                                                                                                                      gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                               >  1:1391404/1‑150 (MQ=255)
                                                                                                                      gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                               >  2:536972/1‑150 (MQ=255)
                                                                                                                          gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGACCGCCGGaa                           >  1:1640570/1‑150 (MQ=255)
                                                                                                                          gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa                           >  1:1021656/1‑150 (MQ=255)
                                                                                                                           ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCGCCAGTACCGGACAACAGAAAACCACTTTATGGAACGCCGGAAc                          >  2:887215/1‑150 (MQ=255)
                                                                                                                            gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGCACGCCGGAAcc                         >  1:851197/1‑150 (MQ=255)
                                                                                                                               aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATc                      >  1:1875027/1‑150 (MQ=255)
                                                                                                                                  acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGACGTACAGCACCTAAAATATTAAAAAACCTAACTGGGTGTGAAGCATGCAGGTCGGCTCCAGTGCCGGCGTAGAGAAAAACACTTTATCGCTCGCCGGACGCATCTGa                   >  2:224864/1‑150 (MQ=255)
                                                                                                                                   cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAc                  <  2:1315748/150‑1 (MQ=255)
                                                                                                                                     cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAACACCACTTTATCGAACGCCGGACCCATCTGAccc                >  2:1639071/1‑150 (MQ=255)
                                                                                                                                     cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                >  2:1657405/1‑150 (MQ=255)
                                                                                                                                     cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                <  2:652529/150‑1 (MQ=255)
                                                                                                                                      tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt               <  2:191947/150‑1 (MQ=255)
                                                                                                                                      gcAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt               >  2:157749/2‑150 (MQ=255)
                                                                                                                                       gaaGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCGAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTg              <  2:617464/149‑1 (MQ=255)
                                                                                                                                         aGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACCGCCCCGTAATGATGTGATATCCGAACTGGGTGTGCAGCGGTCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACGTTAGCGAACGCCGGAACACTCTGCCCCTGGc            >  2:1163443/1‑150 (MQ=255)
                                                                                                                                         aGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTTCCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGc            >  2:86440/1‑150 (MQ=255)
                                                                                                                                         aGGCCTTCTCCTGGAGAAGGACTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGc            <  2:851962/150‑1 (MQ=255)
                                                                                                                                          ggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCg           >  2:475593/1‑150 (MQ=255)
                                                                                                                                           gCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGTCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCgg          <  1:93493/150‑1 (MQ=255)
                                                                                                                                              ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGATAAAACCACTTTTTCGAACGCCGGAACCATCTGACCCTGTCGGAAt       <  1:50458/150‑1 (MQ=255)
                                                                                                                                              ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAAt       <  1:1200427/150‑1 (MQ=255)
                                                                                                                                                caccTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTc     <  2:733471/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACC  >  NC_000913/3959367‑3959663

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