Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,077,754 G→A A62A (GCG→GCA cbtA → CP4‑44 prophage; toxin of the YeeV‑YeeU toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,077,7540GA98.4% 228.9 / ‑4.5 63A62A (GCG→GCAcbtACP4‑44 prophage; toxin of the YeeV‑YeeU toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/39);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CGCCCGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAC  >  NC_000913/2077611‑2077895
                                                                                                                                               |                                                                                                                                             
cGCCCGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGaa                                                                                                                                          <  2:813513/149‑1 (MQ=255)
 gCCCGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACtt                                                                                                                                       >  2:594340/1‑151 (MQ=255)
   ccGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTtc                                                                                                                                     <  2:194695/151‑1 (MQ=255)
   ccGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTtc                                                                                                                                     <  2:891047/151‑1 (MQ=255)
       ctcCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGt                                                                                                                                 >  2:841037/1‑151 (MQ=255)
         ccTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCg                                                                                                                                  <  1:122711/148‑1 (MQ=255)
         ccTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTgg                                                                                                                               >  1:137714/1‑151 (MQ=255)
              tGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGaaaa                                                                                                                           <  1:96176/150‑1 (MQ=255)
                      ggCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACgcgc                                                                                                                   <  1:375808/150‑1 (MQ=255)
                        cAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTgg                                                                                                                <  1:546978/151‑1 (MQ=255)
                              ctgctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGt                                                                                                           >  2:908964/1‑150 (MQ=255)
                              ctgctgTCCCGACTGCGGGACCAGCACTATGGCCTCACACTGACTGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCACGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGt                                                                                                           >  2:101778/1‑150 (MQ=255)
                                gctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTa                                                                                                          <  1:884282/149‑1 (MQ=255)
                                gctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTAcc                                                                                                        >  2:781197/1‑151 (MQ=255)
                                gctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTAcc                                                                                                        <  1:89051/151‑1 (MQ=255)
                                 ctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTtc                                                                                                                                     >  1:807558/1‑121 (MQ=255)
                                 ctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTtc                                                                                                                                     <  2:807493/121‑1 (MQ=255)
                                        gACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGcc                                                                                                <  1:694510/151‑1 (MQ=255)
                                        gACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGcc                                                                                                >  2:92916/1‑151 (MQ=255)
                                               ggACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAtt                                                                                          <  1:841104/150‑1 (MQ=255)
                                                      cACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCt                                                                                  <  2:823847/151‑1 (MQ=255)
                                                      cACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCt                                                                                  <  1:759516/151‑1 (MQ=255)
                                                         tATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGTTGGGGGGTACGGACCGCCCGGGATTCAgcgc                                                                                    >  2:594195/1‑146 (MQ=255)
                                                         tATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                >  2:644217/1‑150 (MQ=255)
                                                          aTGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                <  1:456029/149‑1 (MQ=255)
                                                          aTGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                <  2:703427/149‑1 (MQ=255)
                                                             gCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCg                                                                           <  1:15386/151‑1 (MQ=255)
                                                                 cacaCTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGctc                                                                        <  1:563925/150‑1 (MQ=255)
                                                                   cacTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCtctc                                                                      >  2:98963/1‑150 (MQ=255)
                                                                         aTGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaat                                                               >  1:203448/1‑151 (MQ=255)
                                                                            acacacCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaataa                                                             <  1:707733/150‑1 (MQ=255)
                                                                                acCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGc                                                        >  2:898883/1‑151 (MQ=255)
                                                                                 ccTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGc                                                        <  2:213260/150‑1 (MQ=255)
                                                                                      tGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGa                                                   <  2:431582/150‑1 (MQ=255)
                                                                                       gCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGata                                                 <  1:850256/151‑1 (MQ=255)
                                                                                        ccGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGatat                                                <  1:557979/151‑1 (MQ=255)
                                                                                          gATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtcc                                              <  2:884355/151‑1 (MQ=255)
                                                                                           aTGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtcct                                             <  2:533812/151‑1 (MQ=255)
                                                                                           aTGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtcct                                             <  1:874424/151‑1 (MQ=255)
                                                                                            tGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtcct                                             <  2:886122/150‑1 (MQ=255)
                                                                                             gAACGTGAGATTGAGCAGCATATTGAGGCAGGCATTTCACTGAGTGATGCGGTGAATTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACACGCGCTCAGTTAATCAACAGCATCGCTAtcctcc                                           >  2:851448/1‑151 (MQ=255)
                                                                                              aaCGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCg                                          <  2:896159/151‑1 (MQ=255)
                                                                                              aaCGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCg                                          <  2:46098/151‑1 (MQ=255)
                                                                                                      aTTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCa                                  <  1:315973/151‑1 (MQ=255)
                                                                                                        tGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAgg                                >  1:381451/1‑151 (MQ=255)
                                                                                                        tGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAg                                 <  2:474145/150‑1 (MQ=255)
                                                                                                         gAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGTCCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAgg                                <  2:579235/150‑1 (MQ=255)
                                                                                                         gAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAggg                               >  2:287795/1‑151 (MQ=255)
                                                                                                               caTATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCg                         <  2:102426/151‑1 (MQ=255)
                                                                                                                atatTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCgg                        >  2:89743/1‑151 (MQ=255)
                                                                                                                atatTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCgg                        <  2:13406/151‑1 (MQ=255)
                                                                                                                   ttGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCt                     >  2:337115/1‑151 (MQ=255)
                                                                                                                    tGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtg                    >  2:136798/1‑151 (MQ=255)
                                                                                                                    tGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCt                     <  1:8465/150‑1 (MQ=255)
                                                                                                                           ggcaTTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGAccc             <  2:24278/151‑1 (MQ=255)
                                                                                                                              aTTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCcgcg          >  1:79546/1‑151 (MQ=255)
                                                                                                                                 tCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCcgcg          >  2:848404/1‑148 (MQ=255)
                                                                                                                                    cTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAAtt    >  2:865511/1‑151 (MQ=255)
                                                                                                                                      gtgtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAAtt    <  2:467110/149‑1 (MQ=255)
                                                                                                                                      gtgtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAAtt    <  2:322205/149‑1 (MQ=255)
                                                                                                                                      gtgtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTa   <  1:150214/150‑1 (MQ=255)
                                                                                                                                       tgtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAc  >  2:697241/1‑150 (MQ=255)
                                                                                                                                        gtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAc  >  2:527534/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CGCCCGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAC  >  NC_000913/2077611‑2077895

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