Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,969 C→G 6.9% Q259H (CAG→CAC) ‡ yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,681,9690CG 6.9% 329.8 / 12.5 102Q259H (CAG→CAC) ‡yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base C (50/45);  new base G (2/5);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 104 total aligned reads displayed.

CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCG  >  NC_000913/2681820‑2682117
                                                                                                                                                     |                                                                                                                                                     
ccGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTAc                                                                                                                                                       <  1:883384/150‑1 (MQ=255)
ccGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTAc                                                                                                                                                       <  2:1601467/150‑1 (MQ=255)
 cGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAAATCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACt                                                                                                                                                      <  1:473155/150‑1 (MQ=255)
       tGTCGGTCAGTCTCTTATCCCACGCATGGCCAAATTCTCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtt                                                                                                                                                <  1:1174087/150‑1 (MQ=255)
          cGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttg                                                                                                                                             <  2:1662426/150‑2 (MQ=255)
            gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCTCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAATGCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCg                                                                                                                                          <  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)
              cAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGAAAAATCGCAGTACGATTTTGTCGGc                                                                                                                                        <  2:1049670/150‑1 (MQ=255)
                gACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCa                                                                                                                                      >  2:66185/1‑150 (MQ=255)
                    ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCatata                                                                                                                                  >  1:705743/1‑150 (MQ=255)
                      tATCCCACGCCTGGCCAAATTGACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATCTGGTGGGGGGCAACGTGCAGCACACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGAACTTCGATTTTTCGGCCATATATg                                                                                                                                <  2:865448/150‑1 (MQ=255)
                      tATCCCACGCCTGGCCAAATTCACTGAGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATg                                                                                                                                <  2:1399272/150‑1 (MQ=255)
                       aTCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGa                                                                                                                               >  1:425985/1‑150 (MQ=255)
                       aTCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGa                                                                                                                               <  1:165809/150‑1 (MQ=255)
                         cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAgg                                                                                                                             <  1:1282485/150‑1 (MQ=255)
                         cccACGCCTGGCCAAATTCACTGTGACCCCAACTCCACTGTTTTTTACCTTGCGCAATATGGTGGTTGGAAACGTGCAGCAAACCGCCATCTTCATCGTCACACCACCCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAgg                                                                                                                             <  1:280648/150‑1 (MQ=255)
                          ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGt                                                                                                                            >  2:883121/1‑150 (MQ=255)
                            aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTg                                                                                                                          <  1:1665685/150‑1 (MQ=255)
                              gCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGa                                                                                                                        >  1:927817/1‑150 (MQ=255)
                               ccTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGCCAAAATCGTACGGCGATTTTTCGGCCATATATGAGGTTGGaa                                                                                                                       >  1:1469608/1‑150 (MQ=255)
                               ccTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGaa                                                                                                                       >  2:1418189/1‑150 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAgg                                                                                                                   <  1:455571/150‑1 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAg                                                                                                                    >  1:242956/1‑149 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAg                                                                                                                    <  2:242956/149‑1 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTGTTTTACCTGGCGCAATATGGTGGTTGGCACCGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCTTACTGCGATTATTCGGCCATATATGAGGTTGGAACatg                                                                                                                   >  1:859222/1‑148 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAAc                                                                                                                      <  1:66758/146‑1 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAAc                                                                                                                      >  2:66758/1‑146 (MQ=255)
                                         ttCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAtt                                                                                                             >  1:1608672/1‑150 (MQ=255)
                                            aCTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCCCGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAttttt                                                                                                          >  1:535226/1‑150 (MQ=255)
                                             cTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTta                                                                                                         <  1:351566/150‑1 (MQ=255)
                                              tgtgTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTtat                                                                                                        <  2:715291/150‑1 (MQ=255)
                                               gtgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTtata                                                                                                       >  1:556816/1‑150 (MQ=255)
                                                 gtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGc                                                                                                     >  2:13512/1‑150 (MQ=255)
                                                   ccccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                   >  1:1611620/1‑150 (MQ=255)
                                                      cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaa                                                                                                <  2:577664/150‑1 (MQ=255)
                                                      cAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaa                                                                                                <  2:584854/150‑1 (MQ=255)
                                                       aaCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGaaa                                                                                               <  2:1457902/150‑1 (MQ=255)
                                                        aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAAt                                                                                              <  1:790995/150‑1 (MQ=255)
                                                         cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                             <  2:1400295/150‑1 (MQ=255)
                                                         cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                             <  1:1669980/150‑1 (MQ=255)
                                                          tCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTTCAGCAAAGCGCCATCTTGATCGTGACAACACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGt                                                                                            <  1:1034044/150‑1 (MQ=255)
                                                           ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCCTCTTCCTCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCTCATTTTTATAGCGAGAAATGTc                                                                                           >  2:1004418/1‑150 (MQ=255)
                                                           ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTc                                                                                           >  2:840565/1‑150 (MQ=255)
                                                             aCTGTTTTTTACCTGGCGCAATATGGTGGGTGGCAACGTGCAGCAAACCGCCGTCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCa                                                                                         <  2:927110/150‑1 (MQ=255)
                                                                 ttttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTcc                                                                                     >  2:1373912/1‑150 (MQ=255)
                                                                   ttttACGTGGCGCAATATGGTGGTTGCCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAg                                                                                   <  1:1135829/150‑1 (MQ=255)
                                                                   ttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAg                                                                                   <  2:1681742/150‑1 (MQ=255)
                                                                      tACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCAATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgg                                                                                >  1:859654/1‑150 (MQ=255)
                                                                        ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                              >  1:1362196/1‑150 (MQ=255)
                                                                        ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                              <  1:1427661/150‑1 (MQ=255)
                                                                        ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                              <  2:1401565/150‑1 (MQ=255)
                                                                         cTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagt                                                                             >  1:901486/1‑150 (MQ=255)
                                                                           ggCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGagtag                                                                           <  1:78917/150‑1 (MQ=255)
                                                                            gcgcAATATGGGGGCTGGCTTCGTGCAGCAAAGCGCCATCTTCATCGTGACACCACGCGCCGACATAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGATATGTCCACTGCAGCGGagtagt                                                                          <  2:419419/150‑1 (MQ=255)
                                                                                 atatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                     <  1:1049744/150‑1 (MQ=255)
                                                                                 atatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                     <  1:1206385/150‑1 (MQ=255)
                                                                                 atatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTGGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACt                                                                     <  1:614052/150‑1 (MQ=255)
                                                                                     ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATAGTGACACCACGCGCCGACAAAATCGTACTGGGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                 >  1:1125286/1‑150 (MQ=255)
                                                                                        ggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtagt                                                              >  1:134962/1‑150 (MQ=255)
                                                                                             gCAACGTGCCGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTg                                                         >  1:78947/1‑150 (MQ=255)
                                                                                             gCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTg                                                         >  1:504727/1‑150 (MQ=255)
                                                                                                 cGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGAGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCgg                                                     >  1:110670/1‑150 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                 >  2:127258/1‑150 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                 >  2:1656088/1‑150 (MQ=255)
                                                                                                       gcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGa                                               >  1:603195/1‑150 (MQ=255)
                                                                                                            ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                          >  2:1388935/1‑150 (MQ=255)
                                                                                                            ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                          >  2:672769/1‑150 (MQ=255)
                                                                                                              gCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGaa                                        >  2:544262/1‑150 (MQ=255)
                                                                                                               ccATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGaaa                                       <  2:134962/150‑1 (MQ=255)
                                                                                                                cATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAg                                      >  2:56734/1‑150 (MQ=255)
                                                                                                                   cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                   <  1:1201645/150‑1 (MQ=255)
                                                                                                                      cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                >  1:337686/1‑150 (MQ=255)
                                                                                                                         cGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACgg                             >  1:1334303/1‑150 (MQ=255)
                                                                                                                          gTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                            >  2:645735/1‑150 (MQ=255)
                                                                                                                          gTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGGCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                            <  1:883121/150‑1 (MQ=255)
                                                                                                                              caccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                        >  2:1569239/1‑150 (MQ=255)
                                                                                                                              caccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                        >  2:477851/1‑150 (MQ=255)
                                                                                                                                ccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                      >  2:587730/1‑150 (MQ=255)
                                                                                                                                ccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGttt                      >  2:1059802/1‑150 (MQ=255)
                                                                                                                                 cacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTg                     >  2:1118013/1‑150 (MQ=255)
                                                                                                                                   ctcgcCGTAAAATTCGTACTGCGTTTTTTCGGCCATAAATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                   <  2:864081/148‑1 (MQ=255)
                                                                                                                                   cgcgCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                   <  1:47087/150‑1 (MQ=255)
                                                                                                                                   cgcgCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGcc                   <  1:243390/150‑1 (MQ=255)
                                                                                                                                       gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                >  1:1521724/6‑150 (MQ=255)
                                                                                                                                       ccGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGa               >  2:255117/1‑150 (MQ=255)
                                                                                                                                        ccgaaaaaTCGCAGTACGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt              >  1:448210/5‑150 (MQ=255)
                                                                                                                                        cGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt              <  2:425985/150‑1 (MQ=255)
                                                                                                                                        cGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt              <  2:676285/150‑1 (MQ=255)
                                                                                                                                         gACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc             >  2:418643/1‑150 (MQ=255)
                                                                                                                                         gACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc             >  2:17828/1‑150 (MQ=255)
                                                                                                                                         gACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATc             <  2:110670/150‑1 (MQ=255)
                                                                                                                                          aCAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCa            <  1:1134566/150‑1 (MQ=255)
                                                                                                                                             aaaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGATAAATGTCCACTCCAGAGGATTAGTCCACTTTGTAGTAAGTGCCGTTAGCGATGTGGAAAGCGGAGACGGCCCGTTTTCCGTGAtacacc         >  1:1449863/1‑145 (MQ=255)
                                                                                                                                              aaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAaca        <  1:13512/150‑1 (MQ=255)
                                                                                                                                               aTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAacac       <  2:51464/150‑1 (MQ=255)
                                                                                                                                               aTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAacac       >  2:195117/1‑150 (MQ=255)
                                                                                                                                                    aCTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAACCAccgccg  >  2:813648/1‑150 (MQ=255)
                                                                                                                                                    aCTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAccgccg  >  2:651072/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCG  >  NC_000913/2681820‑2682117

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: