Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,400 T→A 6.6% intergenic (+24/+19) 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,4000TA 6.6% 429.9 / 14.5 136intergenic (+24/+19)corA/yigFNi(2(+))/Co(2(+))/Mg(2(+)) transporter/DUF2628 domain‑containing protein YigF
Reads supporting (aligned to +/- strand):  ref base T (63/64);  new base A (5/4);  total (68/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 138 total aligned reads displayed.

TATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGAT  >  NC_000913/4002251‑4002548
                                                                                                                                                     |                                                                                                                                                    
tATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCt                                                                                                                                                      >  1:512537/1‑150 (MQ=255)
   ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc                                                                                                                                                   <  1:1573785/150‑1 (MQ=255)
   ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCGGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc                                                                                                                                                   <  1:818048/150‑1 (MQ=255)
    gCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCt                                                                                                                                                  >  1:650140/1‑150 (MQ=255)
       tGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAg                                                                                                                                               <  1:1059753/150‑1 (MQ=255)
       tGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAg                                                                                                                                               >  1:298922/1‑150 (MQ=255)
       tGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAg                                                                                                                                               <  1:59955/150‑1 (MQ=255)
        gAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGc                                                                                                                                              >  2:672359/1‑150 (MQ=255)
         aaCTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGcc                                                                                                                                             <  1:115291/150‑1 (MQ=255)
           cTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAc                                                                                                                                           <  2:680256/150‑1 (MQ=255)
                 gTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa                                                                                                                                     <  1:62746/150‑1 (MQ=255)
                        ccAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAg                                                                                                                              <  2:29079/150‑1 (MQ=255)
                           gAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCt                                                                                                                           >  2:1445996/1‑150 (MQ=255)
                              cTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCa                                                                                                                        <  1:1507995/150‑1 (MQ=255)
                               tGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAg                                                                                                                       <  2:1049035/150‑1 (MQ=255)
                                gAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAgg                                                                                                                      <  1:1216161/150‑1 (MQ=255)
                                    tGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTAtt                                                                                                                  <  1:1467921/150‑1 (MQ=255)
                                     ggAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTc                                                                                                                                  <  2:426054/133‑1 (MQ=255)
                                     ggAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTc                                                                                                                                  >  1:426054/1‑133 (MQ=255)
                                             ggCTACCCTGGCGCGATTATCTTTAAGATACTCGCGGGCCTGGCACCGAATCTGTACTTTAACCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCttt                                                                                                         >  2:1263025/1‑150 (MQ=255)
                                              gCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCtttt                                                                                                        <  2:235831/150‑1 (MQ=255)
                                               cTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGGTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGGTATTCAAAGGCTCCAGGTATTTAACCCTTTTa                                                                                                       <  2:56627/150‑1 (MQ=255)
                                                tACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTAc                                                                                                      >  1:1563613/1‑150 (MQ=255)
                                                  cccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACct                                                                                                    >  1:1385359/1‑150 (MQ=255)
                                                  cccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACct                                                                                                    <  1:30876/150‑1 (MQ=255)
                                                   ccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctc                                                                                                   >  1:132774/1‑150 (MQ=255)
                                                      ggCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCttt                                                                                                >  2:1068360/1‑150 (MQ=255)
                                                       gcgcGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCGGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTtc                                                                                               <  1:776117/150‑1 (MQ=255)
                                                        cgcgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTtct                                                                                              <  1:548794/150‑1 (MQ=255)
                                                          cgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCa                                                                                            >  1:917045/1‑150 (MQ=255)
                                                            aTTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATa                                                                                          <  1:169673/150‑1 (MQ=255)
                                                             ttATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAg                                                                                         <  2:1353651/150‑1 (MQ=255)
                                                              tATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCCCCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGa                                                                                        >  1:1615183/1‑150 (MQ=255)
                                                               aTCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGaa                                                                                       >  2:637360/1‑150 (MQ=255)
                                                                tCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTTGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAAc                                                                                      <  1:732775/150‑1 (MQ=255)
                                                                   ttATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAt                                                                                   >  1:343442/1‑150 (MQ=255)
                                                                    tATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAtt                                                                                  >  1:1183971/1‑150 (MQ=255)
                                                                     aTGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAttt                                                                                 >  1:11220/1‑150 (MQ=255)
                                                                      tGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTg                                                                                >  2:499949/1‑150 (MQ=255)
                                                                       gATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGt                                                                               >  2:306737/1‑150 (MQ=255)
                                                                       gATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGt                                                                               >  1:1500518/1‑150 (MQ=255)
                                                                          ccTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCg                                                                            <  1:1539191/150‑1 (MQ=255)
                                                                            tCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCGGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCgtg                                                                          <  2:351059/150‑1 (MQ=255)
                                                                             cgcgGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCgtgt                                                                         >  1:464632/1‑150 (MQ=255)
                                                                              gcgGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGtt                                                                        >  2:563909/1‑150 (MQ=255)
                                                                                gggCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTaa                                                                      >  1:1364210/1‑150 (MQ=255)
                                                                                 ggCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAAc                                                                     >  2:562816/1‑150 (MQ=255)
                                                                                  gCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACa                                                                    <  1:781879/150‑1 (MQ=255)
                                                                                  gCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACa                                                                    <  2:1270669/150‑1 (MQ=255)
                                                                                    cTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGc                                                                  <  1:1254673/150‑1 (MQ=255)
                                                                                    cTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGc                                                                  <  1:1260524/150‑1 (MQ=255)
                                                                                    cTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGc                                                                  >  1:1363390/1‑150 (MQ=255)
                                                                                      ggCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCaa                                                                >  2:1108775/1‑150 (MQ=255)
                                                                                      ggCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCaa                                                                >  2:167997/1‑150 (MQ=255)
                                                                                          ccGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATa                                                            >  1:1419824/1‑150 (MQ=255)
                                                                                          ccGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATa                                                            <  2:1205710/150‑1 (MQ=255)
                                                                                              aTCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgct                                                        <  2:404899/150‑1 (MQ=255)
                                                                                              aTCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgct                                                        >  1:805319/1‑150 (MQ=255)
                                                                                               tCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctg                                                       <  1:1000162/150‑1 (MQ=255)
                                                                                                 tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgct                                                     >  1:1231029/1‑150 (MQ=255)
                                                                                                 tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgct                                                     >  1:426824/1‑150 (MQ=255)
                                                                                                 tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgct                                                     >  1:396121/1‑150 (MQ=255)
                                                                                                 tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgct                                                     <  2:1324704/150‑1 (MQ=255)
                                                                                                     cTTTAAGCTGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACgg                                                 >  1:527746/1‑150 (MQ=255)
                                                                                                     cTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACgg                                                 >  1:1285026/1‑150 (MQ=255)
                                                                                                      tttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                                <  2:3822/150‑1 (MQ=255)
                                                                                                      tttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                                >  2:222292/1‑150 (MQ=255)
                                                                                                      tttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                                <  2:1363390/150‑1 (MQ=255)
                                                                                                        tAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCaa                                              <  1:306737/150‑1 (MQ=255)
                                                                                                               aagaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAACGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAg                                       <  2:1217735/150‑1 (MQ=255)
                                                                                                               aagaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAg                                       <  2:1385359/150‑1 (MQ=255)
                                                                                                                   aCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATg                                   <  1:1584367/150‑1 (MQ=255)
                                                                                                                   aCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATg                                   >  1:1411137/1‑150 (MQ=255)
                                                                                                                   aCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTCATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATg                                   <  1:455540/150‑1 (MQ=255)
                                                                                                                      ggTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTa                                <  1:672359/150‑1 (MQ=255)
                                                                                                                         tGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACg                             <  2:1500518/150‑1 (MQ=255)
                                                                                                                          gTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGc                            >  1:396288/1‑150 (MQ=255)
                                                                                                                            aaaaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCa                          >  2:1473887/1‑150 (MQ=255)
                                                                                                                             aaaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCaa                         >  2:766474/1‑150 (MQ=255)
                                                                                                                              aaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAt                        >  1:561719/1‑150 (MQ=255)
                                                                                                                               aaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAtt                       <  1:1287234/150‑1 (MQ=255)
                                                                                                                                 aCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAg                     <  1:1115437/150‑1 (MQ=255)
                                                                                                                                 aCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAg                     <  2:396288/150‑1 (MQ=255)
                                                                                                                                   gagaGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTg                   <  2:917045/150‑1 (MQ=255)
                                                                                                                                    agagCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGc                  <  1:637360/150‑1 (MQ=255)
                                                                                                                                    agagCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGc                  >  1:1467330/1‑150 (MQ=255)
                                                                                                                                      agCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAAc                                                                     >  1:309621/1‑97 (MQ=255)
                                                                                                                                      agCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAAc                                                                     <  2:309621/97‑1 (MQ=255)
                                                                                                                                      acaggTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAAAGTTAACGCCAATTAGTGCAc                >  1:841787/4‑150 (MQ=255)
                                                                                                                                       caggTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCAcc               >  1:575278/3‑150 (MQ=255)
                                                                                                                                          gTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctg                                                       >  1:91216/1‑107 (MQ=255)
                                                                                                                                          gTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctg                                                       <  2:91216/107‑1 (MQ=255)
                                                                                                                                           tGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaaca           <  2:298922/150‑1 (MQ=255)
                                                                                                                                           tGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaaca           <  2:343442/150‑1 (MQ=255)
                                                                                                                                            ggCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaacaa          >  1:1618028/1‑150 (MQ=255)
                                                                                                                                            ggCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaacaa          >  1:930163/1‑150 (MQ=255)
                                                                                                                                                 aGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATaga     <  1:1068360/150‑1 (MQ=255)
                                                                                                                                                 aGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGACAATTAGTGCACCAACAAATaga     >  1:744192/1‑150 (MQ=255)
                                                                                                                                                  gTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATagag    <  2:650140/150‑1 (MQ=255)
                                                                                                                                                    cTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGAt  >  2:507593/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGAT  >  NC_000913/4002251‑4002548

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: