Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,193,709 T→A 6.7% M154L (ATG→TTG)  mrp ← antiporter inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,193,7090TA 6.7% 317.1 / 9.2 90M154L (ATG→TTG) mrpantiporter inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/44);  new base A (4/2);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TGCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACAC  >  NC_000913/2193565‑2193855
                                                                                                                                                |                                                                                                                                                  
tGCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAttgtt                                                                                                                                                <  4:744501/149‑1 (MQ=255)
tGCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAttgt                                                                                                                                                 <  2:403054/148‑1 (MQ=255)
tGCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAttgt                                                                                                                                                 <  2:76457/148‑1 (MQ=255)
tGCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAttgt                                                                                                                                                 <  1:473416/148‑1 (MQ=255)
 gCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAttgttg                                                                                                                                               <  3:570143/149‑1 (MQ=255)
 gCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAttgtt                                                                                                                                                <  2:822893/148‑1 (MQ=255)
   tCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAATGTTggg                                                                                                                                             >  4:326551/1‑149 (MQ=255)
   tCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACAATGTTggg                                                                                                                                             <  3:326551/149‑1 (MQ=255)
     gCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGAt                                                                                                                                           >  3:690663/1‑149 (MQ=255)
        aTCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGa                                                                                                                                        <  3:137460/149‑1 (MQ=255)
           ggTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATgg                                                                                                                                     <  3:890161/149‑1 (MQ=255)
            gTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATggg                                                                                                                                    <  3:223512/149‑1 (MQ=255)
               cACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCg                                                                                                                                 >  3:613347/1‑149 (MQ=255)
                aCGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGa                                                                                                                                            >  1:874383/1‑137 (MQ=255)
                aCGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGa                                                                                                                                            <  2:874383/137‑1 (MQ=255)
                 cGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGt                                                                                                                                <  2:860784/148‑1 (MQ=255)
                 cGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGt                                                                                                                                <  1:875066/148‑1 (MQ=255)
                 cGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGt                                                                                                                                <  2:431224/148‑1 (MQ=255)
                 cGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTa                                                                                                                               <  4:803278/149‑1 (MQ=255)
                  gCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAg                                                                                                                              <  4:838434/149‑1 (MQ=255)
                   ccACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGa                                                                                                                             >  3:660681/1‑149 (MQ=255)
                   ccACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGa                                                                                                                             <  4:704583/149‑1 (MQ=255)
                       acCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATc                                                                                                                         >  4:646320/1‑149 (MQ=255)
                        ccATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATc                                                                                                                         <  2:671586/148‑1 (MQ=255)
                        ccATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATc                                                                                                                         <  1:896199/148‑1 (MQ=255)
                        ccATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATc                                                                                                                         <  2:31878/148‑1 (MQ=255)
                        ccATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATc                                                                                                                         >  2:142989/1‑148 (MQ=255)
                         cATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCgg                                                                                                                       >  4:674626/1‑149 (MQ=255)
                             gCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTc                                                                                                                   <  3:674626/149‑1 (MQ=255)
                               aTTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTcc                                                                                                                  >  1:120321/1‑148 (MQ=255)
                                ttGTCCTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCa                                                                                                                 >  2:648177/1‑148 (MQ=255)
                                  gtcgtcAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGaa                                                                                                              <  3:190796/149‑1 (MQ=255)
                                  gtcgtcAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGa                                                                                                               <  2:75452/148‑1 (MQ=255)
                                  gtcgtcAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGa                                                                                                               <  2:763743/148‑1 (MQ=255)
                                    cgtcAGTGACCAGATAACCAATTGAAGTGGGTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATCTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATg                                                                                                            <  3:549230/149‑1 (MQ=255)
                                      tcAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGcc                                                                                                          <  3:885983/149‑1 (MQ=255)
                                         gTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGAc                                                                                                       <  3:847901/149‑1 (MQ=255)
                                          ggACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCTGGCGAAGTTGGACGCTGATTTTCTGCGCCCAGCATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGAc                                                                                                       >  2:600844/2‑148 (MQ=255)
                                            aCCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAAt                                                                                                             >  4:902607/1‑140 (MQ=255)
                                            aCCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAAt                                                                                                             <  3:902607/140‑1 (MQ=255)
                                               aGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGc                                                                                                 <  4:861147/149‑1 (MQ=255)
                                                gATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATgg                                                                                                                                     <  3:878917/112‑1 (MQ=255)
                                                gATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATgg                                                                                                                                     >  4:878917/1‑112 (MQ=255)
                                                 aTAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGcc                                                                                                <  1:38657/148‑1 (MQ=255)
                                                   aaCCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGATGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGCTGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGcccc                                                                                              <  1:666952/148‑1 (MQ=255)
                                                    aCCAAGTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCt                                                                                             <  2:2327/148‑1 (MQ=255)
                                                      cAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTc                                                                                           <  2:786638/148‑1 (MQ=255)
                                                      cAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCt                                                                                          >  3:677306/1‑149 (MQ=255)
                                                       aattgaattgGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTg                                                                                         >  4:870727/1‑149 (MQ=255)
                                                        attgaattgGTTGCCAGGACATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTg                                                                                         >  1:895335/1‑148 (MQ=255)
                                                         ttgaattgGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTgcc                                                                                       <  3:18662/149‑1 (MQ=255)
                                                            aattgGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTgccgcc                                                                                    >  3:411980/1‑149 (MQ=255)
                                                            aattgGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTgccgcc                                                                                    >  3:358187/1‑149 (MQ=255)
                                                             attgGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCcgcca                                                                                   >  3:567648/1‑149 (MQ=255)
                                                              ttgGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCcgcca                                                                                   >  2:350246/1‑148 (MQ=255)
                                                                    gCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCcgccagcgccag                                                                            <  4:135615/149‑1 (MQ=255)
                                                                       aGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTg                                                                          >  1:26076/1‑148 (MQ=255)
                                                                       aGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTg                                                                          >  2:813483/1‑148 (MQ=255)
                                                                           cATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAgg                                                                     >  4:496010/1‑149 (MQ=255)
                                                                           cATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAg                                                                      <  1:648177/148‑1 (MQ=255)
                                                                              gagaCATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAAt                                                                                                             <  4:426978/106‑1 (MQ=255)
                                                                              gagaCATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAAt                                                                                                             >  3:426978/1‑106 (MQ=255)
                                                                                    tAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCATAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCa                                                                       <  1:66276/138‑1 (MQ=255)
                                                                                    tAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCATAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCa                                                                       >  2:66276/1‑138 (MQ=255)
                                                                                     aaTCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTAccgccg                                                           >  3:461601/1‑149 (MQ=255)
                                                                                      aTCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTAccgccg                                                           >  2:99349/1‑148 (MQ=255)
                                                                                             ccATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGa                                                    <  1:350246/148‑1 (MQ=255)
                                                                                              cATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGAtt                                                  >  3:734420/1‑149 (MQ=255)
                                                                                                tatGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGAttt                                                 >  1:229521/1‑148 (MQ=255)
                                                                                                 atGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTAc                                               <  4:613347/149‑1 (MQ=255)
                                                                                                  tGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTAcc                                              <  4:91325/149‑1 (MQ=255)
                                                                                                   gagTACCGTCCGGCGAAGTTGGACGCTGATGTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTAccc                                             <  3:700538/147‑1 (MQ=255)
                                                                                                      tACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCAcg                                          >  4:832458/1‑149 (MQ=255)
                                                                                                         cGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCAcgccg                                       <  4:411980/149‑1 (MQ=255)
                                                                                                          gTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATGGGCGTCCAGAATGCCGACTATCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCAcgccgc                                      >  3:556124/1‑149 (MQ=255)
                                                                                                               gCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCttt                                  >  1:86021/1‑148 (MQ=255)
                                                                                                                cGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTa                                 >  1:28086/1‑148 (MQ=255)
                                                                                                                cGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTAGGGATTGATGGGCCGTAGATATCGGCTTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTAGCGCCGTGGACGATTTACCCACGCCGCCTTTAc                                >  3:563719/1‑149 (MQ=255)
                                                                                                                     ttGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGa                            <  1:699608/148‑1 (MQ=255)
                                                                                                                           gCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCAcc                     >  4:44921/1‑149 (MQ=255)
                                                                                                                                 tttCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGAt                <  1:813483/148‑1 (MQ=255)
                                                                                                                                       cgcCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttc         <  3:870727/149‑1 (MQ=255)
                                                                                                                                       atcccAACAATGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttc         <  4:903171/147‑1 (MQ=255)
                                                                                                                                       atcccAACAATGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttc         >  3:903171/3‑149 (MQ=255)
                                                                                                                                         cccAACAATGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttct        >  1:234215/1‑148 (MQ=255)
                                                                                                                                         cccAACAATGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttct        >  2:416806/1‑148 (MQ=255)
                                                                                                                                          ccAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttcttc      >  4:782542/1‑149 (MQ=255)
                                                                                                                                          ccAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttctt       <  3:258063/148‑1 (MQ=255)
                                                                                                                                          ccAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttctt       >  4:258063/1‑148 (MQ=255)
                                                                                                                                           cAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTca     <  3:782542/149‑1 (MQ=255)
                                                                                                                                            aaCATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTcac    >  4:819951/1‑149 (MQ=255)
                                                                                                                                            aaCATTGTTGGAATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTcac    >  3:529024/1‑149 (MQ=255)
                                                                                                                                             aCATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTcac    >  1:891705/1‑148 (MQ=255)
                                                                                                                                              cATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCa                                            >  1:899830/1‑107 (MQ=255)
                                                                                                                                              cATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCa                                            <  2:899830/107‑1 (MQ=255)
                                                                                                                                              cATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCacac  <  4:358187/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TGCTCGCCATCGGTCCACGCCACACCATCGCATTGTCGTCAGTGACCAGATAACCAATTGAATTGGTTGCCAGGCCATGAGACATAATCGGTGCCATATGGGTACCGTCCGGCGAAGTTGGACGCTGATTTTCTGCGCCCAACATTGTTGGGATTGATGGGCCGTAGATATCGGCGTCCAGAATGCCGACTTTCGCCCCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACAC  >  NC_000913/2193565‑2193855

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: