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

GTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  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)
        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:324737/1‑150 (MQ=255)
                                  ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   <  1:204425/150‑1 (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                                                                                                         <  2:1296516/150‑1 (MQ=255)
                                            tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  1:1296837/150‑1 (MQ=255)
                                            tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                         <  1:1349627/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:971955/1‑150 (MQ=255)
                                                                     aCCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGa                                                                                >  1:978552/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                                                                  <  1:985708/150‑1 (MQ=255)
                                                                                   ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  2:1632667/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:1308550/150‑1 (MQ=255)
                                                                                               cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                      <  2:1148008/150‑1 (MQ=255)
                                                                                                 tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGTCCATATATGAGGTTGGAACAGGCACAATTTAATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGTTa                                                    >  2:919324/1‑150 (MQ=255)
                                                                                                 tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTa                                                    >  2:515173/1‑150 (MQ=255)
                                                                                                 tGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTa                                                    >  2:205885/1‑150 (MQ=255)
                                                                                                   cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     >  1:1258358/1‑131 (MQ=255)
                                                                                                   cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     <  2:1258358/131‑1 (MQ=255)
                                                                                                       aaaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATg                                              <  1:43214/150‑1 (MQ=255)
                                                                                                       aaaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATg                                              >  1:763731/1‑150 (MQ=255)
                                                                                                         aCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATggg                                            <  2:665575/150‑1 (MQ=255)
                                                                                                                tCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCg                                     >  1:758369/1‑150 (MQ=255)
                                                                                                                 cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    <  1:1117016/150‑1 (MQ=255)
                                                                                                                 cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    <  1:5275/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                           <  1:977538/150‑1 (MQ=255)
                                                                                                                          gACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGccc                           >  2:226174/1‑150 (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‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  2:978552/150‑1 (MQ=255)
                                                                                                                                                gTACTGCGATTTT‑TCGGCCATATATGACGGTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgc     <  1:860991/150‑1 (MQ=255)
                                                                                                                                                   cTGCGATTTT‑TCGGCAATATATGAGGGGGGACCAGGCACATTTTTATAGCGACAAATGTCCACTCCAGCGGAGTAGTCCAATTTGTAGTAAGGGCCGGTAGCGATGGGGAACGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  <  1:1532124/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  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: ‑