Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,502,106 A→C 5.1% A442A (GCT→GCG) ‡ insG ← KpLE2 phage‑like element; IS4 putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,502,1060AC 5.1% 560.2 / 15.7 172A442A (GCT→GCG) ‡insGKpLE2 phage‑like element; IS4 putative transposase
Reads supporting (aligned to +/- strand):  ref base A (75/88);  new base C (6/3);  total (81/91)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 174 total aligned reads displayed.

GGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCG  >  NC_000913/4501957‑4502254
                                                                                                                                                     |                                                                                                                                                    
ggCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTaa                                                                                                                                                      >  2:714020/1‑150 (MQ=255)
  cATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGc                                                                                                                                                    >  1:523225/1‑150 (MQ=255)
    tGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCaa                                                                                                                                                  >  1:2050/1‑150 (MQ=255)
     gACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAac                                                                                                                                                 >  2:385279/1‑150 (MQ=255)
     gACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAac                                                                                                                                                 <  1:660443/150‑1 (MQ=255)
     gACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAac                                                                                                                                                 <  1:691895/150‑1 (MQ=255)
      aCAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAact                                                                                                                                                >  2:880367/1‑150 (MQ=255)
      aCAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAact                                                                                                                                                >  1:1263611/1‑150 (MQ=255)
          tGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactgact                                                                                                                                            >  1:158784/1‑150 (MQ=255)
           gTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactgactg                                                                                                                                           >  1:31588/1‑150 (MQ=255)
            tCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTGCTTAACTGACTgg                                                                                                                                          <  1:724355/150‑1 (MQ=255)
              aaaCAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGct                                                                                                                                        >  1:1138700/1‑150 (MQ=255)
                   gTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCtttt                                                                                                                                   >  2:244324/1‑150 (MQ=255)
                   gTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCtttt                                                                                                                                   >  1:710635/1‑150 (MQ=255)
                   gTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCtttt                                                                                                                                   >  1:1420579/1‑150 (MQ=255)
                     cAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTc                                                                                                                                 >  2:843185/1‑150 (MQ=255)
                         gATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCgggg                                                                                                                             <  2:1091780/150‑1 (MQ=255)
                         gATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTAAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCgggg                                                                                                                             <  2:610508/150‑1 (MQ=255)
                          aTTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGc                                                                                                                            <  1:1325972/150‑1 (MQ=255)
                            tgtgGCGTTGACGGTCAGTCTTGGGGCAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCtgt                                                                                                                         <  2:1290949/150‑1 (MQ=255)
                                    tGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTAt                                                                                                                  <  2:1268135/150‑1 (MQ=255)
                                         gTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCa                                                                                                             >  1:111811/1‑150 (MQ=255)
                                          tCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAg                                                                                                            <  2:296407/150‑1 (MQ=255)
                                           cAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAgg                                                                                                           >  2:1239866/1‑150 (MQ=255)
                                           cAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAgg                                                                                                           <  2:311363/150‑1 (MQ=255)
                                            aGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAggg                                                                                                          <  1:144990/150‑1 (MQ=255)
                                            aGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAggg                                                                                                          <  2:1263611/150‑1 (MQ=255)
                                              tCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGcc                                                                                                        <  1:880367/150‑1 (MQ=255)
                                               cTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCct                                                                                                       >  2:717638/1‑150 (MQ=255)
                                                   gggACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTcc                                                                                                   <  1:461373/150‑1 (MQ=255)
                                                       cAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                               >  2:32188/1‑150 (MQ=255)
                                                       cAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                               >  1:1389437/1‑150 (MQ=255)
                                                         aaaaaGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTAcc                                                                                             >  1:1001321/1‑150 (MQ=255)
                                                         aaaaaGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTAcc                                                                                             <  1:32188/150‑1 (MQ=255)
                                                             aGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCActc                                                                                         >  1:951249/1‑150 (MQ=255)
                                                                 cAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCgg                                                                                     >  2:284078/1‑150 (MQ=255)
                                                                    tATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGaa                                                                                  <  2:158784/150‑1 (MQ=255)
                                                                      ttCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAgg                                                                                >  2:1282276/1‑150 (MQ=255)
                                                                      ttCTCCTTGGCAACCCGATAATTTAACACTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAgg                                                                                <  2:1036891/150‑1 (MQ=255)
                                                                       tctcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGc                                                                               <  2:768024/150‑1 (MQ=255)
                                                                         tcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGccc                                                                             <  2:195423/150‑1 (MQ=255)
                                                                         tcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGccc                                                                             <  1:843185/150‑1 (MQ=255)
                                                                          cGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCt                                                                            >  2:1325389/1‑150 (MQ=255)
                                                                             tGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTc                                                                         <  2:1138700/150‑1 (MQ=255)
                                                                                cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTggg                                                                                                                      <  1:1007796/102‑1 (MQ=255)
                                                                                cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGTGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCct                                                                      <  2:826338/150‑1 (MQ=255)
                                                                                cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCct                                                                      <  2:1001321/150‑1 (MQ=255)
                                                                                cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCct                                                                      <  2:448265/150‑1 (MQ=255)
                                                                                cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCct                                                                      <  1:1325389/150‑1 (MQ=255)
                                                                                  aCCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttc                                                                    >  1:615647/1‑150 (MQ=255)
                                                                                   ccAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttct                                                                   <  2:2050/150‑1 (MQ=255)
                                                                                   ccAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttct                                                                   <  2:405359/150‑1 (MQ=255)
                                                                                   ccAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttct                                                                   <  2:743291/150‑1 (MQ=255)
                                                                                     aGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTg                                                                 <  2:544524/150‑1 (MQ=255)
                                                                                       aTAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGTGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTc                                                               <  2:1121537/150‑1 (MQ=255)
                                                                                       aTAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTc                                                               <  1:894566/150‑1 (MQ=255)
                                                                                        tAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCg                                                              <  1:244324/150‑1 (MQ=255)
                                                                                              aaCAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAAtt                                                        <  2:1163906/150‑1 (MQ=255)
                                                                                              aaCAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAAtt                                                        >  2:690981/1‑150 (MQ=255)
                                                                                               aCAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAAttt                                                       <  2:338439/150‑1 (MQ=255)
                                                                                                cAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTc                                                      >  2:617734/1‑150 (MQ=255)
                                                                                                  gTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTcac                                                    <  1:682746/150‑1 (MQ=255)
                                                                                                   ttGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTcaca                                                   <  2:502213/150‑1 (MQ=255)
                                                                                                   ttGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTcaca                                                   <  2:1363331/150‑1 (MQ=255)
                                                                                                    tGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACaa                                                  <  2:523225/150‑1 (MQ=255)
                                                                                                    tGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACaa                                                  <  2:171167/150‑1 (MQ=255)
                                                                                                     gATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAg                                                 >  2:1020437/1‑150 (MQ=255)
                                                                                                      aTGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGt                                                <  1:461828/150‑1 (MQ=255)
                                                                                                           gAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTcc                                           >  1:777482/1‑150 (MQ=255)
                                                                                                           gAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTcc                                           <  2:924058/150‑1 (MQ=255)
                                                                                                           gAAACGAAATAGGCATCTTGGTTCTCTAATGCCACTCAGTTAAGCAACTCACTGGCTCTTTTTCGGGGCTATGGGGTATTTAAAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTATCCCTTCTTGTCGGTAATTTCACAATTAGTcc                                           >  1:99868/1‑150 (MQ=255)
                                                                                                               cGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATa                                       >  2:883828/1‑150 (MQ=255)
                                                                                                                   aTAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTg                                   <  1:558672/150‑1 (MQ=255)
                                                                                                                   aTAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTg                                   <  2:554697/150‑1 (MQ=255)
                                                                                                                     aGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCa                                 >  1:1185901/1‑150 (MQ=255)
                                                                                                                     aGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCa                                 <  2:343318/150‑1 (MQ=255)
                                                                                                                     aGGCATCTTGGTTCTGTAATGCCAGTCAGTAAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCa                                 >  2:677103/1‑150 (MQ=255)
                                                                                                                      ggCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCaa                                >  2:793/1‑150 (MQ=255)
                                                                                                                       gCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGGATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAg                               <  2:357140/150‑1 (MQ=255)
                                                                                                                          tCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGCCTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATc                            >  2:638290/1‑150 (MQ=255)
                                                                                                                           cTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATcg                           <  1:377608/150‑1 (MQ=255)
                                                                                                                           cTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATcg                           <  1:252216/150‑1 (MQ=255)
                                                                                                                              ggTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCgcgc                        <  1:116338/150‑1 (MQ=255)
                                                                                                                              ggTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCgcgc                        <  1:1409764/150‑1 (MQ=255)
                                                                                                                              ggTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTAACACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCgcgc                        <  2:945807/150‑1 (MQ=255)
                                                                                                                                 tCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATc                     >  1:588369/1‑150 (MQ=255)
                                                                                                                                 tCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCATGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATc                     >  2:1408935/1‑150 (MQ=255)
                                                                                                                                    gTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGc                  <  1:352139/150‑1 (MQ=255)
                                                                                                                                      aaTGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTc                >  1:1205145/1‑150 (MQ=255)
                                                                                                                                         gCCAGTCAGTTGCTTAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCgg             >  1:1253381/1‑150 (MQ=255)
                                                                                                                                         gCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGt            >  2:514821/1‑150 (MQ=255)
                                                                                                                                          ccAGTCAGTTGCTTAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGt            >  1:726602/1‑150 (MQ=255)
                                                                                                                                           cagtcagtTGCTTAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTATCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGta           >  1:432106/1‑150 (MQ=255)
                                                                                                                                           cagtcagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCACATACTTGCAAGATCGCGCATCAGCTCCGGta           >  1:506839/1‑150 (MQ=255)
                                                                                                                                            agtcagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCt                                                                            >  1:715400/1‑84 (MQ=255)
                                                                                                                                            agtcagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAAt                                                         >  2:106464/1‑103 (MQ=255)
                                                                                                                                            agtcagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAAt                                                         <  1:106464/103‑1 (MQ=255)
                                                                                                                                            agtcagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGtat          >  1:883774/1‑150 (MQ=255)
                                                                                                                                                 gtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTc     >  2:1377536/1‑150 (MQ=255)
                                                                                                                                                   tAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTccc   >  1:1093557/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCG  >  NC_000913/4501957‑4502254

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