Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,398 T→C 7.2% intergenic (+22/+21) corA → / ← yigF Ni(2(+))/Co(2(+))/Mg(2(+)) transporter/DUF2628 domain‑containing protein YigF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,3980TC 7.2% 329.7 / 11.9 111intergenic (+22/+21)corA/yigFNi(2(+))/Co(2(+))/Mg(2(+)) transporter/DUF2628 domain‑containing protein YigF
Reads supporting (aligned to +/- strand):  ref base T (53/50);  new base C (4/4);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 112 total aligned reads displayed.

GCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAG  >  NC_000913/4002249‑4002546
                                                                                                                                                     |                                                                                                                                                    
gttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGt                                                                                                                                                      >  1:1494705/5‑150 (MQ=255)
   ctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTg                                                                                                                                                   >  1:707915/2‑150 (MQ=255)
   ctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTg                                                                                                                                                   >  1:1361846/2‑150 (MQ=255)
   ctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGCGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTg                                                                                                                                                   >  2:1036252/2‑150 (MQ=255)
   aTGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTg                                                                                                                                                   <  1:334215/150‑1 (MQ=255)
    tGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTgg                                                                                                                                                  >  1:725687/1‑150 (MQ=255)
    tGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTgg                                                                                                                                                  <  2:1361846/150‑1 (MQ=255)
      gCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCt                                                                                                                                                >  2:422110/1‑150 (MQ=255)
             cTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAc                                                                                                                                         <  1:974361/150‑1 (MQ=255)
             cTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAc                                                                                                                                         <  1:310752/150‑1 (MQ=255)
              tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAcc                                                                                                                                        <  1:733192/150‑1 (MQ=255)
               ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCt                                                                                                                                       >  1:928449/1‑150 (MQ=255)
               ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCt                                                                                                                                       >  2:649852/1‑150 (MQ=255)
                tGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccgc                                                                                                                                      <  1:847214/150‑3 (MQ=255)
                  aGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGtt                                                                                                                                    >  2:1303978/1‑150 (MQ=255)
                  aGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGtt                                                                                                                                    >  1:160016/1‑150 (MQ=255)
                   gTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa                                                                                                                                   >  2:671401/1‑150 (MQ=255)
                   gTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccgctct                                                                                                                                   <  1:103022/150‑6 (MQ=255)
                       aTGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCa                                                                                                                               <  2:425622/150‑1 (MQ=255)
                        tGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCaa                                                                                                                              >  1:1343780/1‑150 (MQ=255)
                         gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCaaa                                                                                                                             <  2:1293203/150‑1 (MQ=255)
                          ccAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAg                                                                                                                            >  2:60439/1‑150 (MQ=255)
                          ccAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAg                                                                                                                            >  1:989077/1‑150 (MQ=255)
                          ccAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAg                                                                                                                            >  1:728177/1‑150 (MQ=255)
                            aGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGc                                                                                                                          <  1:491016/150‑1 (MQ=255)
                               aCTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTcc                                                                                                                       >  1:652657/1‑150 (MQ=255)
                                cTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCa                                                                                                                      >  1:1263356/1‑150 (MQ=255)
                                  gAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAgg                                                                                                                    >  1:1296355/1‑150 (MQ=255)
                                    aGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTa                                                                                                                  <  2:567439/150‑1 (MQ=255)
                                     gTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTAt                                                                                                                 <  2:725687/150‑1 (MQ=255)
                                            ttCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAccc                                                                                                          <  1:1391845/150‑1 (MQ=255)
                                            ttCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAccc                                                                                                          <  1:553989/150‑1 (MQ=255)
                                                gCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTCACCCtttt                                                                                                      <  2:886608/150‑1 (MQ=255)
                                                    cccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACct                                                                                                  >  2:607570/1‑150 (MQ=255)
                                                     ccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctc                                                                                                 >  2:589679/1‑150 (MQ=255)
                                                     ccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctc                                                                                                 <  2:1436519/150‑1 (MQ=255)
                                                     ccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctc                                                                                                 >  1:1135636/1‑150 (MQ=255)
                                                      cTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctct                                                                                                >  2:1534840/1‑150 (MQ=255)
                                                      cTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctct                                                                                                >  1:954099/1‑150 (MQ=255)
                                                       gggCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCtt                                                                                               <  2:204412/149‑1 (MQ=255)
                                                          cgcgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTtct                                                                                            <  2:492306/150‑1 (MQ=255)
                                                            cgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCa                                                                                          <  1:973010/150‑1 (MQ=255)
                                                             gATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCAt                                                                                         <  1:1303978/150‑1 (MQ=255)
                                                              aTTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATa                                                                                        <  1:371879/150‑1 (MQ=255)
                                                                tATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGa                                                                                      <  1:1270227/150‑1 (MQ=255)
                                                                   cTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAAcc                                                                                   >  2:1309833/1‑150 (MQ=255)
                                                                   cTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAAcc                                                                                   <  1:1010949/150‑1 (MQ=255)
                                                                    tttATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCa                                                                                  >  2:1293246/1‑150 (MQ=255)
                                                                       aTGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAttt                                                                               >  1:111303/1‑150 (MQ=255)
                                                                         gATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGt                                                                             >  2:30412/1‑150 (MQ=255)
                                                                          aTCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGtt                                                                            >  1:517842/1‑150 (MQ=255)
                                                                          aTCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGtt                                                                            <  2:652657/150‑1 (MQ=255)
                                                                           tCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTc                                                                           >  1:769287/1‑150 (MQ=255)
                                                                           tCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTc                                                                           <  1:1036252/150‑1 (MQ=255)
                                                                            ccTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCg                                                                          >  1:902528/1‑150 (MQ=255)
                                                                             cTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCgt                                                                         >  1:1490878/1‑150 (MQ=255)
                                                                                gcgGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCATGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGtt                                                                      <  2:928449/150‑1 (MQ=255)
                                                                                  gggCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTaa                                                                    <  1:649852/150‑1 (MQ=255)
                                                                                       tGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCATTTTAACTCTTTCTCATAGAACCATATGTTCGTGTTAACAGCa                                                               >  2:840095/1‑150 (MQ=255)
                                                                                         gCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAAc                                                             >  1:1029627/1‑150 (MQ=255)
                                                                                           aCCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACAt                                                           >  1:741410/1‑150 (MQ=255)
                                                                                               tATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgc                                                       <  2:1494705/150‑1 (MQ=255)
                                                                                                  cTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgc                                                    >  1:1428547/1‑150 (MQ=255)
                                                                                                       cTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACgg                                               >  1:475672/1‑150 (MQ=255)
                                                                                                       cTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCGTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTATTCGTGTTAACAGCAACATAGGCTGCTACgg                                               <  2:858717/150‑1 (MQ=255)
                                                                                                        tttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                              <  2:1135636/150‑1 (MQ=255)
                                                                                                        tttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                              <  2:1235761/150‑1 (MQ=255)
                                                                                                        tttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                              <  1:607570/150‑1 (MQ=255)
                                                                                                          tAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCaa                                            >  2:907743/1‑150 (MQ=255)
                                                                                                             gCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTc                                         <  1:1034023/150‑1 (MQ=255)
                                                                                                               ggAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTccc                                       >  1:246375/1‑150 (MQ=255)
                                                                                                                 aagaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAg                                     <  2:317599/150‑1 (MQ=255)
                                                                                                                  agaaCTGGTTGTAAAAAACGAGAGCGGTGTCTTAGTCTGGCTAAGCCACATGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAAACATTAGTTCGTGTTAACAGCAACAAAGGCTGCTAAGGCAATTCCCAGa                                    >  2:428200/1‑150 (MQ=255)
                                                                                                                  agaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGa                                    <  2:707915/150‑1 (MQ=255)
                                                                                                                     aCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATg                                 >  1:1039715/1‑150 (MQ=255)
                                                                                                                        ggTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTa                              <  1:633415/150‑1 (MQ=255)
                                                                                                                        ggTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTa                              <  2:246375/150‑1 (MQ=255)
                                                                                                                         gTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTaa                             >  1:1206934/1‑150 (MQ=255)
                                                                                                                         gTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTaa                             >  1:773108/1‑150 (MQ=255)
                                                                                                                          ttGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAAc                            >  1:1030767/1‑150 (MQ=255)
                                                                                                                           tGTAAAAACCGAGAGCGGTGGCTTAGTCAGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACg                           <  2:860946/150‑1 (MQ=255)
                                                                                                                           tGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAAcc                           >  2:1150886/1‑149 (MQ=255)
                                                                                                                            gTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGc                          >  1:447309/1‑150 (MQ=255)
                                                                                                                             tAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGcc                         >  1:707379/1‑150 (MQ=255)
                                                                                                                             tAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGcc                         >  1:810832/1‑150 (MQ=255)
                                                                                                                              aaaaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGCACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCa                        >  2:1105912/1‑150 (MQ=255)
                                                                                                                              aaaaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCa                        <  1:671401/150‑1 (MQ=255)
                                                                                                                                 aaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAtt                     <  2:1500577/150‑1 (MQ=255)
                                                                                                                                    cGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGt                  >  2:807430/1‑150 (MQ=255)
                                                                                                                                    cGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGt                  >  1:399221/1‑150 (MQ=255)
                                                                                                                                      taacaggTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGc                >  1:1229663/6‑150 (MQ=255)
                                                                                                                                           ggTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa           >  1:692802/1‑150 (MQ=255)
                                                                                                                                           ggTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTACCCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa           >  1:231441/1‑150 (MQ=255)
                                                                                                                                           ggTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa           >  2:441477/1‑150 (MQ=255)
                                                                                                                                           ggTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa           >  2:1367129/1‑150 (MQ=255)
                                                                                                                                               gCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACaaa       >  1:1309491/1‑150 (MQ=255)
                                                                                                                                                 ttAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATa     <  2:1030767/150‑1 (MQ=255)
                                                                                                                                                   aGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATaga   <  2:160016/150‑1 (MQ=255)
                                                                                                                                                   aGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATaga   <  1:907743/150‑1 (MQ=255)
                                                                                                                                                    gTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATagag  <  2:517842/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAG  >  NC_000913/4002249‑4002546

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: