Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,455,422 T→A 5.3% Q63L (CAG→CTG)  fbp ← fructose‑1,6‑bisphosphatase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,455,4220TA 5.3% 352.4 / 4.6 120Q63L (CAG→CTG) fbpfructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base T (63/50);  new base A (5/2);  total (68/52)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAACTCG  >  NC_000913/4455273‑4455571
                                                                                                                                                     |                                                                                                                                                     
ccAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTct                                                                                                                                                       >  2:60906/1‑150 (MQ=255)
 cAGGGGGTCCATCAGCACGAAGTATTTTGCGTTTTCACAGCCTTCATAGACGACAATATCATCTTCTTCTTCAGAGGGAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTctg                                                                                                                                                      <  1:154571/150‑1 (MQ=255)
 cAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGCGGCAATGCCCGCCACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTctg                                                                                                                                                      >  1:1084034/1‑150 (MQ=255)
 cAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCTAGTTTctg                                                                                                                                                      <  1:1676065/150‑1 (MQ=255)
  aGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTctgc                                                                                                                                                     <  2:1631760/150‑1 (MQ=255)
   gggggTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTctgct                                                                                                                                                    <  2:274903/150‑1 (MQ=255)
      ggTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGaa                                                                                                                                                 >  2:1333244/1‑150 (MQ=255)
          cATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGACGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTc                                                                                                                                             >  2:319602/1‑150 (MQ=255)
               gCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGcc                                                                                                                                          <  1:569319/148‑1 (MQ=255)
               gCACCACGTAGTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCt                                                                                                                                        <  1:1668536/150‑1 (MQ=255)
                caccacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGGTTCAGCAGAAACTCgacttg                                                                                                                                       <  2:2211/150‑6 (MQ=255)
                 accacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGc                                                                                                                                      >  1:1231081/1‑150 (MQ=255)
                 accacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGc                                                                                                                                      >  1:189138/1‑150 (MQ=255)
                   cacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGGGTTCAGGTTTTCATTAGCGAACAAGTCGAGTGTCTGCTGAACCTCGCCCTGCAc                                                                                                                                    >  2:118003/1‑150 (MQ=255)
                   cacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCAc                                                                                                                                    <  1:1650041/150‑1 (MQ=255)
                   cacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCAc                                                                                                                                    >  1:1221405/1‑150 (MQ=255)
                   caagTATTTTGCGGGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGACCGCAACGAAAACGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCAc                                                                                                                                    <  1:843072/147‑1 (MQ=255)
                    acGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACg                                                                                                                                   >  1:1302474/1‑150 (MQ=255)
                        tttttGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCATCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGATTTTCTGCTGAACCTCGCCCTGCACGTtct                                                                                                                               <  1:936644/149‑1 (MQ=255)
                        aTTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTtct                                                                                                                               <  2:275886/150‑1 (MQ=255)
                          tttGCGTGTTCACAGCCTTCAAAGACGACAATCTTATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCa                                                                                                                             <  1:1493154/150‑1 (MQ=255)
                          tttGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCa                                                                                                                             <  1:530223/150‑1 (MQ=255)
                           ttGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAg                                                                                                                            <  2:820126/150‑1 (MQ=255)
                           ttGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCGTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAg                                                                                                                            <  2:116097/150‑1 (MQ=255)
                              cGTGTTCACAGCCTTCAAAGATGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCAc                                                                                                                         >  2:1648993/1‑150 (MQ=255)
                                    cacaGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTgg                                                                                                                   >  2:1506071/1‑150 (MQ=255)
                                    cacaGCCTTCAAAGACGACAATCTAATATTGTTTTTCAGTGGCAATGCCCGCAACGATATCGCGAGCGTTACGTGCGGCTTTATTTTTTTCATTTGCGAACAAGTCGAGTTACTGGTGAACCTAGCCCTGAAGGTTCTCAGCCCCGCttg                                                                                                                   >  2:1210298/1‑148 (MQ=255)
                                     acaGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGc                                                                                                                  >  2:83280/1‑150 (MQ=255)
                                     acaGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTAAGTTTTTCATTAGCGAACAAGTCGAGTTTATGCTTAACCTCGCCCTGCACGTTCTCAGCACCGCTGGc                                                                                                                  >  1:419207/1‑150 (MQ=255)
                                        gCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAcc                                                                                                               <  1:1333244/150‑1 (MQ=255)
                                        gCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAcc                                                                                                               >  2:269590/1‑150 (MQ=255)
                                        gCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCGTGCAGTGTTTCATGAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCCCGTTCTCAGCACCGCTTGCAcc                                                                                                               >  2:1088790/1‑150 (MQ=255)
                                            tCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAgg                                                                                                           >  2:1291328/1‑150 (MQ=255)
                                              aaaGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGat                                                                                                         <  2:600184/150‑1 (MQ=255)
                                                aGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGatat                                                                                                       <  2:152323/150‑1 (MQ=255)
                                                  acgacAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCa                                                                                                     <  1:696937/150‑1 (MQ=255)
                                                     acAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAAcc                                                                                                  <  2:1703754/150‑1 (MQ=255)
                                                          ctcATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGCACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCACCCAGTcc                                                                                             >  2:104994/1‑150 (MQ=255)
                                                          ctcATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGCACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAGCCAGTcc                                                                                             >  2:693981/1‑150 (MQ=255)
                                                          ctcATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTcc                                                                                             <  2:631908/150‑1 (MQ=255)
                                                               cttcttcttCAGAGGCAATGCCCGCCACGATATGGCGTGCTTTCAGTGCGGCTTTCAGTGTTTCTTTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCtt                                                                                        <  2:304124/150‑1 (MQ=255)
                                                                    tcttcAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTttgttg                                                                                   <  1:60906/150‑1 (MQ=255)
                                                                    tcttcAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTttgttg                                                                                   >  1:103187/1‑150 (MQ=255)
                                                                    tcttcAGAGGCAATGCCCGCAACGATAGCGCGGGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTttgttg                                                                                   <  1:309655/150‑1 (MQ=255)
                                                                       tcAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACACGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGata                                                                                >  2:1693676/1‑150 (MQ=255)
                                                                           agGCAATGCCCGCAACGTTATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATcgc                                                                            <  1:1128652/150‑1 (MQ=255)
                                                                           agGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATcgc                                                                            <  2:211004/150‑1 (MQ=255)
                                                                            ggCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGATCCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATcgcg                                                                           >  2:1179390/1‑150 (MQ=255)
                                                                                   cccGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCACTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATa                                                                    >  1:1127195/1‑150 (MQ=255)
                                                                                    ccGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATaa                                                                   >  1:1670614/1‑150 (MQ=255)
                                                                                      gccaCGATAGCGCGGGCTGTCAGTGCGGCGTTCCGATATTCATTAGCCAACAAGGCGACTTTCTGCTGAACCTCGCCCTGCACGGTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTGGATATCGCGATGGATAATc                                                                 <  2:1034735/147‑1 (MQ=255)
                                                                                      gCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATc                                                                 >  2:1272362/1‑150 (MQ=255)
                                                                                        aaCGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCtt                                                               <  2:411429/150‑1 (MQ=255)
                                                                                           gATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGgc                                                            <  2:798855/150‑1 (MQ=255)
                                                                                            atatCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGgcg                                                           >  2:388710/1‑150 (MQ=255)
                                                                                                cgcgTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTACTGCTTTGTTGAAATCGCGATGGATAATCTTGGCGGCCa                                                       >  1:1078983/1‑150 (MQ=255)
                                                                                                  cgTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGCACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAACCAGGATATCAACCAGTCCGGCGTTGTTGAGATCGCGATGGATAATCTTGGCGCCCAGt                                                     >  2:1604351/1‑150 (MQ=255)
                                                                                                  cgTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGt                                                     >  2:1181362/1‑150 (MQ=255)
                                                                                                    tGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGATGAACCTCGCCCTGAACGTTCTCAGCACCGATGTCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGttt                                                   >  1:313799/1‑150 (MQ=255)
                                                                                                    ggCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACCAGTCGAGTTTCAGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGGCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGttt                                                   <  1:1088790/149‑1 (MQ=255)
                                                                                                       tttCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTTCTGAACCTCGCCCTGCACGTTCTCAGCACCGCAGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTAt                                                >  1:178462/1‑150 (MQ=255)
                                                                                                       tttCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTAt                                                <  2:1302474/150‑1 (MQ=255)
                                                                                                       tttCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTAt                                                >  1:1648834/1‑150 (MQ=255)
                                                                                                         tCAGTTCGGCTTTCAGTTTTTCTTTAGCGAAGAAGGCGAGGTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGTATATCAACCACTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTg                                              <  2:1443421/150‑1 (MQ=255)
                                                                                                         tCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATTTCGCGATGGATAATCTTGGCGCCCAGTTTTAttt                                              >  1:1283004/1‑149 (MQ=255)
                                                                                                         tCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTg                                              <  2:178972/150‑1 (MQ=255)
                                                                                                          cAGAGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTGTGCTGAACCTCGCACTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGc                                             <  2:50486/150‑1 (MQ=255)
                                                                                                            gTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCg                                           >  1:605991/1‑150 (MQ=255)
                                                                                                            gTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCg                                           >  1:19302/1‑150 (MQ=255)
                                                                                                              gCGTCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGAc                                         <  1:1125524/150‑1 (MQ=255)
                                                                                                               cGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGAcc                                        >  1:1617214/1‑149 (MQ=255)
                                                                                                                ggCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAg                                       >  2:1706685/1‑150 (MQ=255)
                                                                                                                     tCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAg                                  <  1:2327/150‑1 (MQ=255)
                                                                                                                      cAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGc                                 <  2:1630486/150‑1 (MQ=255)
                                                                                                                      cAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGc                                 <  1:644635/150‑1 (MQ=255)
                                                                                                                       aGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTAGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCa                                >  2:1447/1‑150 (MQ=255)
                                                                                                                           attCATTAGCGAACAAGTCGAGTATCTGATGAGCATAGACCAGCACGATCTCACCACCCAAGGCACCCAGGATATAAACCAGTCCTGATTTGATGATATCGCGATGCATAAACTAGGCGCCCAATATTAATGCAGACAACAACTCAGTGa                            >  2:428395/2‑150 (MQ=255)
                                                                                                                            ttCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAg                           >  1:442324/1‑150 (MQ=255)
                                                                                                                             tCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGc                          >  1:617521/1‑150 (MQ=255)
                                                                                                                                  aGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCAcc                     >  2:130852/1‑150 (MQ=255)
                                                                                                                                    cGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCgg                   >  1:1696271/1‑150 (MQ=255)
                                                                                                                                      tgcagggcGAGGTTCAGCAGAAACTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTa                 >  1:859839/8‑150 (MQ=255)
                                                                                                                                       aCAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAg                <  1:900968/150‑1 (MQ=255)
                                                                                                                                         agggcGAGGTTCAGCAGAAACTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCa              >  1:246423/5‑150 (MQ=255)
                                                                                                                                         aaGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGCCAGCAAAGCAGTGAGCTCACCGGTAGCa              >  1:792410/1‑150 (MQ=255)
                                                                                                                                          aGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCAt             >  2:191329/1‑150 (MQ=255)
                                                                                                                                           ggcGAGGTTCAGCAGAAACTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATg            >  1:1499599/3‑150 (MQ=255)
                                                                                                                                           gTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATg            >  1:1180191/1‑150 (MQ=255)
                                                                                                                                           gTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATg            >  2:485550/1‑150 (MQ=255)
                                                                                                                                           gTCGAGTTACTGCTGAACCTCGCAATTCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGCATAATCTTGGCGCCCTGTTTTATTGCCGACAGCTAAGCAGTGAGCTCACCGGTAGCATg            <  2:330550/150‑1 (MQ=255)
                                                                                                                                            tCGAGTTACTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAGCCAGGAAATCTAACAGACCTGCTTTCTGGATATCGCGATGCATAATCTGGGCGAGCAGTTTTATAGCCGACGGCAAAGCAGATAGCTCACAGGTAGCATga           >  2:308779/1‑150 (MQ=255)
                                                                                                                                             cGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATgag          >  2:1197963/1‑150 (MQ=255)
                                                                                                                                               aGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGaa        >  2:447205/1‑150 (MQ=255)
                                                                                                                                                gTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGaaa       >  2:1722942/1‑150 (MQ=255)
                                                                                                                                                  ttCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAAc     >  2:466480/1‑150 (MQ=255)
                                                                                                                                                  ttCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAAc     <  1:319602/150‑1 (MQ=255)
                                                                                                                                                  ttCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAAc     <  1:1315517/150‑1 (MQ=255)
                                                                                                                                                   tCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAACt    <  1:1128294/150‑1 (MQ=255)
                                                                                                                                                   tCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGAGTGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAACt    <  1:269590/150‑1 (MQ=255)
                                                                                                                                                     tgctgAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAACTCg  <  1:889544/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGCCCAGTTTTATTGCCGACAGCAAAGCAGTGAGCTCACCGGTAGCATGAGAAAACTCG  >  NC_000913/4455273‑4455571

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