Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,474,485 G→A A23V (GCG→GTG)  rpsL ← 30S ribosomal subunit protein S12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,474,4850GA98.6% 281.3 / ‑3.8 72A23V (GCG→GTG) rpsL30S ribosomal subunit protein S12
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (30/41);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

ACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTT  >  NC_000913/3474340‑3474630
                                                                                                                                                 |                                                                                                                                                 
aCCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAgg                                                                                                                                               <  1:420155/150‑1 (MQ=255)
aCCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGc                                                                                                                                              <  1:920789/151‑1 (MQ=255)
    tCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCAc                                                                                                                                            <  1:540688/149‑1 (MQ=255)
    tCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACg                                                                                                                                           <  1:341808/150‑1 (MQ=255)
     caccacCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGtt                                                                                                                                         >  2:864021/1‑151 (MQ=255)
      accaccGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTg                                                                                                                                        >  1:5816/1‑151 (MQ=255)
      accaccGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTg                                                                                                                                        <  2:22128/151‑1 (MQ=255)
       ccaccGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGc                                                                                                                                       <  1:503552/151‑1 (MQ=255)
          ccGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCttt                                                                                                                                    >  2:924035/1‑151 (MQ=255)
                tAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCa                                                                                                                               <  2:364070/150‑1 (MQ=255)
                 aGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCaa                                                                                                                              <  1:430202/150‑1 (MQ=255)
                  ggAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACt                                                                                                                            <  2:613532/151‑1 (MQ=255)
                  ggAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACt                                                                                                                            <  1:1025243/151‑1 (MQ=255)
                   gAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACtt                                                                                                                           <  2:526539/151‑1 (MQ=255)
                    aaGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACttt                                                                                                                          <  2:569106/151‑1 (MQ=255)
                        cACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCg                                                                                                                       <  2:1041732/150‑1 (MQ=255)
                         aCTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATTCACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAg                                                                                                                     <  1:741784/151‑1 (MQ=255)
                          cTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAg                                                                                                                     <  1:1103540/150‑1 (MQ=255)
                              gAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGt                                                                                                                >  2:1066030/1‑151 (MQ=255)
                               aaaCCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTg                                                                                                               <  2:244047/151‑1 (MQ=255)
                                aaCCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTgg                                                                                                              >  1:16227/1‑151 (MQ=255)
                                   cGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGttt                                                                                                           <  2:984682/151‑1 (MQ=255)
                                   cGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGttt                                                                                                           >  1:382937/1‑151 (MQ=255)
                                    gTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTg                                                                                                          <  1:337960/151‑1 (MQ=255)
                                    gTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTg                                                                                                          >  1:956805/1‑151 (MQ=255)
                                       aGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGt                                                                                                       <  1:661685/151‑1 (MQ=255)
                                          cAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTAcc                                                                                                    >  1:755192/1‑151 (MQ=255)
                                            gACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTAcc                                                                                                    >  2:518191/1‑149 (MQ=255)
                                               gAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTg                                                                                               >  2:146183/1‑151 (MQ=255)
                                               gAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTg                                                                                               <  2:69640/151‑1 (MQ=255)
                                                aaCACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTgg                                                                                              <  2:780891/151‑1 (MQ=255)
                                                  cacGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGtt                                                                                            <  1:990215/151‑1 (MQ=255)
                                                       cATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTg                                                                                       >  2:302645/1‑151 (MQ=255)
                                                       cATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTg                                                                                       >  1:1057981/1‑151 (MQ=255)
                                                       cATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTg                                                                                       <  2:1081298/151‑1 (MQ=255)
                                                         tACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGtt                                                                                     >  2:236845/1‑151 (MQ=255)
                                                          aCTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGtt                                                                                     <  2:607760/150‑1 (MQ=255)
                                                          aCTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTg                                                                                    <  2:1074600/151‑1 (MQ=255)
                                                            tttACGCCGCGCGGAGGTCGGCTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACATGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTg                                                                                    >  1:743978/1‑149 (MQ=255)
                                                              tACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCa                                                                                 <  2:1017974/150‑1 (MQ=255)
                                                              tACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCAt                                                                                >  2:153178/1‑151 (MQ=255)
                                                                cGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTa                                                                              <  2:203179/151‑1 (MQ=255)
                                                                cGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCAtt                                                                               >  1:681720/1‑150 (MQ=255)
                                                                cGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCAt                                                                                <  1:398283/149‑1 (MQ=255)
                                                                      gcgGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAg                                                                         <  1:732514/150‑1 (MQ=255)
                                                                         gAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTcc                                                                     <  1:406829/151‑1 (MQ=255)
                                                                             tCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGt                                                                 >  2:530749/1‑151 (MQ=255)
                                                                                 ttttttAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGtttt                                                              >  1:551296/1‑150 (MQ=255)
                                                                                  tttttAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTa                                                             <  2:673558/150‑1 (MQ=255)
                                                                                    tttAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGc                                                           <  2:725785/150‑1 (MQ=255)
                                                                                    tttAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCt                                                          >  1:508589/1‑151 (MQ=255)
                                                                                       aGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCtttt                                                       >  1:293601/1‑151 (MQ=255)
                                                                                       aGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCtttt                                                       <  1:864961/151‑1 (MQ=255)
                                                                                            tGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTc                                                  >  1:1012314/1‑151 (MQ=255)
                                                                                            tGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTc                                                  >  1:825416/1‑151 (MQ=255)
                                                                                             ggTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTc                                                  <  1:22965/150‑1 (MQ=255)
                                                                                             ggTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCg                                                 <  2:48847/151‑1 (MQ=255)
                                                                                             ggTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCg                                                 <  1:280964/151‑1 (MQ=255)
                                                                                                  tatataCACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAac                                            >  1:678637/1‑151 (MQ=255)
                                                                                                   atataCACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAac                                            >  2:846685/1‑150 (MQ=255)
                                                                                                     ataCACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACg                                         <  1:781363/151‑1 (MQ=255)
                                                                                                        cacGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGta                                       >  1:367978/1‑150 (MQ=255)
                                                                                                         acGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGtaat                                     >  1:597450/1‑151 (MQ=255)
                                                                                                           gAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGtaataa                                   <  1:1021626/151‑1 (MQ=255)
                                                                                                            aGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGtaataa                                   >  1:646866/1‑150 (MQ=255)
                                                                                                                 aTACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTc                             <  1:683103/151‑1 (MQ=255)
                                                                                                                 aTACGCCACGTTTTTGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTc                             >  2:950533/1‑151 (MQ=255)
                                                                                                                             tttGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAatat                 >  2:434991/1‑151 (MQ=255)
                                                                                                                              ttGCGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAatat                 <  2:109724/150‑1 (MQ=255)
                                                                                                                                 cGGGCATGCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAgg             >  2:425079/1‑151 (MQ=255)
                                                                                                                                        gCTTCCAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGa      <  2:1043453/151‑1 (MQ=255)
                                                                                                                                            ccAGCACAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAAttt  <  1:30925/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTT  >  NC_000913/3474340‑3474630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: