Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,971 G→A 5.2% Q259* (CAG→TAG) ‡ yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,681,9710GA 5.2% 377.3 / 6.6 116Q259* (CAG→TAG) ‡yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (56/53);  new base A (2/4);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 119 total aligned reads displayed.

TTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681823‑2682118
                                                                                                                                                    |                                                                                                                                                   
ttATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTg                                                                                                                                                     >  3:864802/1‑149 (MQ=255)
 tATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGc                                                                                                                                                    <  2:67051/149‑1 (MQ=255)
  aTTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCg                                                                                                                                                   <  4:691125/149‑1 (MQ=255)
     gTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtt                                                                                                                                                <  4:326204/149‑1 (MQ=255)
     gTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAtt                                                                                                                                                <  1:164731/149‑1 (MQ=255)
        ggTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAttttt                                                                                                                                             >  3:731946/1‑149 (MQ=255)
         gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGAAAAATCGCAGTACGAttttgtc                                                                                                                                           <  2:517929/149‑4 (MQ=255)
          tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCg                                                                                                                                           >  1:795729/1‑149 (MQ=255)
          tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCg                                                                                                                                           <  2:56403/149‑1 (MQ=255)
              aCTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCa                                                                                                                                       >  2:352851/1‑149 (MQ=255)
              aCTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCa                                                                                                                                       >  4:784067/1‑149 (MQ=255)
                tctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCata                                                                                                                                     >  4:438053/1‑149 (MQ=255)
                tctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCata                                                                                                                                     >  2:401856/1‑149 (MQ=255)
                tctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCata                                                                                                                                     >  3:432332/1‑149 (MQ=255)
                tctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCata                                                                                                                                     >  2:596587/1‑149 (MQ=255)
                actTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCata                                                                                                                                     >  1:170863/2‑149 (MQ=255)
                  ttATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCata                                                                                                                                     >  4:797582/1‑147 (MQ=255)
                  ttATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCata                                                                                                                                     <  3:797582/147‑1 (MQ=255)
                    aTCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATg                                                                                                                                 >  1:620210/1‑149 (MQ=255)
                    aTCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATg                                                                                                                                 >  1:806851/1‑149 (MQ=255)
                     tCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCATCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGa                                                                                                                                >  1:47798/1‑149 (MQ=255)
                      cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAg                                                                                                                               <  2:461153/149‑1 (MQ=255)
                        cACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGt                                                                                                                             >  1:648733/1‑149 (MQ=255)
                         gcGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACGCCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGtt                                                                                                                            <  3:89052/148‑1 (MQ=255)
                         aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGtt                                                                                                                            <  1:274270/149‑1 (MQ=255)
                         aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGtt                                                                                                                            >  3:218637/1‑149 (MQ=255)
                             cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGaa                                                                                                                        <  1:308101/149‑1 (MQ=255)
                             cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGaa                                                                                                                        <  4:218637/149‑1 (MQ=255)
                             cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGaa                                                                                                                        <  4:731946/149‑1 (MQ=255)
                                 ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGAAAAATCGCAGTACGAttttgtc                                                                                                                                           >  4:361777/1‑122 (MQ=255)
                                 ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGAAAAATCGCAGTACGAttttgtc                                                                                                                                           <  3:361777/125‑4 (MQ=255)
                                   aaaTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGca                                                                                                                  <  4:597446/149‑1 (MQ=255)
                                   aaaTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGca                                                                                                                  >  1:233723/1‑149 (MQ=255)
                                     aTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAtttt                                                                                                                                              <  4:483453/119‑1 (MQ=255)
                                     aTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAtttt                                                                                                                                              >  3:483453/1‑119 (MQ=255)
                                         aCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAtttt                                                                                                            >  4:100256/1‑149 (MQ=255)
                                           tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTta                                                                                                          <  3:784067/149‑1 (MQ=255)
                                             tgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTtata                                                                                                        >  3:844628/1‑149 (MQ=255)
                                               tCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGc                                                                                                      <  3:484077/149‑1 (MQ=255)
                                                 cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    >  3:706057/1‑149 (MQ=255)
                                                 cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    <  3:744096/149‑1 (MQ=255)
                                                 cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    <  2:233723/149‑1 (MQ=255)
                                                 cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    >  2:31114/1‑149 (MQ=255)
                                                 cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    >  2:589554/1‑149 (MQ=255)
                                                 cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    >  4:744096/1‑149 (MQ=255)
                                                  ccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   >  3:809241/1‑149 (MQ=255)
                                                  ccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   >  4:582458/1‑149 (MQ=255)
                                                  ccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   >  3:671866/1‑149 (MQ=255)
                                                  ccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGACCATATATGAGGTTGGAACAGGCACATTTTTATATCgag                                                                                                   >  4:551068/1‑149 (MQ=255)
                                                    aaCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaa                                                                                                 >  1:1573/1‑149 (MQ=255)
                                                     aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaaa                                                                                                <  3:14408/149‑1 (MQ=255)
                                                         cACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTc                                                                                            <  3:452072/149‑1 (MQ=255)
                                                          aCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTcc                                                                                           >  1:170176/1‑149 (MQ=255)
                                                            tGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTc                                                                                            >  4:852560/1‑146 (MQ=255)
                                                            tGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCAc                                                                                         >  2:349126/1‑149 (MQ=255)
                                                            tGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCAc                                                                                         >  1:643593/1‑149 (MQ=255)
                                                            tGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCAc                                                                                         >  3:422070/1‑149 (MQ=255)
                                                                ttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATCGCGAGAAATGTCCACTCCa                                                                                     <  3:695013/149‑1 (MQ=255)
                                                                ttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCa                                                                                     <  3:163435/149‑1 (MQ=255)
                                                                ttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCa                                                                                     <  4:855600/149‑1 (MQ=255)
                                                                    aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgg                                                                                 <  4:517477/149‑1 (MQ=255)
                                                                      cTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                                <  3:20871/148‑1 (MQ=255)
                                                                      cTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                                >  4:20871/1‑148 (MQ=255)
                                                                           gcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTc                                                                          <  2:456170/149‑1 (MQ=255)
                                                                                   gtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  4:671866/149‑1 (MQ=255)
                                                                                    tggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgta                                                                 >  3:328189/1‑149 (MQ=255)
                                                                                    tggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgta                                                                 <  1:708904/149‑1 (MQ=255)
                                                                                        tGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTaa                                                             <  1:330277/149‑1 (MQ=255)
                                                                                         ggtaaCGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAg                                                            >  4:159391/4‑149 (MQ=255)
                                                                                          gCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGt                                                           >  1:359273/1‑149 (MQ=255)
                                                                                          gCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGt                                                           >  1:38717/1‑149 (MQ=255)
                                                                                          gCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGt                                                           >  1:180660/1‑149 (MQ=255)
                                                                                            aaCGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTGTGTAGTAAGTGc                                                         >  1:70412/1‑149 (MQ=255)
                                                                                              cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCg                                                       <  1:349126/149‑1 (MQ=255)
                                                                                                 gcagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTa                                                    <  4:422070/149‑1 (MQ=255)
                                                                                                   agcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                  <  3:87498/149‑1 (MQ=255)
                                                                                                        aCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgg                                             <  1:274665/149‑1 (MQ=255)
                                                                                                                cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATggg                                            >  2:778032/1‑142 (MQ=255)
                                                                                                                cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATggg                                            <  1:778032/142‑1 (MQ=255)
                                                                                                                 ttCATCGTGACACCACGCGCCGACAAAATCGTGCTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGATAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    <  3:551068/149‑1 (MQ=255)
                                                                                                                  tCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  1:614696/149‑1 (MQ=255)
                                                                                                                   cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgag                                  <  2:795729/149‑1 (MQ=255)
                                                                                                                    aTCGTGACACCCCGCGCCGACAAATTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGTGAAATGTCCACTCCAGCGGAGTAGTCCCCTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                 <  2:346764/149‑1 (MQ=255)
                                                                                                                    aTCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgg                                             <  3:738035/137‑1 (MQ=255)
                                                                                                                              cacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                       <  4:706057/149‑1 (MQ=255)
                                                                                                                               acGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  2:11755/136‑1 (MQ=255)
                                                                                                                               acGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   >  1:11755/1‑136 (MQ=255)
                                                                                                                               acGCGCCGACAAAATCGTACTGCGATTTTTCCGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTg                      >  2:857955/1‑149 (MQ=255)
                                                                                                                                         aaaaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa            >  2:336444/1‑149 (MQ=255)
                                                                                                                                         aaaaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa            >  1:820068/1‑149 (MQ=255)
                                                                                                                                         aaaaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa            >  3:653298/1‑149 (MQ=255)
                                                                                                                                          aaaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaaa           <  1:352851/149‑1 (MQ=255)
                                                                                                                                           aaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   >  1:95918/1‑124 (MQ=255)
                                                                                                                                           aaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  2:95918/124‑1 (MQ=255)
                                                                                                                                            aTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca         >  1:586416/1‑149 (MQ=255)
                                                                                                                                            aTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca         <  2:806851/149‑1 (MQ=255)
                                                                                                                                            aTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca         >  3:782672/1‑149 (MQ=255)
                                                                                                                                                   tGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  >  1:137064/1‑149 (MQ=255)
                                                                                                                                                    gCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   >  4:661083/1‑115 (MQ=255)
                                                                                                                                                    gCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  3:661083/115‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681823‑2682118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: