Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,248,399 T→A 5.5% Q45L (CAG→CTG)  lysP ← lysine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,248,3990TA 5.5% 221.7 / 3.4 73Q45L (CAG→CTG) lysPlysine transporter
Reads supporting (aligned to +/- strand):  ref base T (40/29);  new base A (2/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

AAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATC  >  NC_000913/2248252‑2248546
                                                                                                                                                   |                                                                                                                                                   
aaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGACCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGCcctgcct                                                                                                                                                     >  6:120920/1‑148 (MQ=255)
     ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGaaa                                                                                                                                               >  7:179556/1‑149 (MQ=255)
        ttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCg                                                                                                                                            >  5:162480/1‑149 (MQ=255)
         tcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGt                                                                                                                                           >  7:315645/1‑149 (MQ=255)
           aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg                                                                                                                                         >  8:257146/1‑149 (MQ=255)
           aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg                                                                                                                                         <  2:113859/149‑1 (MQ=255)
               tAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc                                                                                                                                      >  4:158760/1‑148 (MQ=255)
                aGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACATGCCTGAGAAATCGTTgcgc                                                                                                                                      >  3:279821/1‑147 (MQ=255)
                 gTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCag                                                                                                                                   <  5:57015/149‑1 (MQ=255)
                 gTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCa                                                                                                                                    <  7:79584/148‑1 (MQ=255)
                          gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTgag                                                                                                                                                  >  2:42437/1‑125 (MQ=255)
                          gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTgag                                                                                                                                                  <  1:42437/125‑1 (MQ=255)
                              gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAg                                                                                                                      <  4:345421/149‑1 (MQ=255)
                                 gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                   <  6:252929/149‑1 (MQ=255)
                                  cAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCt                                                                                                                  <  5:329706/149‑1 (MQ=255)
                                   aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                 >  2:338863/1‑149 (MQ=255)
                                   aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                 >  7:64761/1‑149 (MQ=255)
                                   aaaCGAACAGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCACGAAGTACACCATCACGCCAATCAGCATATACGAGAGCAATGCCCCGCCCTGACCTGCCTGAGAAATCGTTGCCCCAGAGGCAACAAAAAGACCTg                                                                                                                 >  6:320016/1‑149 (MQ=255)
                                     aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTa                                                                                                               <  3:110360/149‑1 (MQ=255)
                                                  cGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGc                                                                                                  >  5:196156/1‑149 (MQ=255)
                                                     cATCTATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCaa                                                                                               <  8:307722/149‑1 (MQ=255)
                                                     cATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCaa                                                                                               <  1:338863/149‑1 (MQ=255)
                                                         tatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGGACCGATGGAACCGCCCATgg                                                                                            >  5:333696/1‑148 (MQ=255)
                                                          atGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGc                                                                                           >  3:73441/1‑148 (MQ=255)
                                                                ccAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAatcatc                                                                                    <  8:64761/149‑1 (MQ=255)
                                                                       aCCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGt                                                                             >  3:185074/1‑149 (MQ=255)
                                                                       aCCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGt                                                                             >  8:330374/1‑149 (MQ=255)
                                                                       aCCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGt                                                                             >  8:352998/1‑149 (MQ=255)
                                                                         cGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGa                                                                           <  1:209223/149‑1 (MQ=255)
                                                                               tGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAcc                                                                                                    >  2:57462/1‑118 (MQ=255)
                                                                               tGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAcc                                                                                                    <  1:57462/118‑1 (MQ=255)
                                                                               tGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGcc                                                                     >  2:39627/1‑149 (MQ=255)
                                                                                 gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCtt                                                                   >  1:109920/1‑149 (MQ=255)
                                                                                  tcatcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                        >  6:248231/1‑143 (MQ=255)
                                                                                  tcatcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                        <  5:248231/143‑1 (MQ=255)
                                                                                    atcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGCGCAATGCCCCGCCCGGACCTGCCTGAGAAATCTTTGCGCCAGACGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCCGCAGGTGACGCGCCTTTaa                                                                >  5:306000/1‑149 (MQ=255)
                                                                                     tcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                        >  2:306644/1‑140 (MQ=255)
                                                                                     tcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                        <  1:306644/140‑1 (MQ=255)
                                                                                     tcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAt                                                               >  8:43651/1‑149 (MQ=255)
                                                                                      caGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                              >  1:369149/1‑149 (MQ=255)
                                                                                        ggAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCaa                                                                                               >  5:183990/1‑114 (MQ=255)
                                                                                        ggAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCaa                                                                                               <  6:183990/114‑1 (MQ=255)
                                                                                        ggAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCa                                                            >  2:315017/1‑149 (MQ=255)
                                                                                          aaGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACg                                                          <  7:330374/149‑1 (MQ=255)
                                                                                          aaGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACg                                                          <  2:369149/149‑1 (MQ=255)
                                                                                           aGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                         <  3:158760/149‑1 (MQ=255)
                                                                                            gTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTa                                                                                                               <  8:140443/94‑1 (MQ=255)
                                                                                              aaaCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCCCGGCg                                                       >  7:211108/1‑148 (MQ=255)
                                                                                              aaaCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCg                                                       >  2:322440/1‑148 (MQ=255)
                                                                                                aCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCAGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAc                                                                                                                    >  6:292880/1‑85 (MQ=255)
                                                                                                aCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCAGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAc                                                                                                                    >  6:287955/1‑85 (MQ=255)
                                                                                                aCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCAGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAc                                                                                                                    <  5:292880/85‑1 (MQ=255)
                                                                                                aCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCAGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAc                                                                                                                    <  5:287955/85‑1 (MQ=255)
                                                                                                        gCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgcg                                            >  5:259969/1‑149 (MQ=255)
                                                                                                          cAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCt                                          >  8:51440/1‑149 (MQ=255)
                                                                                                           aaTCAGCATATACGAGAGCAATGACCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGGCAGGTTACGCGCAGTTAATTCCCGGCGTAAGCCCGGCTCtt                                         >  5:159004/1‑149 (MQ=255)
                                                                                                            aTCAGCAGATACGAGAGCAATGCCCCGCCCGGACCTGCCTGCGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                        <  6:79485/149‑1 (MQ=255)
                                                                                                                   tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                 >  8:231482/1‑149 (MQ=255)
                                                                                                                   tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                 >  1:308127/1‑149 (MQ=255)
                                                                                                                   tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTCCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                 >  7:291279/1‑149 (MQ=255)
                                                                                                                     taCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                     >  5:246827/1‑79 (MQ=255)
                                                                                                                     taCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                     <  6:246827/79‑1 (MQ=255)
                                                                                                                     taCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                <  5:156527/148‑1 (MQ=255)
                                                                                                                       cGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                              <  7:352998/148‑1 (MQ=255)
                                                                                                                             cAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCg                        >  2:191658/1‑148 (MQ=255)
                                                                                                                              aaTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCgg                       >  3:163685/1‑148 (MQ=255)
                                                                                                                               aTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                     <  8:315645/149‑1 (MQ=255)
                                                                                                                                tGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaca                    >  7:324262/1‑147 (MQ=255)
                                                                                                                                  ccccGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAc                   <  4:279821/148‑1 (MQ=255)
                                                                                                                                   cccgcccgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                  <  6:316341/148‑1 (MQ=255)
                                                                                                                                   cccgcccgGACCTGCATGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                  >  5:316341/1‑148 (MQ=255)
                                                                                                                                     cgcccgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATa               <  4:185074/149‑1 (MQ=255)
                                                                                                                                         cgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa            >  2:296900/1‑148 (MQ=255)
                                                                                                                                          ggACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaa          <  6:196156/149‑1 (MQ=255)
                                                                                                                                             cctgcctgAGAAATCGGTGCGCCAGAGGCAAAAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg       <  5:120920/149‑1 (MQ=255)
                                                                                                                                                gcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGGCGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTa    >  6:210713/1‑149 (MQ=255)
                                                                                                                                                 cctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCTTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCt     <  7:257146/147‑1 (MQ=255)
                                                                                                                                                  ctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc  >  1:40032/1‑149 (MQ=255)
                                                                                                                                                   tgAGAAATCGTTTAGCAATAGGCAACCAAAAGAAATTTACCGATTGAACCGCAAATGGCAAGCATCGTAAGGTGACGCGCCTTTAATTCACTGCGTAAGACCTGAGCTTCTGTTGTTTTAGTTTCGGAAACCAGACAAAAATGCTATc  >  8:267247/1‑148 (MQ=255)
                                                                                                                                                   tgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc  >  3:71841/1‑148 (MQ=255)
                                                                                                                                                   tgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc  >  5:241167/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATC  >  NC_000913/2248252‑2248546

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