Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,971 G→A 5.0% 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.0% 320.7 / 3.2 99Q259* (CAG→TAG) ‡yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (40/54);  new base A (3/2);  total (43/56)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

GTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681822‑2682118
                                                                                                                                                     |                                                                                                                                                    
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)
        cGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAttttt                                                                                                                                             >  2:1521984/1‑150 (MQ=255)
        cGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttg                                                                                                                                              <  1:1648921/150‑2 (MQ=255)
          gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGACATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCg                                                                                                                                           >  2:1118232/1‑150 (MQ=255)
           tCAGACTCTTATCCCACGCCTGGCCTAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCgg                                                                                                                                          <  1:266200/150‑1 (MQ=255)
             aGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGcc                                                                                                                                        >  2:141749/1‑150 (MQ=255)
               aCTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCat                                                                                                                                      <  1:642107/150‑1 (MQ=255)
                  ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTGTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCatata                                                                                                                                   <  1:1702961/150‑1 (MQ=255)
                    tATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATg                                                                                                                                 <  1:885551/150‑1 (MQ=255)
                    tATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATg                                                                                                                                 >  1:927240/1‑150 (MQ=255)
                        ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCg                                                                                                                                           <  1:758319/136‑1 (MQ=255)
                        ccACGCCTGGCCAAAATCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCg                                                                                                                                           >  2:758319/1‑136 (MQ=255)
                          aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTg                                                                                                                           >  2:1646368/1‑150 (MQ=255)
                            gCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATAGTACTGCGATTTTTCGGCCATATATGAGGTTGGa                                                                                                                         <  1:1549698/150‑1 (MQ=255)
                             ccTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGaa                                                                                                                        <  2:1247712/150‑1 (MQ=255)
                              cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAAc                                                                                                                       >  1:1387327/1‑150 (MQ=255)
                                ggCCAAATTCACTGTGTACCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAg                                                                                                                     >  1:1081345/1‑150 (MQ=255)
                                 gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAgg                                                                                                                    <  1:141749/150‑1 (MQ=255)
                                  ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   <  2:1072853/150‑1 (MQ=255)
                                  ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   >  1:324737/1‑150 (MQ=255)
                                  ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   <  1:204425/150‑1 (MQ=255)
                                    aaaTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGcac                                                                                                                 >  1:1382272/1‑150 (MQ=255)
                                      aTTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAt                                                                                                               >  1:665575/1‑150 (MQ=255)
                                          aCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAttttt                                                                                                           <  1:556781/150‑1 (MQ=255)
                                            tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  2:1296516/150‑1 (MQ=255)
                                            tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  1:1296837/150‑1 (MQ=255)
                                            tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  1:1349627/150‑1 (MQ=255)
                                             gtgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGGGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTtata                                                                                                        <  1:434919/150‑1 (MQ=255)
                                                tCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCg                                                                                                     >  1:303970/1‑150 (MQ=255)
                                                tCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCg                                                                                                     >  1:617118/1‑150 (MQ=255)
                                                 ccccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    <  2:1580380/150‑1 (MQ=255)
                                                  cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  2:1444292/150‑1 (MQ=255)
                                                    cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAATCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaa                                                                                                 >  2:1117016/1‑150 (MQ=255)
                                                    cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCTACAAAATCGTAATGCGATTTTTCGGCAATTTATGAGGTTTGAACAGGCACATTTTTATAGCGAGaa                                                                                                 >  2:1532124/1‑150 (MQ=255)
                                                     acctaaaCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCTATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGATaaa                                                                                                <  2:347314/144‑1 (MQ=255)
                                                      aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAAt                                                                                               <  2:1382272/150‑1 (MQ=255)
                                                         ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTc                                                                                            >  2:43214/1‑150 (MQ=255)
                                                          cACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTcc                                                                                           >  2:1416112/1‑150 (MQ=255)
                                                               ttttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTcc                                                                                      <  1:929605/150‑1 (MQ=255)
                                                                tttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCa                                                                                     >  1:1479475/1‑150 (MQ=255)
                                                                    tACCTGGCGCAATATGGTGGGTGGCAACGTGCAGCCAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgg                                                                                 <  1:225913/150‑1 (MQ=255)
                                                                     aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                                >  1:971955/1‑150 (MQ=255)
                                                                     aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                                >  1:978552/1‑150 (MQ=255)
                                                                       cTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagt                                                                              >  1:1253729/1‑150 (MQ=255)
                                                                           cgcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTc                                                                          <  2:556848/150‑1 (MQ=255)
                                                                               atatGGGGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                      <  2:324737/150‑1 (MQ=255)
                                                                                  tggtggTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTg                                                                   <  2:1253729/150‑1 (MQ=255)
                                                                                   ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  1:985708/150‑1 (MQ=255)
                                                                                   ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  2:1632667/150‑1 (MQ=255)
                                                                                     tggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtag                                                                <  1:1434751/150‑1 (MQ=255)
                                                                                       ggtGGCCACGTGCAGCAAACCGCCATCTTCATCGTGACACCCCGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtagta                                                              <  1:698900/148‑1 (MQ=255)
                                                                                           gCAACTTGCAGCAAACCGCAATCTTCATCCTGACACCCCGCGCCGACAAAATCCTACTGCGATTTTTCGGCCATATATGATGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCAACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTg                                                          <  2:508664/150‑1 (MQ=255)
                                                                                               cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                      <  2:733187/150‑1 (MQ=255)
                                                                                               cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                      <  2:1308550/150‑1 (MQ=255)
                                                                                               cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                      <  2:1148008/150‑1 (MQ=255)
                                                                                                 tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGTCCATATATGAGGTTGGAACAGGCACAATTTAATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGTTa                                                    >  2:919324/1‑150 (MQ=255)
                                                                                                 tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTa                                                    >  2:515173/1‑150 (MQ=255)
                                                                                                 tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTa                                                    >  2:205885/1‑150 (MQ=255)
                                                                                                   cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     >  1:1258358/1‑131 (MQ=255)
                                                                                                   cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     <  2:1258358/131‑1 (MQ=255)
                                                                                                       aaaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATg                                              <  1:43214/150‑1 (MQ=255)
                                                                                                       aaaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATg                                              >  1:763731/1‑150 (MQ=255)
                                                                                                         aCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATggg                                            <  2:665575/150‑1 (MQ=255)
                                                                                                                tCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCg                                     >  1:758369/1‑150 (MQ=255)
                                                                                                                 cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    <  1:1117016/150‑1 (MQ=255)
                                                                                                                 cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    <  1:5275/150‑1 (MQ=255)
                                                                                                                  ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  2:1400381/150‑1 (MQ=255)
                                                                                                                  ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  1:1405681/150‑1 (MQ=255)
                                                                                                                  ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGGCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  1:453561/150‑1 (MQ=255)
                                                                                                                    cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                 >  2:224320/1‑150 (MQ=255)
                                                                                                                        gTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                             >  2:477739/1‑150 (MQ=255)
                                                                                                                          gACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGccc                           <  1:977538/150‑1 (MQ=255)
                                                                                                                          gACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGccc                           >  2:226174/1‑150 (MQ=255)
                                                                                                                             accaaGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTT‑‑TAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCAGTTTg                      <  1:1118232/150‑1 (MQ=255)
                                                                                                                              ccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                       >  1:721928/1‑150 (MQ=255)
                                                                                                                              ccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                       >  1:1653718/1‑150 (MQ=255)
                                                                                                                                 cgcgCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                    >  1:922723/1‑150 (MQ=255)
                                                                                                                                  gcgcCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCg                   >  1:775811/1‑150 (MQ=255)
                                                                                                                                  gcgcCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCg                   <  2:927240/150‑1 (MQ=255)
                                                                                                                                   cgcCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGt                  <  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)
                                                                                                                                     ccGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGa                >  1:1469116/1‑150 (MQ=255)
                                                                                                                                      cGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt               >  1:247755/1‑150 (MQ=255)
                                                                                                                                       cgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt               <  2:1389524/147‑1 (MQ=255)
                                                                                                                                         cAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                           >  2:52864/1‑119 (MQ=255)
                                                                                                                                         cAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                           <  1:52864/119‑1 (MQ=255)
                                                                                                                                          aaaaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaaa           <  1:552089/150‑1 (MQ=255)
                                                                                                                                            aaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca         <  1:1521984/150‑1 (MQ=255)
                                                                                                                                            aaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca         >  2:966359/1‑150 (MQ=255)
                                                                                                                                             aTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAacac        <  1:1646368/150‑1 (MQ=255)
                                                                                                                                                gTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  2:1387327/150‑1 (MQ=255)
                                                                                                                                                gTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  2:971955/150‑1 (MQ=255)
                                                                                                                                                gTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  2:978552/150‑1 (MQ=255)
                                                                                                                                                gTACTGCGATTTTTCGGCCATATATGACGGTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  1:860991/150‑1 (MQ=255)
                                                                                                                                                   cTGCGATTTTTCGGCAATATATGAGGGGGGACCAGGCACATTTTTATAGCGACAAATGTCCACTCCAGCGGAGTAGTCCAATTTGTAGTAAGGGCCGGTAGCGATGGGGAACGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  <  1:1532124/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681822‑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: