Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,971 G→A 6.9% 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 6.9% 315.0 / 7.6 102Q259* (CAG→TAG) ‡yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (51/43);  new base A (2/5);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 102 total aligned reads displayed.

TGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTT  >  NC_000913/2681827‑2682119
                                                                                                                                                |                                                                                                                                                     
tGTCGGTCAGTCTCTTATCCCACGCATGGCCAAATTCTCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtt                                                                                                                                                  <  1:1174087/150‑1 (MQ=255)
   cGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttg                                                                                                                                               <  2:1662426/150‑2 (MQ=255)
     gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCTCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAATGCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCg                                                                                                                                            <  1:746951/150‑1 (MQ=255)
     gTCAGACTCTTATACCACGCCTGGCCAAATTCACTGTGTCCACAACTCCACTGTTTTTTACCTGGCGCAAAATGGTGGTTGGCAACGTGCAGCAAACCTCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttgtc                                                                                                                                             <  2:335947/150‑4 (MQ=255)
      tCAGACTCTTATCCCACGCGTGGACTTTTTCTCTGTGTGCCCAAGTCCACTGGTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttgtcg                                                                                                                                            <  2:1068096/150‑5 (MQ=255)
      gcAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACATGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttgtcg                                                                                                                                            <  2:1143770/149‑5 (MQ=255)
       cAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGA‑AAAATCGCAGTACGATTTTGTCGGc                                                                                                                                          <  2:1049670/150‑1 (MQ=255)
         gACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCa                                                                                                                                        >  2:66185/1‑150 (MQ=255)
             ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCatata                                                                                                                                    >  1:705743/1‑150 (MQ=255)
               tATCCCACGCCTGGCCAAATTGACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATCTGGTGGGGGGCAACGTGCAGCACACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGAACTTCGATTTT‑TCGGCCATATATg                                                                                                                                  <  2:865448/150‑1 (MQ=255)
               tATCCCACGCCTGGCCAAATTCACTGAGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATg                                                                                                                                  <  2:1399272/150‑1 (MQ=255)
                aTCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGa                                                                                                                                 <  1:165809/150‑1 (MQ=255)
                aTCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGa                                                                                                                                 >  1:425985/1‑150 (MQ=255)
                  cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAgg                                                                                                                               <  1:1282485/150‑1 (MQ=255)
                  cccACGCCTGGCCAAATTCACTGTGACCCCAACTCCACTGTTTTTTACCTTGCGCAATATGGTGGTTGGAAACGTGCAGCAAACCGCCATCTTCATCGTCACACCACCCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAgg                                                                                                                               <  1:280648/150‑1 (MQ=255)
                   ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGt                                                                                                                              >  2:883121/1‑150 (MQ=255)
                     aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                            <  1:1665685/150‑1 (MQ=255)
                       gCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGa                                                                                                                          >  1:927817/1‑150 (MQ=255)
                        ccTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGCCAAAATCGTACGGCGATTTT‑TCGGCCATATATGAGGTTGGaa                                                                                                                         >  1:1469608/1‑150 (MQ=255)
                        ccTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGaa                                                                                                                         >  2:1418189/1‑150 (MQ=255)
                            gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAgg                                                                                                                     <  1:455571/150‑1 (MQ=255)
                            gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAg                                                                                                                      <  2:242956/149‑1 (MQ=255)
                            gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAg                                                                                                                      >  1:242956/1‑149 (MQ=255)
                            gCCAAATTCACTGTGTCCCCAACTCCACTGTGTTTTACCTGGCGCAATATGGTGGTTGGCACCGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCTTACTGCGATTAT‑TCGGCCATATATGAGGTTGGAACatg                                                                                                                     >  1:859222/1‑148 (MQ=255)
                            gCCAAATTCACTGTGTCCCCAACTCCACTG‑TTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAAc                                                                                                                        >  2:66758/1‑146 (MQ=255)
                            gCCAAATTCACTGTGTCCCCAACTCCACTG‑TTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAAc                                                                                                                        <  1:66758/146‑1 (MQ=255)
                                  ttCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACAtt                                                                                                               >  1:1608672/1‑150 (MQ=255)
                                     aCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCCCGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACAttttt                                                                                                            >  1:535226/1‑150 (MQ=255)
                                      cTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTta                                                                                                           <  1:351566/150‑1 (MQ=255)
                                       tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                          <  2:715291/150‑1 (MQ=255)
                                        gtgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTtata                                                                                                         >  1:556816/1‑150 (MQ=255)
                                          gtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGc                                                                                                       >  2:13512/1‑150 (MQ=255)
                                            ccccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                     >  1:1611620/1‑150 (MQ=255)
                                               cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaa                                                                                                  <  2:577664/150‑1 (MQ=255)
                                               cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaa                                                                                                  <  2:584854/150‑1 (MQ=255)
                                                aaCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaaa                                                                                                 <  2:1457902/150‑1 (MQ=255)
                                                 aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAAt                                                                                                <  1:790995/150‑1 (MQ=255)
                                                  cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                               <  2:1400295/150‑1 (MQ=255)
                                                  cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                               <  1:1669980/150‑1 (MQ=255)
                                                   tCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTTCAGCAAAGCGCCATCTTGATCGTGACAACACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGt                                                                                              <  1:1034044/150‑1 (MQ=255)
                                                    ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCCTCTTCCTCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCTCATTTTTATAGCGAGAAATGTc                                                                                             >  2:1004418/1‑150 (MQ=255)
                                                    ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTc                                                                                             >  2:840565/1‑150 (MQ=255)
                                                      aCTGTTTTTTACCTGGCGCAATATGGTGGGTGGCAACGTGCAGCAAACCGCCGTCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCa                                                                                           <  2:927110/150‑1 (MQ=255)
                                                          ttttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTcc                                                                                       >  2:1373912/1‑150 (MQ=255)
                                                            ttttACGTGGCGCAATATGGTGGTTGCCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAg                                                                                     <  1:1135829/150‑1 (MQ=255)
                                                            ttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAg                                                                                     <  2:1681742/150‑1 (MQ=255)
                                                               tACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCAATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgg                                                                                  >  1:859654/1‑150 (MQ=255)
                                                                 ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                                >  1:1362196/1‑150 (MQ=255)
                                                                 ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                                <  1:1427661/150‑1 (MQ=255)
                                                                 ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                                <  2:1401565/150‑1 (MQ=255)
                                                                  cTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagt                                                                               >  1:901486/1‑150 (MQ=255)
                                                                    ggCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagtag                                                                             <  1:78917/150‑1 (MQ=255)
                                                                     gcgcAATATGGGGGCTGGCTTCGTGCAGCAAAGCGCCATCTTCATCGTGACACCACGCGCCGACATAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGATATGTCCACTGCAGCGGagtagt                                                                            <  2:419419/150‑1 (MQ=255)
                                                                      cgcAATATGGGGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGACGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTc                                                                           <  1:25524/150‑1 (MQ=255)
                                                                          atatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                       <  1:1049744/150‑1 (MQ=255)
                                                                          atatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                       <  1:1206385/150‑1 (MQ=255)
                                                                          atatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTGGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                       <  1:614052/150‑1 (MQ=255)
                                                                              ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATAGTGACACCACGCGCCGACAAAATCGTACTGGGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                   >  1:1125286/1‑150 (MQ=255)
                                                                                 ggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtagt                                                                >  1:134962/1‑150 (MQ=255)
                                                                                      gCAACGTGCCGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTg                                                           >  1:78947/1‑150 (MQ=255)
                                                                                      gCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTg                                                           >  1:504727/1‑150 (MQ=255)
                                                                                       cAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGc                                                          >  2:243390/1‑150 (MQ=255)
                                                                                          cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGAGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                       >  1:110670/1‑150 (MQ=255)
                                                                                              cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                   >  2:127258/1‑150 (MQ=255)
                                                                                              cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                   >  2:1656088/1‑150 (MQ=255)
                                                                                                gcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGa                                                 >  1:603195/1‑150 (MQ=255)
                                                                                                     ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                            >  2:1388935/1‑150 (MQ=255)
                                                                                                     ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                            >  2:672769/1‑150 (MQ=255)
                                                                                                       gCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGaa                                          >  2:544262/1‑150 (MQ=255)
                                                                                                        ccATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGaaa                                         <  2:134962/150‑1 (MQ=255)
                                                                                                         cATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAg                                        >  2:56734/1‑150 (MQ=255)
                                                                                                            cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                     <  1:1201645/150‑1 (MQ=255)
                                                                                                               cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                  >  1:337686/1‑150 (MQ=255)
                                                                                                                  cGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACgg                               >  1:1334303/1‑150 (MQ=255)
                                                                                                                   gTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                              >  2:645735/1‑150 (MQ=255)
                                                                                                                   gTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGGCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                              <  1:883121/150‑1 (MQ=255)
                                                                                                                       caccacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                          >  2:1569239/1‑150 (MQ=255)
                                                                                                                       caccacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                          >  2:477851/1‑150 (MQ=255)
                                                                                                                         ccacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                        >  2:587730/1‑150 (MQ=255)
                                                                                                                         ccacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                        >  2:1059802/1‑150 (MQ=255)
                                                                                                                          cacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTg                       >  2:1118013/1‑150 (MQ=255)
                                                                                                                            ctcgcCGTAAAATTCGTACTGCGTTTTT‑TCGGCCATAAATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                     <  2:864081/148‑1 (MQ=255)
                                                                                                                            cgcgCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                     <  1:243390/150‑1 (MQ=255)
                                                                                                                            cgcgCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                     <  1:47087/150‑1 (MQ=255)
                                                                                                                                gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                  >  1:1521724/6‑150 (MQ=255)
                                                                                                                                ccGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGa                 >  2:255117/1‑150 (MQ=255)
                                                                                                                                 ccgaaaaaTCGCAGTACGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt                >  1:448210/5‑150 (MQ=255)
                                                                                                                                 cGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt                <  2:425985/150‑1 (MQ=255)
                                                                                                                                 cGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt                <  2:676285/150‑1 (MQ=255)
                                                                                                                                  gACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc               >  2:418643/1‑150 (MQ=255)
                                                                                                                                  gACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc               <  2:110670/150‑1 (MQ=255)
                                                                                                                                  gACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc               >  2:17828/1‑150 (MQ=255)
                                                                                                                                   aCAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCa              <  1:1134566/150‑1 (MQ=255)
                                                                                                                                      aaaTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGATAAATGTCCACTCCAGAGGATTAGTCCACTTTGTAGTAAGTGCCGTTAGCGATGTGGAAAGCGGAGACGGCCCGTTTTCCGTGAtacacc           >  1:1449863/1‑145 (MQ=255)
                                                                                                                                       aaTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca          <  1:13512/150‑1 (MQ=255)
                                                                                                                                        aTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAacac         <  2:51464/150‑1 (MQ=255)
                                                                                                                                        aTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAacac         >  2:195117/1‑150 (MQ=255)
                                                                                                                                             aCTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAACCAccgccg    >  2:813648/1‑150 (MQ=255)
                                                                                                                                             aCTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgccg    >  2:651072/1‑150 (MQ=255)
                                                                                                                                               tGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGtt  >  1:476448/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTT  >  NC_000913/2681827‑2682119

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