Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 207,483 A→C 6.4% V786V (GTA→GTC dnaE → DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913207,4830AC 6.4% 231.9 / 7.5 78V786V (GTA→GTCdnaEDNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base C (2/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

ACGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATG  >  NC_000913/207335‑207625
                                                                                                                                                    |                                                                                                                                              
aCGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTaa                                                                                                                                               >  2:1018799/1‑150 (MQ=255)
    tGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                                                                           >  1:1613579/1‑150 (MQ=255)
     gAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAcc                                                                                                                                          >  1:1234943/1‑150 (MQ=255)
             gATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                                                                  <  2:520221/150‑5 (MQ=255)
             gATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                                                                  <  2:399315/150‑5 (MQ=255)
              aTGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccata                                                                                                                                 <  2:1382524/150‑6 (MQ=255)
                gAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                                                               >  1:1658645/1‑150 (MQ=255)
                gAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                                                               >  1:1421535/1‑150 (MQ=255)
                gAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                                                               >  1:1415774/1‑150 (MQ=255)
                 aaaaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGac                                                                                                                              <  2:1602117/150‑1 (MQ=255)
                             cTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAg                                                                                                                  >  2:132457/1‑150 (MQ=255)
                              tggtggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAgg                                                                                                                 >  1:615934/1‑150 (MQ=255)
                               ggtggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACTTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGTTAATGACCGCCGATATGGACAACACCGAGAAggt                                                                                                                <  1:179708/150‑1 (MQ=255)
                               ggtggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggt                                                                                                                <  1:1018799/150‑1 (MQ=255)
                                  ggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCga                                                                                                                      >  2:1512090/1‑141 (MQ=255)
                                  ggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCga                                                                                                                      <  1:1512090/141‑1 (MQ=255)
                                  ggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtggt                                                                                                             <  2:53573/150‑1 (MQ=255)
                                  ggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggt                                                                                                                >  2:1627665/1‑147 (MQ=255)
                                  ggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggt                                                                                                                <  1:1627665/147‑1 (MQ=255)
                                     taaaTTCGCTGGTTACTGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATTGCTTAAAGCGCACTATCCTGCGGAGTTTATGGCGTCGGTAATGACCGCCGATATGGACAACACCGAGAATGTGGTggg                                                                                                          <  2:1387107/149‑1 (MQ=255)
                                           cGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggt                                                                                                    >  2:1115486/1‑150 (MQ=255)
                                            gCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggag                                                                                                   <  2:1259797/150‑3 (MQ=255)
                                            gCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTGTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtg                                                                                                   <  2:764529/150‑1 (MQ=255)
                                             cTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                                                                  >  2:588177/1‑150 (MQ=255)
                                             cTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                                                                  <  2:689137/150‑1 (MQ=255)
                                             cTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                                                                  <  2:1133553/150‑1 (MQ=255)
                                                gTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATg                                                                                               >  2:958606/1‑150 (MQ=255)
                                                gTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATg                                                                                               <  1:1211167/150‑1 (MQ=255)
                                                    cGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTg                                                                                           >  1:1182735/1‑150 (MQ=255)
                                                      gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATCGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCt                                                                                         >  1:870488/1‑150 (MQ=255)
                                                          tAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCg                                                                                     >  2:681289/1‑150 (MQ=255)
                                                           aaCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtg                                                                                                   >  1:1183834/1‑135 (MQ=255)
                                                           aaCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtg                                                                                                   <  2:1183834/135‑1 (MQ=255)
                                                              aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATg                                                                                 <  1:1285302/150‑1 (MQ=255)
                                                                 tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgggg                                                                              <  2:1613579/150‑1 (MQ=255)
                                                                  cGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGgtggtg                                                                                                            <  2:1383997/119‑1 (MQ=255)
                                                                  cGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGgtggtg                                                                                                            >  1:1383997/1‑119 (MQ=255)
                                                                  cGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                                                             <  1:364765/150‑1 (MQ=255)
                                                                     aCTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGa                                                                          >  1:1474678/1‑150 (MQ=255)
                                                                     aCTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGa                                                                          <  1:288320/150‑1 (MQ=255)
                                                                       tctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGaaa                                                                        >  1:1030548/1‑150 (MQ=255)
                                                                        ctGCGGCCTATGCTTTTGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGTTGGGGCTTaaaa                                                                       <  1:915578/150‑1 (MQ=255)
                                                                         tGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAAt                                                                      >  2:920457/1‑150 (MQ=255)
                                                                           cGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATcc                                                                    <  2:862914/150‑1 (MQ=255)
                                                                            ggCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCt                                                                   >  1:1192642/1‑150 (MQ=255)
                                                                               cTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGTAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGcc                                                                <  2:1493624/150‑1 (MQ=255)
                                                                                tATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGcca                                                               >  1:422216/1‑150 (MQ=255)
                                                                                  tGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacc                                                             >  2:1264548/1‑150 (MQ=255)
                                                                                   gCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacca                                                            <  2:813378/150‑1 (MQ=255)
                                                                                    cTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg                                                           >  2:1375275/1‑150 (MQ=255)
                                                                                    cTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg                                                           <  1:503334/150‑1 (MQ=255)
                                                                                        ggTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGatat                                                       <  2:390628/150‑1 (MQ=255)
                                                                                         gtgtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata                                                      <  1:95424/150‑1 (MQ=255)
                                                                                          tgtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATaa                                                     >  2:1529840/1‑150 (MQ=255)
                                                                                            tCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAAc                                                   >  1:59428/1‑150 (MQ=255)
                                                                                             cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACt                                                  >  1:808397/1‑150 (MQ=255)
                                                                                               tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTcc                                                >  1:280094/1‑150 (MQ=255)
                                                                                                 tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg                                              >  1:1080024/1‑150 (MQ=255)
                                                                                                 tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg                                              <  2:280094/150‑1 (MQ=255)
                                                                                                 tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGCGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg                                              >  2:392137/1‑150 (MQ=255)
                                                                                                  cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGt                                             >  2:1561153/1‑150 (MQ=255)
                                                                                                  cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGt                                             <  2:1206224/150‑1 (MQ=255)
                                                                                                  cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGt                                             >  2:826010/1‑150 (MQ=255)
                                                                                                            ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                   <  1:767828/150‑1 (MQ=255)
                                                                                                              cTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                   >  1:1238549/1‑148 (MQ=255)
                                                                                                              cTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                   <  2:1238549/148‑1 (MQ=255)
                                                                                                                gAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCa                               <  2:808397/150‑1 (MQ=255)
                                                                                                                  aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACg                             <  1:1115486/150‑1 (MQ=255)
                                                                                                                  aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACg                             >  2:1095868/1‑150 (MQ=255)
                                                                                                                    gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTc                           <  1:826010/150‑1 (MQ=255)
                                                                                                                       cACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAac                        >  2:343461/1‑150 (MQ=255)
                                                                                                                              cTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGAGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacgacgacg                 >  2:166997/1‑150 (MQ=255)
                                                                                                                               tGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGATAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACgg                <  1:588177/150‑1 (MQ=255)
                                                                                                                                       ccatatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATTCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCg        >  2:347079/7‑150 (MQ=255)
                                                                                                                                       ccatatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCg        >  1:1010797/7‑150 (MQ=255)
                                                                                                                                         aTGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtg      <  2:362453/150‑1 (MQ=255)
                                                                                                                                           ggcggcGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTa    <  1:343461/150‑1 (MQ=255)
                                                                                                                                             cggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTTGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATg  <  1:569591/150‑1 (MQ=255)
                                                                                                                                             cggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATg  >  1:262806/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATG  >  NC_000913/207335‑207625

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: