Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,515 T→A 6.5% 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 6.5% 213.6 / 3.0 86intergenic (+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 (32/44);  new base A (3/3);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTC  >  NC_000913/3959366‑3959660
                                                                                                                                                     |                                                                                                                                                 
gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCAAACCCCGACCCACTCCCACAGGGAGATGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                                   >  1:1223385/1‑150 (MQ=255)
 cAGCCATGCGATTGAGCAGGTTGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                  <  1:392744/150‑1 (MQ=255)
 cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                  <  1:913962/150‑1 (MQ=255)
     cATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTACTGgaga                                                                                                                                              >  2:111522/1‑150 (MQ=255)
       ttcGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                            <  2:725030/148‑1 (MQ=255)
       tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                            <  1:389680/150‑1 (MQ=255)
          taTTGAGCAGGTAGGGAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                                         <  1:1105511/149‑1 (MQ=255)
          gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAACCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                                         <  2:822717/150‑1 (MQ=255)
           aTTGCGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                        <  2:871960/150‑1 (MQ=255)
           aTTGAGCAGGTAGGTAAGACACTGCGCGGCTATATGACAGATATGAAACGTATGGCTGTTGCGGCTGAAGTGCGCGCGGATGCACTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCACCAGGAGAAGGCCt                                                                                                                                        <  1:623308/150‑1 (MQ=38)
           aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                        <  2:768128/150‑1 (MQ=255)
           aTTGAGCAGGTAGGGAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                                        <  1:29232/150‑1 (MQ=255)
              gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                                                     <  2:151584/150‑1 (MQ=255)
               aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                                                    <  1:1164674/150‑1 (MQ=255)
               aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                                                    <  1:1194387/150‑1 (MQ=255)
                   ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCGTCAGCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCGGGAGAAGGCGGTGCTTTTa                                                                                                                                >  2:698593/1‑150 (MQ=38)
                   ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCAGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCGAACACACAAGGAAATGGATAAAAAACAAAAACCCATAACCAGGAGAAGGCCtttgaatta                                                                                                                                >  2:999995/1‑143 (MQ=255)
                       ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGTGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCACAGGGTGATGGAGAAAACCCTCAAGGCCGTAACCGGGAGACCGCCTTGCTATTATTTg                                                                                                                            >  1:549966/1‑150 (MQ=255)
                       ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                            >  1:874997/1‑150 (MQ=255)
                        gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTTAGTGCGCTCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAACCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                           <  1:78916/150‑1 (MQ=255)
                         taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCACTACCACAGCGACTGGGGGAAAAAACAAAAGGCAATCTCGTGGAGATGGCCGTCGTATTAgatcgg                                                                                                                          >  2:1115967/1‑145 (MQ=255)
                         taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                          <  1:568165/150‑1 (MQ=255)
                           aGAAACTGCGCGTCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                        <  2:326178/150‑1 (MQ=255)
                           aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCAGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                        <  1:897466/150‑1 (MQ=255)
                               aCTTCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                                                    <  2:881143/150‑1 (MQ=255)
                               aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                                                    <  2:1092316/150‑1 (MQ=255)
                                     cgtcTAGATGACACATATCAAACGTATTGCTGGTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGCCCCTCACCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                                              <  2:248089/147‑1 (MQ=38)
                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                      <  1:793130/150‑1 (MQ=255)
                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                      >  1:143184/1‑150 (MQ=255)
                                                  gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTAACCCCGACCCACACCCACAAGGTGGGGGTTGAAACACACAAAGCATTAAACTGGAGAAGGGCGTTCTATTATTTTCGGGAAATCACCTTtattattgtttat                                                                                                 >  2:1149660/1‑143 (MQ=255)
                                                   atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                                >  1:691305/1‑150 (MQ=255)
                                                   atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                                <  2:691305/150‑1 (MQ=255)
                                                    tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                               <  2:843657/150‑1 (MQ=255)
                                                    tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                               >  1:151767/1‑150 (MQ=255)
                                                      tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAAAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa                                                                                             <  2:742479/150‑1 (MQ=255)
                                                       gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                                                                                            >  1:20554/1‑150 (MQ=255)
                                                        aaaCGTATTGCTGTGGCGGGTTAAGTGCGCGCTGATGCCCTCTCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                                           <  1:307973/150‑1 (MQ=255)
                                                        aaaCGAATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCACCGACCCTCTACACCAGGGAGAGGGAGAAAACACTAAAGGCCGTCTCAGGGAGAAGGCCATGCTATTAGTTGCGTTACAGGACATTAAAGCTGAGATATGCGAAc                                                                                           >  1:1250022/1‑150 (MQ=255)
                                                           cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCACCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                        >  1:1234402/1‑150 (MQ=255)
                                                              aTTGCTCTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTTCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                                                                                     <  2:346140/150‑1 (MQ=255)
                                                               ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCAGAGGGTGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt                                                                                    <  2:821211/150‑1 (MQ=255)
                                                                 gttGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGTGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGGCGTGCTATTTGTTGCGGTACAGCCCCTGAATGATGTGATATCCGAACTGGGTGTGc                                                                                  >  2:619186/3‑150 (MQ=38)
                                                                      tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                                                                             <  2:603613/150‑1 (MQ=255)
                                                                          ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                                                         <  2:254580/150‑1 (MQ=255)
                                                                             tAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                                      <  1:1285923/150‑1 (MQ=255)
                                                                              aaGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCCCAGGGTGAGGGAGAAAACACACAAGGCCTTCTCCTGGAGAATGCCTTGCTATTATTAGGGGTTAAGGACCGTTATGATGTGATTTGCGAACTGGGTGTTCATCCGGCCGGTCg                                                                     >  2:644577/1‑150 (MQ=255)
                                                                               aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGTGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGGCTTGCTATTTGTTGCGGTACAGCACCGTAATGCTGTGATATCCGAACTGGGTGTGCAGCGGGCCGGGCgg                                                                    >  2:295932/1‑150 (MQ=38)
                                                                                   ggcgcTGATGCCCTCACCCCGACCGTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTTCTGGAGAAGGCGTTGCTATTAGTTGCGGTACAGCACCTTAATTATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                                                <  2:1100667/149‑1 (MQ=255)
                                                                                   cgcgcTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGGGATAACACACAAGGCCGTCGGCGGGGGCAGGCCTTGCGATTAGTGGCGGTACAGCCGCGTAGTGAGGGGATATCCGTACTGGGGGTGGAGCGGGCCGGGCGGCTcc                                                                >  2:246793/1‑150 (MQ=19)
                                                                                        ggATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                                                           <  1:1156910/149‑1 (MQ=255)
                                                                                         gATGCCCTCCCCCCGACCCTCGCCCACAGGGAGAGGGAGATAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTGAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGGAcc                                                          >  2:164870/1‑150 (MQ=255)
                                                                                             cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAc                                                      <  2:143184/150‑1 (MQ=255)
                                                                                              ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                                     >  1:928828/1‑150 (MQ=255)
                                                                                               cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGCAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                                    >  1:20440/1‑150 (MQ=255)
                                                                                                 cACCCCGACCCTCTCTCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                  >  2:933459/1‑150 (MQ=255)
                                                                                                 cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTCGAGAAGGCCTTGCTATTAGTTGCGGTTCAGCCCCTTAATGATGTGCTATCCGAACTGGGTGTGCAGCGGGCCGGTGGGCGCCAGTTCCGGACAtga                                                  >  1:1202833/1‑150 (MQ=255)
                                                                                                  aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag                                                 <  1:51469/150‑1 (MQ=255)
                                                                                                  aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag                                                 <  2:879697/150‑1 (MQ=255)
                                                                                                   ccccGCCCCTCTCCCACAGGGAGATGGAGAAAACACTCAAGGCCTACTCCTGGAGAAGGCCTGGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGGGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                <  2:892510/150‑1 (MQ=255)
                                                                                                     ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTTCCGGACAAGAGaaa                                              >  1:801272/1‑150 (MQ=255)
                                                                                                       gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTTATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCGCCAGTACCGGACAAGCGaaacc                                            >  2:236379/1‑148 (MQ=255)
                                                                                                       gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                            <  1:479986/150‑1 (MQ=255)
                                                                                                          ccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTAAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCGAGTACCGGACAAGAGAAAACCAc                                         >  1:873605/1‑150 (MQ=255)
                                                                                                           ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATGCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt                                        >  1:222568/1‑150 (MQ=255)
                                                                                                            tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGCAAACCACtt                                       >  2:1011088/1‑150 (MQ=255)
                                                                                                              tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGAAGCGGGCCGGTCGGCTCCAGTACCTGACAAGAGAAAACCACttaa                                     >  1:1285059/1‑148 (MQ=255)
                                                                                                               cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGATGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGCACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAACCCACTTTAt                                    >  2:554766/1‑150 (MQ=255)
                                                                                                                ccACAGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTAGTAGTTGCGGTACAGCCCCTTAATGATGTGATCTCCGCACTGGGTGGGCAGCGGGCCGGTCGGCGCCAGTCCCGGACAAGAGAAAACCCCTTTATc                                   >  1:1186712/1‑150 (MQ=255)
                                                                                                                 cccaGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTGGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCGCCAGTACCGGACAAGAGAAAACCACTTTATCg                                  <  1:763089/148‑1 (MQ=255)
                                                                                                                  acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGa                                 <  1:582538/150‑1 (MQ=255)
                                                                                                                     gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACg                              <  2:874997/150‑1 (MQ=255)
                                                                                                                      ggatagtgagaAAACACTCAATGCCTTCTACTTGAGAAGGCCTTGCTATTAGTTGCGGTACAGAACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc                             <  2:1223385/150‑1 (MQ=255)
                                                                                                                       gatagtgagaAATCTCTCCAGTCCTTCTCCTGGCTAAGGCCTTGCTAATAGTTGCGGTAAAGCACCTAAATGATGTGTTATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                            <  2:189067/144‑1 (MQ=255)
                                                                                                                            ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc                       <  1:1284229/150‑1 (MQ=255)
                                                                                                                             gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTGCCTGACAAGAGCAAACCCCTTTATCGAACGCCGGAAcc                      >  2:89318/1‑150 (MQ=255)
                                                                                                                             gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAGGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGACGGCTCCAGTACCGGACAAGAGAAAACCACGTTATCGAACGCCGGCAcc                      >  1:653727/1‑150 (MQ=255)
                                                                                                                                  aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACATCACCTTAATGATGTGATATCCGAACTGGGTGTGAAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTg                 <  2:907050/150‑1 (MQ=255)
                                                                                                                                  aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGCACCATCTg                 >  1:758056/1‑150 (MQ=255)
                                                                                                                                   acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                <  2:779890/150‑1 (MQ=255)
                                                                                                                                      cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc             >  1:396819/1‑150 (MQ=255)
                                                                                                                                      agcAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGTCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc             >  1:237155/3‑150 (MQ=255)
                                                                                                                                        cAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGGGGGCCGGGCGGCTCGAGTACCGGAAAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTg           <  2:486950/150‑1 (MQ=255)
                                                                                                                                        cAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTg           >  1:757453/1‑150 (MQ=255)
                                                                                                                                              cTTCTGGTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCCCCTTAATGATGTGATAGCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACGGGACAAGTGAAAAGCAGTTTTTCGGACGCCGGAAGCATCGGACCCTGGCggcc     >  2:1042853/1‑148 (MQ=255)
                                                                                                                                               ttCTCCTGGAGAAGGCCTTGCTATTCGTTGCGGAACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAAt    <  1:253349/150‑1 (MQ=255)
                                                                                                                                                tctcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATACGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAAtt   >  2:516122/1‑150 (MQ=255)
                                                                                                                                                 ctcCTGGAGAAGACCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTGTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTc  >  2:799576/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTC  >  NC_000913/3959366‑3959660

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: