Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,969 C→G 5.1% Q259H (CAG→CAC) ‡ 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,9690CG 5.1% 337.4 / 6.7 99Q259H (CAG→CAC) ‡yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base C (40/53);  new base G (3/2);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 101 total aligned reads displayed.

CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681820‑2682118
                                                                                                                                                     |                                                                                                                                                      
ccGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTAc                                                                                                                                                        <  2:1177119/150‑1 (MQ=255)
  gTTATTGTCGGTCAGACTCTTATCCCATGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTg                                                                                                                                                      >  1:556848/1‑150 (MQ=255)
   ttATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGAGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGc                                                                                                                                                     <  1:506899/150‑1 (MQ=255)
       tGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtt                                                                                                                                                 >  1:1051072/1‑150 (MQ=255)
         tCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTTTCCCCAACTCCACTGTTTTTTACCTGGAGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtttt                                                                                                                                               >  1:1072853/1‑150 (MQ=255)
          cGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtttt‑t                                                                                                                                             >  2:1521984/1‑150 (MQ=255)
          cGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttg                                                                                                                                              <  1:1648921/150‑2 (MQ=255)
            gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGACATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCg                                                                                                                                           >  2:1118232/1‑150 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCTAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCgg                                                                                                                                          <  1:266200/150‑1 (MQ=255)
               aGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGcc                                                                                                                                        >  2:141749/1‑150 (MQ=255)
                 aCTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCat                                                                                                                                      <  1:642107/150‑1 (MQ=255)
                    ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTGTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCatata                                                                                                                                   <  1:1702961/150‑1 (MQ=255)
                      tATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATg                                                                                                                                 <  1:885551/150‑1 (MQ=255)
                      tATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATg                                                                                                                                 >  1:927240/1‑150 (MQ=255)
                          ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCg                                                                                                                                           <  1:758319/136‑1 (MQ=255)
                          ccACGCCTGGCCAAAATCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCg                                                                                                                                           >  2:758319/1‑136 (MQ=255)
                            aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                           >  2:1646368/1‑150 (MQ=255)
                              gCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATAGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGa                                                                                                                         <  1:1549698/150‑1 (MQ=255)
                               ccTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGaa                                                                                                                        <  2:1247712/150‑1 (MQ=255)
                                cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAAc                                                                                                                       >  1:1387327/1‑150 (MQ=255)
                                  ggCCAAATTCACTGTGTACCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAg                                                                                                                     >  1:1081345/1‑150 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAgg                                                                                                                    <  1:141749/150‑1 (MQ=255)
                                    ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   <  2:1072853/150‑1 (MQ=255)
                                    ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   <  1:204425/150‑1 (MQ=255)
                                    ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   >  1:324737/1‑150 (MQ=255)
                                      aaaTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGcac                                                                                                                 >  1:1382272/1‑150 (MQ=255)
                                        aTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACAt                                                                                                               >  1:665575/1‑150 (MQ=255)
                                            aCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACAttttt                                                                                                           <  1:556781/150‑1 (MQ=255)
                                              tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  1:1349627/150‑1 (MQ=255)
                                              tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  2:1296516/150‑1 (MQ=255)
                                              tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  1:1296837/150‑1 (MQ=255)
                                               gtgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGGGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtata                                                                                                        <  1:434919/150‑1 (MQ=255)
                                                  tCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCg                                                                                                     >  1:303970/1‑150 (MQ=255)
                                                  tCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCg                                                                                                     >  1:617118/1‑150 (MQ=255)
                                                   ccccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    <  2:1580380/150‑1 (MQ=255)
                                                    cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  2:1444292/150‑1 (MQ=255)
                                                      cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAATCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaa                                                                                                 >  2:1117016/1‑150 (MQ=255)
                                                      cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCTACAAAATCGTAATGCGATTTT‑TCGGCAATTTATGAGGTTTGAACAGGCACATTTTTATAGCGAGaa                                                                                                 >  2:1532124/1‑150 (MQ=255)
                                                       acctaaaCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCTATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGATaaa                                                                                                <  2:347314/144‑1 (MQ=255)
                                                        aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAAt                                                                                               <  2:1382272/150‑1 (MQ=255)
                                                           ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTc                                                                                            >  2:43214/1‑150 (MQ=255)
                                                            cACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTcc                                                                                           >  2:1416112/1‑150 (MQ=255)
                                                                 ttttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTcc                                                                                      <  1:929605/150‑1 (MQ=255)
                                                                  tttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCa                                                                                     >  1:1479475/1‑150 (MQ=255)
                                                                      tACCTGGCGCAATATGGTGGGTGGCAACGTGCAGCCAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgg                                                                                 <  1:225913/150‑1 (MQ=255)
                                                                       aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                                >  1:978552/1‑150 (MQ=255)
                                                                       aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                                >  1:971955/1‑150 (MQ=255)
                                                                         cTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagt                                                                              >  1:1253729/1‑150 (MQ=255)
                                                                             cgcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTc                                                                          <  2:556848/150‑1 (MQ=255)
                                                                                 atatGGGGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                      <  2:324737/150‑1 (MQ=255)
                                                                                    tggtggTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTg                                                                   <  2:1253729/150‑1 (MQ=255)
                                                                                     ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  2:1632667/150‑1 (MQ=255)
                                                                                     ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  1:985708/150‑1 (MQ=255)
                                                                                       tggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtag                                                                <  1:1434751/150‑1 (MQ=255)
                                                                                         ggtGGCCACGTGCAGCAAACCGCCATCTTCATCGTGACACCCCGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtagta                                                              <  1:698900/148‑1 (MQ=255)
                                                                                             gCAACTTGCAGCAAACCGCAATCTTCATCCTGACACCCCGCGCCGACAAAATCCTACTGCGATTTT‑TCGGCCATATATGATGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCAACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTg                                                          <  2:508664/150‑1 (MQ=255)
                                                                                                 cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                      <  2:733187/150‑1 (MQ=255)
                                                                                                 cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                      <  2:1148008/150‑1 (MQ=255)
                                                                                                 cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                      <  2:1308550/150‑1 (MQ=255)
                                                                                                   tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGTCCATATATGAGGTTGGAACAGGCACAATTTAATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGTTa                                                    >  2:919324/1‑150 (MQ=255)
                                                                                                   tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTa                                                    >  2:205885/1‑150 (MQ=255)
                                                                                                   tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTa                                                    >  2:515173/1‑150 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     >  1:1258358/1‑131 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     <  2:1258358/131‑1 (MQ=255)
                                                                                                         aaaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATg                                              >  1:763731/1‑150 (MQ=255)
                                                                                                         aaaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATg                                              <  1:43214/150‑1 (MQ=255)
                                                                                                           aCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATggg                                            <  2:665575/150‑1 (MQ=255)
                                                                                                                  tCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCg                                     >  1:758369/1‑150 (MQ=255)
                                                                                                                   cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    <  1:5275/150‑1 (MQ=255)
                                                                                                                   cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    <  1:1117016/150‑1 (MQ=255)
                                                                                                                    ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  2:1400381/150‑1 (MQ=255)
                                                                                                                    ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  1:1405681/150‑1 (MQ=255)
                                                                                                                    ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGGCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  1:453561/150‑1 (MQ=255)
                                                                                                                      cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                 >  2:224320/1‑150 (MQ=255)
                                                                                                                          gTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                             >  2:477739/1‑150 (MQ=255)
                                                                                                                            gACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGccc                           >  2:226174/1‑150 (MQ=255)
                                                                                                                            gACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGccc                           <  1:977538/150‑1 (MQ=255)
                                                                                                                               accaaGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTT‑‑TAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCAGTTTg                      <  1:1118232/150‑1 (MQ=255)
                                                                                                                                ccacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                       >  1:721928/1‑150 (MQ=255)
                                                                                                                                ccacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                       >  1:1653718/1‑150 (MQ=255)
                                                                                                                                   cgcgCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                    >  1:922723/1‑150 (MQ=255)
                                                                                                                                    gcgcCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCg                   >  1:775811/1‑150 (MQ=255)
                                                                                                                                    gcgcCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCg                   <  2:927240/150‑1 (MQ=255)
                                                                                                                                     cgcCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGt                  <  2:617118/150‑1 (MQ=255)
                                                                                                                                      ggccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGt                  >  1:764754/7‑150 (MQ=255)
                                                                                                                                       gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                 >  2:916682/6‑150 (MQ=255)
                                                                                                                                       gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                 >  1:1389524/6‑150 (MQ=255)
                                                                                                                                       ccGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGa                >  1:1469116/1‑150 (MQ=255)
                                                                                                                                        cGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt               >  1:247755/1‑150 (MQ=255)
                                                                                                                                         cgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt               <  2:1389524/147‑1 (MQ=255)
                                                                                                                                           cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                           >  2:52864/1‑119 (MQ=255)
                                                                                                                                           cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                           <  1:52864/119‑1 (MQ=255)
                                                                                                                                            aaaaTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaaa           <  1:552089/150‑1 (MQ=255)
                                                                                                                                              aaTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca         <  1:1521984/150‑1 (MQ=255)
                                                                                                                                              aaTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca         >  2:966359/1‑150 (MQ=255)
                                                                                                                                               aTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAacac        <  1:1646368/150‑1 (MQ=255)
                                                                                                                                                  gTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  2:1387327/150‑1 (MQ=255)
                                                                                                                                                  gTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  2:971955/150‑1 (MQ=255)
                                                                                                                                                  gTACTGCGATTTT‑TCGGCCATATATGACGGTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  1:860991/150‑1 (MQ=255)
                                                                                                                                                     cTGCGATTTT‑TCGGCAATATATGAGGGGGGACCAGGCACATTTTTATAGCGACAAATGTCCACTCCAGCGGAGTAGTCCAATTTGTAGTAAGGGCCGGTAGCGATGGGGAACGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  <  1:1532124/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681820‑2682118

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