Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 895,932 G→A 7.1% intergenic (+131/+12) metQ → / ← USA300HOU_0850 ABC superfamily ATP binding cassette transporter, binding protein/integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730895,9320GA 7.1% 242.9 / 5.2 83intergenic (+131/+12)metQ/USA300HOU_0850ABC superfamily ATP binding cassette transporter, binding protein/integrase
Reads supporting (aligned to +/- strand):  ref base G (40/36);  new base A (2/4);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCATAATAGCTTTA  >  CP000730/895794‑896067
                                                                                                                                          |                                                                                                                                       
cTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGAt                                                                                                                                       >  1:45502/1‑141 (MQ=255)
   aataatACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTg                                                                                                                                    <  1:51324/141‑1 (MQ=255)
   aataatACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTg                                                                                                                                    <  1:51321/141‑1 (MQ=255)
   aagaagaCCAATTAAATAAATAGACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGAGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTg                                                                                                                                    <  2:313004/135‑1 (MQ=255)
      aataCCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGTGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGAGTTGAgc                                                                                                                                 >  1:473216/1‑141 (MQ=255)
      aataCCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgc                                                                                                                                 >  1:324641/1‑141 (MQ=255)
       ataCCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                                                                >  1:384629/1‑141 (MQ=255)
       ataCCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcg                                                                                                                                <  1:713410/141‑1 (MQ=255)
        taCCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATGCGCTCAAATCTGAgcgt                                                                                                                               <  1:171618/141‑2 (MQ=255)
          ccAATTAAATAAAAATACCCTTTAATAATAAGTATTAAGTGGTTATTTTTGGAGTATTAAAAATAAAAAAAGGCACAAAAGTGACAGTATATACAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAGCGTAt                                                                                                                             <  2:570143/141‑1 (MQ=255)
          ccAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                                                             >  1:532790/1‑141 (MQ=255)
            aaTTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATGCGCTCAAATCTGAGCGTAttt                                                                                                                           <  2:560453/141‑2 (MQ=255)
            aaTTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGAAAAAAAAAAAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTa                                                                                                                           <  2:403592/141‑1 (MQ=255)
              ttAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATc                                                                                                                         >  2:373032/1‑141 (MQ=255)
               taaataaaAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTTTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                                                                        >  2:301775/1‑141 (MQ=255)
               taaataaaAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                                                                        <  1:484007/141‑1 (MQ=255)
               taaataaaAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGAAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                                                                        >  2:475458/1‑141 (MQ=255)
                   taaaAATACCCTTTAATAACAACGATTAAGGGGGTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGt                                                                                                                    >  1:409013/1‑141 (MQ=255)
                     aaaaTACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGGGCATTATACAAGttt                                                                                                                  >  1:437249/1‑141 (MQ=255)
                     aaaaTACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCCTTATCCAAGttt                                                                                                                  >  1:237297/1‑141 (MQ=255)
                      aaaTACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtttt                                                                                                                 >  2:168462/1‑141 (MQ=255)
                      aaaTACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTAAAAATAAAAAAGGGCAAAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtttt                                                                                                                 >  2:123159/1‑141 (MQ=255)
                          aCCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCt                                                                                                             >  2:373242/1‑141 (MQ=255)
                           cccTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                                                            >  1:698332/1‑141 (MQ=255)
                            ccTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                                                           <  2:398147/141‑1 (MQ=255)
                            ccTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGAAGTATATGGAAATTATAGAAATTTAAACAAAAAAATAAGATCAGACTTGAGCGCATTATCCACGATTTACCTGt                                                                                                           >  1:72683/1‑141 (MQ=255)
                             cTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGtt                                                                                                          >  2:397888/1‑141 (MQ=255)
                             cTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGtt                                                                                                          <  1:389252/141‑1 (MQ=255)
                                 aataaCAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                                                                      <  1:296743/141‑1 (MQ=255)
                                  ataaCAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGt                                                                                                     <  1:413179/141‑1 (MQ=255)
                                   taaCAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                                                                    <  1:223448/141‑1 (MQ=255)
                                     acaacGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCTAATTCTCGAATTTTAAACAGAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                                                                  >  1:664278/1‑141 (MQ=255)
                                       aacGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGtt                                                                                                >  1:606496/1‑141 (MQ=255)
                                        acGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTTAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                                                               <  1:567014/141‑1 (MQ=255)
                                        acGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                                                               <  2:776293/141‑1 (MQ=255)
                                          gATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                                                             >  2:155063/1‑141 (MQ=255)
                                           aTTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatc                                                                                            <  1:522583/141‑1 (MQ=255)
                                            ttAAGGGGTTATTTTCGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                                                                           <  1:653495/141‑1 (MQ=255)
                                                ggggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGATCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                                                                       <  1:724345/141‑1 (MQ=255)
                                                     tATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                                                                  <  1:179356/141‑1 (MQ=255)
                                                       tttgtGTGGTATTAAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                                                                <  1:556609/141‑1 (MQ=255)
                                                            gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGCTCTGAGTCCCTTATCAAAGTATTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                                                           >  1:703728/1‑141 (MQ=255)
                                                             tgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCAAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                                                                          <  2:2522/141‑1 (MQ=255)
                                                                tATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCa                                                                       <  1:44772/141‑1 (MQ=255)
                                                                       aaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                                                                <  2:242920/141‑1 (MQ=255)
                                                                       aaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                                                                <  2:366608/141‑1 (MQ=255)
                                                                       aaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                                                                <  1:54290/141‑1 (MQ=255)
                                                                        aTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTg                                                               <  2:592650/141‑1 (MQ=255)
                                                                        aTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTg                                                               <  1:448692/141‑1 (MQ=255)
                                                                           aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgaga                                                            <  1:109277/141‑1 (MQ=255)
                                                                             aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt                                                          >  1:313876/1‑141 (MQ=255)
                                                                             aaaGAGCACAAAAGGAGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt                                                          <  1:448067/141‑1 (MQ=255)
                                                                              aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATa                                                         >  2:280263/1‑141 (MQ=255)
                                                                              aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATa                                                         >  2:585787/1‑141 (MQ=255)
                                                                                gggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGTGCTTTATCGAAGTTTTACCTTTTCAAGTTGGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa                                                       >  1:32673/1‑141 (MQ=255)
                                                                                  gCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCGAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg                                                     <  2:474220/141‑1 (MQ=255)
                                                                                    acaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa                                                   <  1:132465/141‑1 (MQ=255)
                                                                                        aaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgt                                               <  1:373242/141‑1 (MQ=255)
                                                                                         aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgtg                                              >  2:225240/1‑141 (MQ=255)
                                                                                         aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgtg                                              >  1:778825/1‑141 (MQ=255)
                                                                                         aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgtg                                              >  1:778701/1‑141 (MQ=255)
                                                                                           gggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGtgtt                                            >  1:358189/1‑141 (MQ=255)
                                                                                            ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGtgttg                                           <  1:90231/141‑1 (MQ=255)
                                                                                             gCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGtgttgt                                          >  1:606329/1‑141 (MQ=255)
                                                                                              aaatataTGCACATTATCGAATTTTAAACAAAAAAATACGCTCATATTTTAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTTAGAATAAATGGAAAGTGTTGTa                                         <  2:2794/138‑1 (MQ=255)
                                                                                                         aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGa                              <  2:532790/141‑1 (MQ=255)
                                                                                                          aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGaa                             >  1:438997/1‑141 (MQ=255)
                                                                                                             aTCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATgg                          >  1:718852/1‑141 (MQ=255)
                                                                                                             aTCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATgg                          >  1:382523/1‑141 (MQ=255)
                                                                                                              tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGc                         >  1:584648/1‑141 (MQ=255)
                                                                                                               cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGcc                        >  1:755155/1‑141 (MQ=255)
                                                                                                                gAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCt                       <  1:155063/141‑1 (MQ=255)
                                                                                                                 aaTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTa                      <  1:125314/141‑1 (MQ=255)
                                                                                                                  aTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTac                     <  1:739610/141‑1 (MQ=255)
                                                                                                                       aaaCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGAt                >  1:286276/1‑141 (MQ=255)
                                                                                                                        aaCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATc               >  1:399829/1‑141 (MQ=255)
                                                                                                                           aaaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCat            <  1:373032/141‑1 (MQ=255)
                                                                                                                              aaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCataat         <  1:722611/141‑1 (MQ=255)
                                                                                                                                  tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCATAAAAGCt     >  2:55398/1‑141 (MQ=255)
                                                                                                                                   gcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCATAATAGCtt    >  2:303018/2‑141 (MQ=255)
                                                                                                                                   aCGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAACGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCATAATAGCtt    >  2:66144/1‑141 (MQ=255)
                                                                                                                                    cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCATAATAGCttt   >  2:447042/1‑141 (MQ=255)
                                                                                                                                     gCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCATAATAGCTTTa  >  2:174227/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                       
CTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTACGGTGGTCGGAATGGCCTACACGATCCATAATAGCTTTA  >  CP000730/895794‑896067

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