Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,839 T→C 8.7% H57H (CAT→CAC yihW → putative transcriptional regulator for sulphoquinovose utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,074,8390TC 8.7% 272.5 / 15.2 105H57H (CAT→CACyihWputative transcriptional regulator for sulphoquinovose utilization
Reads supporting (aligned to +/- strand):  ref base T (41/52);  new base C (3/6);  total (45/60)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-01

Only 100 of 106 total aligned reads displayed.

GTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAC  >  NC_000913/4074691‑4074988
                                                                                                                                                    |                                                                                                                                                     
gTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATc                                                                                                                                                      <  1:605712/150‑1 (MQ=255)
  aaCCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCAc                                                                                                                                                    <  1:984875/150‑1 (MQ=255)
    cccGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCAGCACgg                                                                                                                                                  >  2:611721/1‑150 (MQ=255)
     ccGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGt                                                                                                                                                 >  2:1696323/1‑150 (MQ=255)
      cgcgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTg                                                                                                                                                >  2:1443121/1‑150 (MQ=255)
      cgcgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTg                                                                                                                                                <  1:973259/150‑1 (MQ=255)
       gcgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTgg                                                                                                                                               <  1:1139525/150‑1 (MQ=255)
        cgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGgc                                                                                                                                              <  1:1282348/150‑1 (MQ=255)
           cACGACCAACTCCTCATGCTGATCGCCGATCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGcgcg                                                                                                                                           <  2:763430/150‑1 (MQ=255)
            aCGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAACTTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgt                                                                                                                                          <  1:1088793/150‑2 (MQ=255)
             cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCggg                                                                                                                                         >  1:1699937/1‑150 (MQ=255)
             cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCggg                                                                                                                                         <  1:1127946/150‑1 (MQ=255)
             cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgta                                                                                                                                         <  2:1659138/150‑3 (MQ=255)
              gACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGt                                                                                                                                        >  1:1342710/1‑150 (MQ=255)
              gACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaa                                                                                                                                        <  1:1377094/150‑4 (MQ=255)
                  aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                    <  1:607969/150‑1 (MQ=255)
                  aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                    <  2:1452306/150‑1 (MQ=255)
                  aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaatcag                                                                                                                                    <  1:1293588/150‑8 (MQ=255)
                   aCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGc                                                                                                                                   <  1:1160180/150‑1 (MQ=255)
                      cctcATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTc                                                                                                                                <  2:358766/150‑1 (MQ=255)
                       ctcATGCTGATCGCCGAGCGTGGGTATATGAATATTGCTGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTcc                                                                                                                               <  1:14590/150‑1 (MQ=255)
                         cATGCTTATCGCCGAGCGGGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGACTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAg                                                                                                                             <  1:117392/150‑1 (MQ=255)
                         cATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAg                                                                                                                             <  1:495159/150‑1 (MQ=255)
                           tGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaatcag                                                                                                                                    <  2:171665/141‑8 (MQ=255)
                           tGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaatcag                                                                                                                                    >  1:171665/1‑134 (MQ=255)
                             cTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtc                                                                                                                         >  1:1630471/1‑150 (MQ=255)
                              tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                        >  2:252671/1‑150 (MQ=255)
                              tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                        <  1:1161043/150‑1 (MQ=255)
                              tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                        >  2:359600/1‑150 (MQ=255)
                                   gCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCCGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCTAGGCCTGAGTACGCGCCCTCACGGTGGCGCGGGTCGGGCTTCCAGCGTAGTTAAt                                                                                                                   >  2:1230968/1‑150 (MQ=255)
                                   gCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAAt                                                                                                                   >  2:1013964/1‑150 (MQ=255)
                                      gAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACg                                                                                                                <  2:762830/150‑1 (MQ=255)
                                       aGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCCCGCAGACGGACCGCCGGGATATTCGTAAATTAAGCGATCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGGCGGGCTTCCAGCGTCGTTATTACgg                                                                                                               >  2:1467294/1‑150 (MQ=255)
                                           tGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGtt                                                                                                           >  2:208126/1‑150 (MQ=255)
                                           tGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGAAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGtt                                                                                                           >  1:1303855/1‑150 (MQ=255)
                                            gggTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTc                                                                                                          <  1:1013964/150‑1 (MQ=255)
                                              gTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGCCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGa                                                                                                        >  2:1396748/1‑150 (MQ=255)
                                               tataTGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGag                                                                                                       >  2:745902/1‑150 (MQ=255)
                                               taaatGAATATTGACGAGCTGGCAAAACTGCTGGATGGCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCCCGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGag                                                                                                       <  1:1578442/147‑1 (MQ=255)
                                                 tatGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGagca                                                                                                     >  2:637897/1‑150 (MQ=255)
                                                   tGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGagcagc                                                                                                   >  1:457126/1‑150 (MQ=255)
                                                    gAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCg                                                                                                  <  1:745902/150‑1 (MQ=255)
                                                      atatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGGTAATACGGCGTTCGAGCAGCGTg                                                                                                >  2:153512/1‑150 (MQ=255)
                                                       tatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCgaga                                                                                               >  2:770864/1‑150 (MQ=255)
                                                       tatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGa                                                                                               >  2:1187110/1‑150 (MQ=255)
                                                       tatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGa                                                                                               >  2:605262/1‑150 (MQ=255)
                                                        tttGATTAGCTGGCAAATATGCTGGAGGTCTCCAGGCAGACGGTCCTCCGGGATACTCGTAAATTAAGCGAGCAAGGCCTGTTTACGCTCCAGCAGTGTTGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAg                                                                                              <  1:544455/149‑1 (MQ=255)
                                                          tgatgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCAGCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGt                                                                                            <  1:1581662/150‑1 (MQ=255)
                                                          tgatgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGCGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCCCGGTGGCGCGGGTCGGGCTTCCAGAGTCGTTAAAACGGCGTTCGAGCCGCGTGAGGt                                                                                            >  1:298705/1‑150 (MQ=255)
                                                           gatgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGtt                                                                                           >  1:609743/1‑150 (MQ=255)
                                                             tgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGTAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTc                                                                                         >  1:836378/1‑150 (MQ=255)
                                                              gaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCg                                                                                        <  2:69552/150‑1 (MQ=255)
                                                              gaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCg                                                                                        >  2:297936/1‑150 (MQ=255)
                                                               cgCTGGCAACTGTGCTGGCTCTCTCCCCGCAGACGGTCCGCCGGGATATTCGTAACTTACGCGAGCAATGCCTGATTCCGCGCCATCACGGCGGCGCGTGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGc                                                                                       <  2:289583/149‑1 (MQ=255)
                                                               aGCTGGCAAATCTGCTGGATGGCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTACGCGAGCAAGGCCTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGc                                                                                       <  1:1292811/150‑1 (MQ=255)
                                                                gCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCa                                                                                      <  2:1472050/150‑1 (MQ=255)
                                                                   ggCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCAATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAAc                                                                                   >  2:307333/1‑150 (MQ=255)
                                                                       aaTCTGCTGGATGTCTCCACGCAGACGGGCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAg                                                                               <  1:1388679/150‑1 (MQ=255)
                                                                         tCTGCGGGATGGCGCCACGCAGACGCTCCTCCGGGATATTCGTAAAATCAGCGAGCCAGGCCTGATTACGCGCCATCTCGGTGGCGCGGGTCGGGCTTCCAGCGTCGTAAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGa                                                                             <  1:431882/150‑1 (MQ=255)
                                                                           tgctgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaa                                                                           <  1:917368/150‑1 (MQ=255)
                                                                                 aTGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAGAAAAAAg                                                                     <  1:1696323/150‑1 (MQ=255)
                                                                                       ccACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTg                                                               <  2:1342710/150‑1 (MQ=255)
                                                                                         aCGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGcc                                                             <  1:1458310/150‑1 (MQ=255)
                                                                                          cGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCg                                                            <  1:605262/150‑1 (MQ=255)
                                                                                           gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                           <  1:208126/150‑1 (MQ=255)
                                                                                           gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                           <  2:530097/150‑1 (MQ=255)
                                                                                           gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGACCAAGGCCTGAGTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGGCGTTAATACGGCGTGCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                           <  2:134248/150‑1 (MQ=255)
                                                                                                 ggTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGt                                                     <  1:297936/150‑1 (MQ=255)
                                                                                                   tCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTgg                                                   <  1:1531933/150‑1 (MQ=255)
                                                                                                     cgccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCa                                                 <  2:1630471/150‑1 (MQ=255)
                                                                                                        ctggATTTTCTTAAATTAAGCGATCAAGGCCTGTTTACGCGCCATCACTGTGGCGCGGGTCGGGCTTCCAGCGTGGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACGGAGGAATAAAAATCGTTTGCCGAAGCGGTGGCAGAc                                              <  2:694515/148‑1 (MQ=255)
                                                                                                           gATATTCGTAAATAAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACtat                                           <  2:867962/150‑1 (MQ=255)
                                                                                                              atTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGGCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATAtt                                        <  2:363007/150‑1 (MQ=255)
                                                                                                                   tAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGa                                   >  1:274686/1‑150 (MQ=255)
                                                                                                                   tAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGa                                   >  1:1582219/1‑150 (MQ=255)
                                                                                                                    aaaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGAt                                  >  1:607138/1‑150 (MQ=255)
                                                                                                                     aaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATg                                 >  2:1205585/1‑150 (MQ=255)
                                                                                                                     aaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATg                                 >  2:1211583/1‑150 (MQ=255)
                                                                                                                          tgCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTca                            <  1:1638065/149‑1 (MQ=255)
                                                                                                                            cGAGCAAGGCCTGATTACGCGCCATCAGGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                          <  2:1640708/150‑1 (MQ=255)
                                                                                                                            cGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                          <  1:829732/150‑1 (MQ=255)
                                                                                                                            cGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                          >  1:525212/1‑150 (MQ=255)
                                                                                                                             gAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaaca                         >  1:701742/1‑150 (MQ=255)
                                                                                                                              aGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaacaa                        >  2:195958/1‑150 (MQ=255)
                                                                                                                                cAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAata                      <  1:359600/150‑1 (MQ=255)
                                                                                                                                cAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAata                      <  1:252671/150‑1 (MQ=255)
                                                                                                                                   ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCACACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAttt                   >  2:162225/1‑150 (MQ=255)
                                                                                                                                   ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAttt                   >  2:941250/1‑150 (MQ=255)
                                                                                                                                    cccTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTa                  <  1:1533384/149‑1 (MQ=255)
                                                                                                                                     ccTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTAt                 >  1:1303390/1‑150 (MQ=255)
                                                                                                                                      cTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATc                >  2:587297/1‑150 (MQ=255)
                                                                                                                                      cTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATc                >  2:760177/1‑150 (MQ=255)
                                                                                                                                       tGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCa               >  2:483368/1‑150 (MQ=255)
                                                                                                                                       ccgacccgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGACGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCa               >  1:1355332/7‑150 (MQ=255)
                                                                                                                                            aCGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCAtt          <  1:143186/150‑1 (MQ=255)
                                                                                                                                             cgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg         >  1:1672749/1‑150 (MQ=255)
                                                                                                                                              gcgcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTgg        <  2:1699937/150‑1 (MQ=255)
                                                                                                                                              gcgcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTgg        <  2:701742/150‑1 (MQ=255)
                                                                                                                                                   aTCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCCATCCGAGGGAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGa   >  1:279495/1‑150 (MQ=255)
                                                                                                                                                    tCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAc  <  1:637897/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAC  >  NC_000913/4074691‑4074988

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