Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,502,108 C→T 5.1% A442T (GCT→ACT) ‡ 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,1080CT 5.1% 367.6 / 5.8 118A442T (GCT→ACT) ‡insGKpLE2 phage‑like element; IS4 putative transposase
Reads supporting (aligned to +/- strand):  ref base C (57/53);  new base T (3/3);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

Only 100 of 121 total aligned reads displayed.

CATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTG  >  NC_000913/4501959‑4502257
                                                                                                                                                     |                                                                                                                                                     
cATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGc                                                                                                                                                       <  2:577061/150‑1 (MQ=255)
cATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACAAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGc                                                                                                                                                       <  2:1605551/150‑1 (MQ=255)
  tGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCaa                                                                                                                                                     >  1:558401/1‑150 (MQ=255)
   gACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTAATCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAac                                                                                                                                                    >  1:1162579/1‑150 (MQ=255)
         gTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactgactg                                                                                                                                              <  1:460803/150‑1 (MQ=255)
          tCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTGCTTAACTGACTgg                                                                                                                                             <  1:287216/150‑1 (MQ=255)
          tCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTGCTTAACTGACTgg                                                                                                                                             <  1:1243309/150‑1 (MQ=255)
          tCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTGCTTAACTGACTgg                                                                                                                                             <  1:666441/150‑1 (MQ=255)
           cAAACAGTCAAAGATTGTGGCGTTTACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGATCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                                                                                                                            <  1:1313230/150‑1 (MQ=255)
           cAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                                                                                                                            >  2:1498405/1‑150 (MQ=255)
           aaaaCAGTCAAAGATTGTGGCGTTGACGGGCAGTCTTGTGGACAAAATAGCTCAGTATTCTCGTTGGCAACCAGATTATTTAAGAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                                                                                                                            <  2:1377152/149‑1 (MQ=255)
             aaCACTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTATCGTTGGCAACCAGATAAATTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTTGcta                                                                                                                                          >  1:1396451/1‑149 (MQ=255)
              aCAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGctct                                                                                                                                         >  2:892121/1‑150 (MQ=255)
               cAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCtt                                                                                                                                        <  1:83938/150‑1 (MQ=255)
                 gTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCtttt                                                                                                                                      >  1:590251/1‑150 (MQ=255)
                   cAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                                                                                                                            >  1:68790/1‑142 (MQ=255)
                   cAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                                                                                                                            <  2:68790/142‑1 (MQ=255)
                   cAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTc                                                                                                                                    >  1:579806/1‑150 (MQ=255)
                   cAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTc                                                                                                                                    >  2:559178/1‑150 (MQ=255)
                    aaaGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCg                                                                                                                                   >  1:1498763/1‑150 (MQ=255)
                     aaGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCgg                                                                                                                                  >  1:813711/1‑150 (MQ=255)
                      aGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCggg                                                                                                                                 <  2:51947/150‑1 (MQ=255)
                        aTTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGc                                                                                                                               <  2:414756/150‑1 (MQ=255)
                        aTTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGc                                                                                                                               <  2:1046931/150‑1 (MQ=255)
                          tgtgGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCtg                                                                                                                             >  1:36589/1‑150 (MQ=255)
                            tgGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCtgtg                                                                                                                           >  2:922665/1‑150 (MQ=255)
                             ggCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTgg                                                                                                                          <  2:891916/150‑1 (MQ=255)
                             ggCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTgg                                                                                                                          >  1:755036/1‑150 (MQ=255)
                              gCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTggg                                                                                                                         <  1:549686/150‑1 (MQ=255)
                                gTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGt                                                                                                                       <  2:1349782/150‑1 (MQ=255)
                                  tGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTAt                                                                                                                     >  1:602340/1‑150 (MQ=255)
                                     cGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTc                                                                                                                  >  1:1409532/1‑150 (MQ=255)
                                     cGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCGGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTc                                                                                                                  >  1:153049/1‑150 (MQ=255)
                                      ggTCAGTCTTGGGGACAAAAAGGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTcc                                                                                                                 >  1:394788/1‑150 (MQ=255)
                                        tCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATCGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAg                                                                                                               <  1:433228/150‑1 (MQ=255)
                                         cAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAgg                                                                                                              >  2:436089/1‑150 (MQ=255)
                                           gTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGc                                                                                                            >  1:1157429/1‑150 (MQ=255)
                                            tCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCGGAAACGAAATAGGCATCTTGGGTCTGTAATGCCAGTCAGTTACGCAACTGACTGGCTCTTTTTCGGGGTTGTGGGGGATTTCCAGGGcc                                                                                                           <  2:1315356/150‑1 (MQ=255)
                                                 gggACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTcc                                                                                                      <  1:264106/150‑1 (MQ=255)
                                                  ggACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCCTCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCt                                                                                                     >  1:620024/1‑150 (MQ=255)
                                                   gACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCtt                                                                                                    >  1:296585/1‑150 (MQ=255)
                                                     cAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                                  <  2:1625511/150‑1 (MQ=255)
                                                     cAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                                  >  1:912778/1‑150 (MQ=255)
                                                     cAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATCTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                                  >  2:102734/1‑150 (MQ=255)
                                                           aGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCActc                                                                                            >  1:1606951/1‑150 (MQ=255)
                                                            gCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCActct                                                                                           >  1:503708/1‑150 (MQ=255)
                                                            gCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCActct                                                                                           >  2:672228/1‑150 (MQ=255)
                                                             cTCAGTATTCTCGTTGGCAAACAGATAATTTAACAGTTGATGATGAAACGAAATAGGTATATTTGTTCTGTAATGCCAGTAAGTTAAGAAACTGAATGGATATTTTTCGTGGCTGTGTGTTATTTCCAGTTCCTCTCCTTTAAAACtctc                                                                                          <  2:818038/150‑1 (MQ=255)
                                                               cAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGTGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCgg                                                                                        <  2:207903/150‑1 (MQ=255)
                                                                 gTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGa                                                                                      <  2:579806/150‑1 (MQ=255)
                                                                   aTTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAg                                                                                    <  1:672228/150‑1 (MQ=255)
                                                                    ttCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAgg                                                                                   <  1:1581187/150‑1 (MQ=255)
                                                                         gTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCtt                                                                              <  1:439494/150‑1 (MQ=255)
                                                                            ggCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTATACCACTCTCGGGAAGGCCCTTTcc                                                                           >  2:1174237/1‑150 (MQ=255)
                                                                            ggCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTAAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTcc                                                                           >  2:863149/1‑150 (MQ=255)
                                                                                aCCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATTACAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttc                                                                       >  1:12458/1‑150 (MQ=255)
                                                                                aCCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAGGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttc                                                                       >  2:1601027/1‑150 (MQ=255)
                                                                                 ccAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttct                                                                      >  1:744115/1‑150 (MQ=255)
                                                                                   aGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttcttt                                                                    >  2:137077/1‑149 (MQ=255)
                                                                                   aGATAATTTAACAGTTGATGCTGAAAAGAAATAGGGATGTTTGTTCTGTAATGCCAGTCTGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTTTGGGGTATTTCCATGGGGTGTAATTTACAACTCTCGGTAAGGGAGTTTCCCttgttg                                                                    <  2:259841/150‑1 (MQ=255)
                                                                                      tAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCg                                                                 >  1:595462/1‑150 (MQ=255)
                                                                                      tAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCg                                                                 <  2:602340/150‑1 (MQ=255)
                                                                                      tAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCg                                                                 >  1:1088105/1‑150 (MQ=255)
                                                                                      tAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCATTTACCACTCTCGGGAAGGCCCTTGCCGTTCTTGtgg                                                                 >  2:860077/1‑148 (MQ=255)
                                                                                       aaTTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCgg                                                                <  2:5158/150‑1 (MQ=255)
                                                                                               aGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCAAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTca                                                        >  2:1144224/1‑150 (MQ=255)
                                                                                                gTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTcac                                                       <  2:979289/150‑1 (MQ=255)
                                                                                                   gATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAg                                                    <  2:1409532/150‑1 (MQ=255)
                                                                                                    aTGGTGAATCGAAATAGGCATCTTGGTTCTGTTATGCCATTCAGTTAAGCAACTGACTGGCTATTTTTCGGTTCTGTGGTGTATTTGCAGGGCCTCTCCTTTTCCACTCTCTGGAAGGCCCTTTCGCTTCTTGTCTTTAATTTCACAAGt                                                   <  2:1157429/150‑1 (MQ=255)
                                                                                                      gCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTg                                                 <  1:514037/150‑1 (MQ=255)
                                                                                                       cTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCct                                                                         >  1:1244738/1‑125 (MQ=255)
                                                                                                       cTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCct                                                                         <  2:1244738/125‑1 (MQ=255)
                                                                                                       cTGAAACGAAATAGGCATCTTGGGTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGt                                                >  2:791452/1‑150 (MQ=255)
                                                                                                        tGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTc                                               >  2:74937/1‑150 (MQ=255)
                                                                                                          aaaCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTATGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTccc                                             <  2:152811/150‑1 (MQ=255)
                                                                                                          aaaCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTccc                                             <  2:848584/150‑1 (MQ=255)
                                                                                                          aaaCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTccc                                             <  1:436089/150‑1 (MQ=255)
                                                                                                           aaCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGATAATTTCACAAGTTGTCCCa                                            <  1:102734/150‑1 (MQ=255)
                                                                                                                 aTAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTg                                      <  2:361734/150‑1 (MQ=255)
                                                                                                                   aGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCa                                    <  2:590251/150‑1 (MQ=255)
                                                                                                                     gCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCaac                                  >  1:1309467/1‑149 (MQ=255)
                                                                                                                            ggTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCgcgc                           <  1:853797/150‑1 (MQ=255)
                                                                                                                             gTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCa                          <  1:49932/150‑1 (MQ=255)
                                                                                                                                cTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCa                       >  2:1446980/1‑150 (MQ=255)
                                                                                                                                     aTGCCAGTCAGTTAATCAACTGACTGGATCTTTTTCGGGGCTGTGGGCTATTTACAGGGCCTCTCATTTACCACTCTCGGGAATGCCAATTCCCATCTAGTCCGTAATTACACAAGTAGTCCCCTACTTGCAAGATCGCGCATCAGCTcc                  >  2:1027080/1‑150 (MQ=255)
                                                                                                                                      tGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCg                 >  2:638938/1‑150 (MQ=255)
                                                                                                                                       gCCAGTCAGTTGCTTAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCgg                >  2:537664/1‑150 (MQ=255)
                                                                                                                                       gCCAGTCAGTTGCTTAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCgg                >  2:451984/1‑150 (MQ=255)
                                                                                                                                       gCCAGTCAGTTGCTTAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCACATACTTGCAATATCGCGCAACAGCTCCgg                >  2:843894/1‑150 (MQ=255)
                                                                                                                                       gCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCgg                <  1:645676/150‑1 (MQ=255)
                                                                                                                                       gCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCgg                >  1:1104405/1‑150 (MQ=255)
                                                                                                                                        ccAGTCAGTTAATCAACTGACTGGCTCTTTTTAGGGGATGTGGGGTATTTCCAGGGCCTCTAATTTAACAATCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGt               >  2:1013524/1‑150 (MQ=255)
                                                                                                                                        ccAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGt               >  1:1430472/1‑150 (MQ=255)
                                                                                                                                          agtcagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGtat             <  2:533266/150‑1 (MQ=255)
                                                                                                                                             cagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACg          >  1:1450331/1‑150 (MQ=255)
                                                                                                                                             cagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACg          >  2:1134182/1‑150 (MQ=255)
                                                                                                                                                 tAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTccc      >  1:1405247/1‑150 (MQ=255)
                                                                                                                                                   aGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGTGCCTCTCCTTTACCACTCTCGGGAAGTCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCgg    <  2:467880/150‑1 (MQ=255)
                                                                                                                                                   aGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCgg    <  2:1498763/150‑1 (MQ=255)
                                                                                                                                                     cAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTATCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTACCGGTg  >  1:1512343/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTG  >  NC_000913/4501959‑4502257

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