Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,590,876 C→T 6.1% W12* (TGG→TGA)  ugpA ← sn‑glycerol‑3‑phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,590,8760CT 6.1% 215.9 / 3.7 66W12* (TGG→TGA) ugpAsn‑glycerol‑3‑phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base T (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CGTGACGAAGTTATCCAGGCCGACAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGAT  >  NC_000913/3590729‑3591024
                                                                                                                                                   |                                                                                                                                                    
cGTGACGAAGTTATCCAGGCCGACAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGcagcc                                                                                                                                                     >  1:22562/1‑149 (MQ=255)
 gTGACGAAGTTATCCAGGCCGACAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGcagcca                                                                                                                                                    >  2:526190/1‑149 (MQ=255)
            aTCCAGGCCGACAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGa                                                                                                                                         >  2:314444/1‑149 (MQ=255)
               cAGGCCGACAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAaca                                                                                                                                      <  3:189088/149‑1 (MQ=255)
                      aCAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACg                                                                                                                               <  1:246850/149‑1 (MQ=255)
                      aCAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACg                                                                                                                               >  2:246850/1‑149 (MQ=255)
                         aaCTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGga                                                                                                                            <  2:309862/149‑1 (MQ=255)
                            tggctggAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGgatga                                                                                                                         <  1:79159/149‑1 (MQ=255)
                             ggctggAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGga                                                                                                                            >  2:629773/1‑145 (MQ=255)
                             ggctggAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGga                                                                                                                            <  1:629773/145‑1 (MQ=255)
                              gctggAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGatgatg                                                                                                                       <  4:654179/149‑1 (MQ=255)
                                    gaAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTa                                                                                                                 <  4:308482/149‑1 (MQ=255)
                                       aCCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACgg                                                                                                              >  4:340123/1‑149 (MQ=255)
                                          aaaCGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGTCAGCGCGAGc                                                                                                                                            <  3:263482/116‑1 (MQ=255)
                                          aaaCGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGTCAGCGCGAGc                                                                                                                                            >  4:263482/1‑116 (MQ=255)
                                           aaCGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTaa                                                                                                          <  1:424385/149‑1 (MQ=255)
                                           aaCGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTaa                                                                                                          <  1:526190/149‑1 (MQ=255)
                                           aaCGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTaa                                                                                                          <  1:529489/149‑1 (MQ=255)
                                              ggATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTc                                                                                                       <  4:298811/149‑1 (MQ=255)
                                                 tCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTg                                                                                                    <  4:42685/149‑1 (MQ=255)
                                                   gACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGaa                                                                                                  >  4:221368/1‑149 (MQ=255)
                                                   gACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGaa                                                                                                  <  1:314444/149‑1 (MQ=255)
                                                   gACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGaa                                                                                                  <  3:132812/149‑1 (MQ=255)
                                                    aCGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTAc                                                                                                                <  1:528906/134‑1 (MQ=255)
                                                    aCGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTAc                                                                                                                >  2:528906/1‑134 (MQ=255)
                                                    aCGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAAc                                                                                                 >  4:377122/1‑149 (MQ=255)
                                                      gCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACtg                                                                                               >  1:38694/1‑149 (MQ=255)
                                                      gCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACtg                                                                                               >  4:572135/1‑149 (MQ=255)
                                                        tttGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACtgtg                                                                                             >  2:379929/1‑149 (MQ=255)
                                                             aGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGa                                                                                        <  3:377122/149‑1 (MQ=255)
                                                               cGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATg                                                                                      <  2:22562/149‑1 (MQ=255)
                                                                 aGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATgcg                                                                                    >  3:304137/1‑149 (MQ=255)
                                                                 aGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATgcg                                                                                    >  1:20337/1‑149 (MQ=255)
                                                                     ccACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGcgcttc                                                                                <  1:654584/149‑1 (MQ=255)
                                                                      cacaACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGTATGcgcttcg                                                                               >  4:440562/1‑149 (MQ=255)
                                                                      cacaACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGcgcttcg                                                                               <  4:355453/149‑1 (MQ=255)
                                                                       acaACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGTCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGcgcttcgc                                                                              >  3:428315/1‑149 (MQ=255)
                                                                       acaACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGcgcttcgc                                                                              >  4:26381/1‑149 (MQ=255)
                                                                       acaACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGcgcttcgc                                                                              >  4:168516/1‑149 (MQ=255)
                                                                        caACGCTGCGCCCGCAGGCCAGATAAAAAAGATAACGGTGAGGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGcgcttcgct                                                                             <  2:468973/149‑1 (MQ=255)
                                                                           cGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTAt                                                                          <  3:668695/149‑1 (MQ=255)
                                                                               tCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCgg                                                                      <  4:304137/149‑1 (MQ=255)
                                                                                   ccGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTa                                                                  <  2:326290/149‑1 (MQ=255)
                                                                                     gCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTaca                                                                <  4:377552/149‑1 (MQ=255)
                                                                                         gCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACAcgcg                                                            >  3:193821/1‑149 (MQ=255)
                                                                                          ccAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGa                                                           >  3:384883/1‑149 (MQ=255)
                                                                                           cAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGAt                                                          >  2:126972/1‑149 (MQ=255)
                                                                                            aGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATg                                                         >  2:247858/1‑149 (MQ=255)
                                                                                                  aaaaGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGt                                                   <  3:26381/149‑1 (MQ=255)
                                                                                                  aaaaGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGt                                                   <  4:254758/149‑1 (MQ=255)
                                                                                                   aaaGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGa                                                                                        <  4:118232/111‑1 (MQ=255)
                                                                                                   aaaGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGa                                                                                        >  3:118232/1‑111 (MQ=255)
                                                                                                     aGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAg                                                <  4:193821/149‑1 (MQ=255)
                                                                                                      gATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGa                                                                                        >  2:84990/1‑108 (MQ=255)
                                                                                                      gATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGa                                                                                        <  1:84990/108‑1 (MQ=255)
                                                                                                         aaCGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATaa                                            <  1:247858/149‑1 (MQ=255)
                                                                                                              tgatgatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGc                                          <  1:313237/146‑1 (MQ=255)
                                                                                                              tgatgatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGc                                          >  2:313237/1‑146 (MQ=255)
                                                                                                                  gatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGc                                   >  3:28797/1‑149 (MQ=255)
                                                                                                                     gaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATc                                <  4:369644/149‑1 (MQ=255)
                                                                                                                             gcgcGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACa                           >  1:247920/1‑146 (MQ=255)
                                                                                                                             gcgcGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACa                           <  2:247920/146‑1 (MQ=255)
                                                                                                                                  aCCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGACGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTg                   <  3:440562/149‑1 (MQ=255)
                                                                                                                                    cagcagATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGAt                 >  4:395297/1‑149 (MQ=255)
                                                                                                                                    cagcagATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGAt                 >  2:204331/1‑149 (MQ=255)
                                                                                                                                    cagcagATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGAt                 <  2:20337/149‑1 (MQ=255)
                                                                                                                                             aGGCAGTCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTtc        <  4:428315/149‑1 (MQ=255)
                                                                                                                                             aGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTtc        >  1:160085/1‑149 (MQ=255)
                                                                                                                                                cagccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTtcgtc     <  3:340123/149‑1 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATaa                                            <  1:368614/107‑1 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATaa                                            <  1:444627/107‑1 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATaa                                            >  2:444627/1‑107 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATaa                                            >  2:368614/1‑107 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTaa              >  3:608877/1‑137 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTaa              <  4:608877/137‑1 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGAt  >  1:325933/1‑149 (MQ=255)
                                                                                                                                                   ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGAt  <  3:168516/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTGACGAAGTTATCCAGGCCGACAAACTGGCTGGAGAAACCAAACGGATCGACGCTTTGTAGCGAGTACCACAACGCTTCGCCCGCAGGCCAGATAAAAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGAT  >  NC_000913/3590729‑3591024

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑