Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,197,226 C→G 6.0% V806L (GTG→CTG)  glnE ← fused deadenylyltransferase/adenylyltransferase for glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,197,2260CG 6.0% 284.9 / 8.1 84V806L (GTG→CTG) glnEfused deadenylyltransferase/adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base C (42/37);  new base G (3/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

TCGGTAATTCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGA  >  NC_000913/3197081‑3197374
                                                                                                                                                 |                                                                                                                                                    
tCGGTAATTCCGCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGc                                                                                                                                                   >  3:934423/1‑149 (MQ=255)
tCGGTAATTCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGc                                                                                                                                                   >  3:491839/1‑149 (MQ=255)
   gTAATTCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTc                                                                                                                                                >  3:885659/1‑149 (MQ=255)
      aTTCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCaaa                                                                                                                                             >  4:693292/1‑149 (MQ=255)
      aTTCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCaaa                                                                                                                                             <  4:105187/149‑1 (MQ=255)
      aTTCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCGTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGGCaaa                                                                                                                                             <  2:912610/149‑1 (MQ=255)
       ttCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCAGTGCGTcgcga                                                                                                                                            <  4:305767/149‑5 (MQ=255)
           cccTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTgcg                                                                                                                                        >  3:286826/1‑149 (MQ=255)
             cTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGcgcg                                                                                                                                      >  4:682782/1‑149 (MQ=255)
                  tcaGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTgcgc                                                                                                                                       <  4:721931/143‑1 (MQ=255)
                  tcaGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTgcgc                                                                                                                                       >  3:721931/1‑143 (MQ=255)
                  tcaGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAg                                                                                                                                 >  3:373128/1‑149 (MQ=255)
                  tcaGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAg                                                                                                                                 >  1:309796/1‑149 (MQ=255)
                  tcaGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAg                                                                                                                                 <  3:993305/149‑1 (MQ=255)
                     gCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTg                                                                                                                              >  3:397099/1‑149 (MQ=255)
                       tttGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCg                                                                                                                            >  2:1035308/1‑149 (MQ=255)
                       tttGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCg                                                                                                                            >  3:407136/1‑149 (MQ=255)
                        ttGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCgg                                                                                                                           >  2:1091035/1‑149 (MQ=255)
                           atatCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATc                                                                                                                        >  4:47630/1‑149 (MQ=255)
                           atatCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATc                                                                                                                        <  1:147345/149‑1 (MQ=255)
                              tCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGcc                                                                                                                     <  1:138194/149‑1 (MQ=255)
                              tCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGcc                                                                                                                     >  2:174075/1‑149 (MQ=255)
                              tCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGACAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGcc                                                                                                                     >  3:655140/1‑149 (MQ=255)
                                aaaGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGt                                                                                                                   >  4:462642/1‑149 (MQ=255)
                                aaaGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGt                                                                                                                   >  3:1060602/1‑149 (MQ=255)
                                aaaGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGt                                                                                                                   >  2:61609/1‑149 (MQ=255)
                                  aGCGATCGCGATGGTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGGTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGtac                                                                                                                 >  2:655420/1‑149 (MQ=255)
                                    cGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGtacac                                                                                                               <  2:887398/149‑1 (MQ=255)
                                      aTCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGtacacta                                                                                                             >  1:207840/1‑149 (MQ=255)
                                        cgcgATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGtacactaca                                                                                                           <  4:918817/149‑1 (MQ=255)
                                           gATGTTTATCGCGGAGATGAGCGTGCATTTTCTCTCGGATTCCGCGCACTTCCGTTTGCAGAGTTTTACCTGCACGCGGCAGCGACATAATCTCGCGGCGCACTGCGTCAAAGTGCGCGTTGAGCTACGGACGCCCGTGCACTACAcgc                                                                                                        <  2:38513/149‑1 (MQ=255)
                                                 tATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGcgcacgc                                                                                                  >  1:387535/1‑149 (MQ=255)
                                                     gcggagaTGAGCTCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGAGGCAGCTTCATAATCTCGTGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGTACGCACCa                                                                                              <  4:923222/146‑1 (MQ=255)
                                                     gCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCa                                                                                              >  2:121722/1‑149 (MQ=255)
                                                       cGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAgc                                                                                            <  3:918381/149‑1 (MQ=255)
                                                          gaTGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGcc                                                                                         <  4:89519/149‑1 (MQ=255)
                                                           aTGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCt                                                                                        >  4:208962/1‑149 (MQ=255)
                                                                         tctcGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTc                                                                          <  4:1037847/149‑1 (MQ=255)
                                                                           tcGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGTGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCa                                                                        >  2:746819/1‑149 (MQ=255)
                                                                             gcgcATTTCCCGCACTTCCGTTTGCAGAGCTTTACCTTCACGCGACAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAgg                                                                      >  1:624101/1‑149 (MQ=255)
                                                                                cATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGa                                                                                                                          <  3:511691/94‑1 (MQ=255)
                                                                                cATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGa                                                                                                                          >  4:511691/1‑94 (MQ=255)
                                                                                cATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCt                                                                   <  1:1035308/149‑1 (MQ=255)
                                                                                   ttCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGTGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGt                                                                <  1:746819/149‑1 (MQ=255)
                                                                                    tCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGtt                                                               <  1:649587/149‑1 (MQ=255)
                                                                                       cGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGttttt                                                            <  3:838314/149‑1 (MQ=255)
                                                                                        gCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTc                                                           <  3:611630/149‑1 (MQ=255)
                                                                                        gCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTc                                                           <  3:855064/149‑1 (MQ=255)
                                                                                         cACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCt                                                          <  4:529452/149‑1 (MQ=255)
                                                                                           cTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGa                                                        <  2:91648/149‑1 (MQ=255)
                                                                                            ttCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGAt                                                       <  2:374068/149‑1 (MQ=255)
                                                                                            ttCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGAt                                                       >  1:748898/1‑149 (MQ=255)
                                                                                               cGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAAt                                                    <  1:618942/149‑1 (MQ=255)
                                                                                                    gCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCa                                               <  4:30462/149‑1 (MQ=255)
                                                                                                      agagTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCaaa                                             >  3:963614/1‑149 (MQ=255)
                                                                                                      agagTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCaaa                                             <  3:273012/149‑1 (MQ=255)
                                                                                                      agagTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCaaa                                             <  2:88209/149‑1 (MQ=255)
                                                                                                          ttttACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCt                                         >  4:29617/1‑149 (MQ=255)
                                                                                                             tACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCt                                      <  2:876650/149‑1 (MQ=255)
                                                                                                             tACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCt                                      <  1:441227/149‑1 (MQ=255)
                                                                                                               ccTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGc                                    <  4:590328/149‑1 (MQ=255)
                                                                                                               ccTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGc                                    <  1:633941/149‑1 (MQ=255)
                                                                                                                cTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCg                                   >  4:411023/1‑149 (MQ=255)
                                                                                                                cTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCg                                   >  1:822669/1‑149 (MQ=255)
                                                                                                                cTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCCCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGAGGTTTCCACGTCCAGGCCTCGCTTTTCTGATAATCGACAAATTCTTCTGCg                                   >  4:50092/1‑149 (MQ=255)
                                                                                                                 ttCACGCGGCAGCGTCATAATCTCGCGATGCACTGCGGCAAAGTGCGAGACGAGCTGCGGATCGCCGTACACTACACCCGGACGCACCAGCGCCTGGTGGTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTcggcgg                                  >  4:499942/1‑149 (MQ=255)
                                                                                                                 ttCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGa                                                        <  1:999125/127‑1 (MQ=255)
                                                                                                                 ttCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGa                                                        >  2:999125/1‑127 (MQ=255)
                                                                                                                 ttCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCgg                                  >  2:1056559/1‑149 (MQ=255)
                                                                                                                 ttCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCgg                                  >  1:958386/1‑149 (MQ=255)
                                                                                                                  tCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGa                                 <  3:176302/149‑1 (MQ=255)
                                                                                                                   cACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGAt                                <  1:660902/149‑1 (MQ=255)
                                                                                                                      gcgGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTc                             <  4:286826/149‑1 (MQ=255)
                                                                                                                      gcgGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTc                             >  4:658745/1‑149 (MQ=255)
                                                                                                                      gcgGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTc                             >  1:334678/1‑149 (MQ=255)
                                                                                                                       cgGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCa                            <  3:658745/149‑1 (MQ=255)
                                                                                                                           aGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAg                        >  3:448688/1‑149 (MQ=255)
                                                                                                                                   aaTCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGt                              >  2:236059/1‑135 (MQ=255)
                                                                                                                                   aaTCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGt                              <  1:236059/135‑1 (MQ=255)
                                                                                                                                     tctcGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCccgcc              >  2:185893/1‑149 (MQ=255)
                                                                                                                                        tttgACGCAGTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGAt                                >  3:1063259/4‑128 (MQ=255)
                                                                                                                                        tttgACGCAGTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGAt                                <  4:1063259/125‑1 (MQ=255)
                                                                                                                                        tttgACGCAGTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCgccgcc           >  2:272406/4‑149 (MQ=255)
                                                                                                                                        tttgACGCAGTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCgccgcc           >  1:391443/4‑149 (MQ=255)
                                                                                                                                                 cTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACGCGCTCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTGTGATAATCGGCAAGTTCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCcggacgga  >  4:518574/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCGGTAATTCCCCCTTCATCAGCTTTGATATCAAAGCGATCGCGATGTTTATTGCCGAGATGAGCGCGCATTTTCTCGCGCATTTCCCGCACTTCCGTTTGCAGAGTTTTACCTTCACGCGGCAGCGTCATAATCTCGCGACGCACTGCGTCAAAGTGCGCGGTGAGCTGCGGATCGCCGTACACTACACGCGCACGCACCAGCGCCTGATGTTCCCACGTCCAGGCCTCGTTTTTCTGATAATCGGCAAATGCTTCTGCGGATGTCACCAGCATTCCCGCCGCCCCGGACGGA  >  NC_000913/3197081‑3197374

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