Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,841 A→G 8.7% H58R (CAC→CGC) ‡ 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,8410AG 8.7% 280.3 / 14.4 105H58R (CAC→CGC) ‡yihWputative transcriptional regulator for sulphoquinovose utilization
Reads supporting (aligned to +/- strand):  ref base A (41/52);  new base G (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.16e-01

Only 100 of 105 total aligned reads displayed.

AACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAC  >  NC_000913/4074693‑4074988
                                                                                                                                                    |                                                                                                                                                   
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)
          aCGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgt                                                                                                                                          <  1:787732/150‑2 (MQ=255)
           cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCggg                                                                                                                                         <  1:1127946/150‑1 (MQ=255)
           cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCggg                                                                                                                                         >  1:1699937/1‑150 (MQ=255)
           cGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgta                                                                                                                                         <  2:1659138/150‑3 (MQ=255)
            gACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGt                                                                                                                                        >  1:1342710/1‑150 (MQ=255)
            gACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaa                                                                                                                                        <  1:1377094/150‑4 (MQ=255)
                aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                    <  2:1452306/150‑1 (MQ=255)
                aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                    <  1:607969/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)
                         tGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcg                                                                                                                           <  1:293691/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:359600/1‑150 (MQ=255)
                            tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                        >  2:252671/1‑150 (MQ=255)
                            tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                        <  1:1161043/150‑1 (MQ=255)
                                 gCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCCGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCTAGGCCTGAGTACGCGCCCTCACGGTGGCGCGGGTCGGGCTTCCAGCGTAGTTAAt                                                                                                                   >  2:1230968/1‑150 (MQ=255)
                                 gCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAAt                                                                                                                   >  2:1013964/1‑150 (MQ=255)
                                   cGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATAc                                                                                                                 >  2:1531933/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:605262/1‑150 (MQ=255)
                                                     tatTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGa                                                                                               >  2:1187110/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:297936/1‑150 (MQ=255)
                                                            gaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCg                                                                                        <  2:69552/150‑1 (MQ=255)
                                                             cgCTGGCAACTGTGCTGGCTCTCTCCCCGCAGACGGTCCGCCGGGATATTCGTAACTTACGCGAGCAATGCCTGATTCCGCGCCATCACGGCGGCGCGTGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGc                                                                                       <  2:289583/149‑1 (MQ=255)
                                                             aGCTGGCAAATCTGCTGGATGGCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTACGCGAGCAAGGCCTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGc                                                                                       <  1:1292811/150‑1 (MQ=255)
                                                              gCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCa                                                                                      <  2:1472050/150‑1 (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)
                                                                                            gACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGc                                                        <  1:1443121/150‑1 (MQ=255)
                                                                                               ggTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGt                                                     <  1:297936/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)
                                                                                                                gTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                    <  2:425374/150‑1 (MQ=255)
                                                                                                                 tAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGa                                   >  1:1582219/1‑150 (MQ=255)
                                                                                                                 tAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGa                                   >  1:274686/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)
                                                                                                                                                    |                                                                                                                                                   
AACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAC  >  NC_000913/4074693‑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: