Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,716,758 C→T 10.8% E77K (GAA→AAA)  USA300HOU_1609 ← possible allophanate hydrolase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,716,7580CT10.8% 182.3 / 3.0 73E77K (GAA→AAA) USA300HOU_1609possible allophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (28/35);  new base T (7/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 8.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTATGTTTTTCAATAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATATTTCTTTTTCTCGTATCCACCGTAC  >  CP000730/1716618‑1716896
                                                                                                                                            |                                                                                                                                          
tGTATGTTTTTCAATAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTc                                                                                                                                            >  1:462401/1‑141 (MQ=255)
tGTATGTTTTTCAATAACTTGCTCAACAGATATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATAGGTTTGTTTCACAGCATGAACAATTTTTTc                                                                                                                                            >  2:595118/1‑141 (MQ=255)
 gTATGTTTTTCAATAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCa                                                                                                                                           >  1:785134/1‑141 (MQ=255)
  tATGTTTTTCAATAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCAt                                                                                                                                          >  2:342228/1‑141 (MQ=255)
   atgctcttCACTACCTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATAGCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCAtt                                                                                                                                         <  1:269983/135‑1 (MQ=255)
       ttttCAATAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAAt                                                                                                                                     <  1:12100/141‑1 (MQ=255)
        tttCAATAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAAtt                                                                                                                                    >  1:210836/1‑141 (MQ=255)
             aTAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTa                                                                                                                               <  1:508354/141‑1 (MQ=255)
              tAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTaa                                                                                                                              <  1:5058/141‑1 (MQ=255)
               aaCTTGCTCAACAGTTATTCGATTACGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACGACCATATAGAACAGGTATTTTAATAACACGCGTGGTTTGATTCACAGCATGAACATTTTTTTTACATAATTTGTTaga                                                                                                                             >  2:192332/1‑139 (MQ=255)
                 cTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTAGATTCACAACATGAACATTTTTTTCATTAAATTTGTTAAGtt                                                                                                                           >  1:735427/1‑141 (MQ=255)
                   tGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCt                                                                                                                         >  2:45888/1‑141 (MQ=255)
                       taataGTTATTCGTTTATGCTTTGCTACTTCTTCAATATTTGGTCCATATGTACCATCTTATTGAACAGGTATTTTAATATTACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAAt                                                                                                                     <  2:549176/137‑1 (MQ=255)
                         aCATTTATTCGATTATGCTTTGCTACCTCTTCAATATCTGGTTCATATGTACCTACATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTTACAGCATGAACACTTTTTTCATTAAATTTGTTAAGTTCTAAAttt                                                                                                                   <  2:706607/141‑1 (MQ=255)
                               aTTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCaatt                                                                                                             >  1:919/1‑141 (MQ=255)
                               aTTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTCTTTTCATTAAATTTGTTAAGTTCTAAATTTTCaatt                                                                                                             >  1:438624/1‑141 (MQ=255)
                                  cGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCaattaat                                                                                                          >  1:675167/1‑141 (MQ=255)
                                    tttATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCATCATATTGAACAGTTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAAttt                                                                                                        <  2:3020/140‑1 (MQ=255)
                                    aTTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAAttt                                                                                                        <  2:684528/141‑1 (MQ=255)
                                     ttctGCTTTTCTACTCCTTCTATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTTATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTa                                                                                                       <  2:27299/138‑1 (MQ=255)
                                     ttATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTa                                                                                                       <  1:625541/141‑1 (MQ=255)
                                      tATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTag                                                                                                      <  2:679134/141‑1 (MQ=255)
                                         gCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTTTTTTTTTTTTTTAAGTTCTAAATTTTCAATTTATTTagaat                                                                                                   >  2:567625/1‑140 (MQ=255)
                                                tACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTCAATTTGTTAAGTTCTAAATTTTCAATTACTTTGGAAGTGGTtat                                                                                            >  2:753728/1‑141 (MQ=255)
                                                    tcttcAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCt                                                                                        >  1:396266/1‑141 (MQ=255)
                                                     tttcAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTTATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATTGCtt                                                                                       <  2:481591/140‑1 (MQ=255)
                                                      ttcAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCTTTATATTTATTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATGGCttt                                                                                      >  2:404843/1‑140 (MQ=255)
                                                      ttcAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTg                                                                                      >  1:349976/1‑141 (MQ=255)
                                                       tcAATATCTGGTCCATATGTCCCCAAACAGTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTTTTTATTTTTGTTAAGTTCTAAATTTTAAATTAATTAAGAAGTTGTTTTGGCTTGt                                                                                     >  2:559641/1‑141 (MQ=255)
                                                           tatCTGGTCCATATGTACCACCATATTGAATAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTATTTAAATTTGTTAATTTCTCAATTTACAATTAATTTATATGATGttttttcttcttctt                                                                                 >  2:69070/1‑128 (MQ=255)
                                                              ttGGTTCATATGTACCACCATCTTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCaa                                                                              <  1:42044/140‑1 (MQ=255)
                                                              ttGGTCCATATGGACCACCAAATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCaa                                                                              <  1:392214/140‑1 (MQ=255)
                                                              cTGGTCCATATGTACCACCATATTGAACATGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCaa                                                                              <  2:82919/141‑1 (MQ=255)
                                                                 gTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAat                                                                           >  2:95874/1‑141 (MQ=255)
                                                                  tCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCTTCAAAata                                                                          >  1:621440/1‑141 (MQ=255)
                                                                  tCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATAAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAata                                                                          >  1:531357/1‑141 (MQ=255)
                                                                          tACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGaa                                                                  >  1:769074/1‑141 (MQ=255)
                                                                             caccaTATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTg                                                               <  1:737907/141‑1 (MQ=255)
                                                                                caTATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTc                                                            >  2:452435/1‑141 (MQ=255)
                                                                                 atatTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTTTAAATTTTCATTTAATTTAGAAGATGTTATATCTTCTTCATCAAAATCTATTAGAACTGCTCg                                                           >  1:503792/1‑141 (MQ=255)
                                                                                  tatTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGa                                                          >  2:336357/1‑141 (MQ=255)
                                                                                   tgtgaggcgagttttcAATAGTACGAGTGGTTGGATTCACATGATGAACTTTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGGTCATCAAAATATATTAGAACTGCTCGAt                                                         <  2:666188/129‑1 (MQ=255)
                                                                                   atTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGAt                                                         <  1:606854/141‑1 (MQ=255)
                                                                                      gAACAGGTATTTTAATTATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACg                                                      <  1:26180/141‑1 (MQ=255)
                                                                                      gAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACg                                                      <  1:699880/141‑1 (MQ=255)
                                                                                        aCAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGaa                                                    >  1:381352/1‑141 (MQ=255)
                                                                                           ggTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTTTTTTTTTTTGGTAAGATCTAAATTTTCAATTAATTTGGAAGATGTGATAGCTTGTTCATCAAAATTTCTTACAACTGTTCGATACGAAggt                                                 >  2:227679/1‑140 (MQ=255)
                                                                                            gTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGa                                                <  1:603321/141‑1 (MQ=255)
                                                                                            gTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGAc                                               >  2:417277/1‑141 (MQ=255)
                                                                                              aTTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACa                                              >  2:48154/1‑141 (MQ=255)
                                                                                               ttttAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACaa                                             >  1:445016/1‑141 (MQ=255)
                                                                                               ttttAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACaa                                             >  1:281974/1‑141 (MQ=255)
                                                                                               ttttAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACaa                                             >  1:388768/1‑141 (MQ=255)
                                                                                                  taagtagtcGATTGGTTTGATTCACATCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAatat                                          <  2:335862/133‑1 (MQ=255)
                                                                                                        taCGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGTTACGAAGGGACAATATCTTGaa                                    >  2:470081/1‑141 (MQ=255)
                                                                                                         aCGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAAt                                   <  2:210836/141‑1 (MQ=255)
                                                                                                          cGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAAtt                                  <  2:340435/141‑1 (MQ=255)
                                                                                                           gATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACGGCTCGATACGAAGGGACAATATCTTGAAttt                                 <  2:190592/141‑1 (MQ=255)
                                                                                                            aTTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTc                                <  1:336357/141‑1 (MQ=255)
                                                                                                            aTTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTc                                <  2:320172/141‑1 (MQ=255)
                                                                                                               ggTTTGATTCACAGCATGAACATTTTTTTTTTTTTTTTTTTTAAGTTATAAGGTTTAAATTATTTTATAAGGTGTTATAGGTTTTTCATCAAAATATATTAGAACTGATCAATATGAATGGACAATATCTTTAATTTCTat                             >  2:6634/1‑140 (MQ=255)
                                                                                                               ggTTTGATTCACAGCATGAACATTTTTTTTTTTTATTTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTaa                             >  2:652239/1‑141 (MQ=255)
                                                                                                               ggTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAGATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATATTGAATTTCTaa                             >  2:665647/1‑141 (MQ=255)
                                                                                                                 tttGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAata                           <  2:171935/141‑1 (MQ=255)
                                                                                                                   tGATTCACAGCATGAACATTTTTTTTTTTTTTTTTTTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGTAGGGAGAGTATCTTTATTTTCTAAtttt                         >  1:684405/1‑141 (MQ=255)
                                                                                                                   tGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATAtt                         <  1:518423/141‑1 (MQ=255)
                                                                                                                    gATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATAttt                        <  1:624964/141‑1 (MQ=255)
                                                                                                                       ttataccATGAACATTTTTTTTATTGAATTTGTTAAGTTCTCAATTCTCCATTAATTTCGAAGATGTTATAGCTTGTCCTCTAAAATATATTAGGACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAAtctttctt                     <  2:400287/135‑1 (MQ=255)
                                                                                                                             cATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATATTTCTTTTtctc               <  1:137341/141‑1 (MQ=255)
                                                                                                                                   cttttttttCATTAAATTTGTTAAGTTCTATATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATCTATTCTATTTCTAATATTTTTTTTTCTCGTATcc         <  2:328500/139‑1 (MQ=255)
                                                                                                                                     tttttgtCCTTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATACTTCTTTTTCTCGTATCCAc       <  1:192332/141‑1 (MQ=255)
                                                                                                                                         tttCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATATTTCTTTTTCTCGTATCCACCGTa   >  1:582178/1‑141 (MQ=255)
                                                                                                                                          ttCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATATTTCTTTTTCTCGTATCCACCGTAc  <  1:95874/141‑1 (MQ=255)
                                                                                                                                          ttCATTAAATTTGTTAAGTTCTAAACTTTCAATTAATTTAGAAGATGTTATATCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATATTTCTTTTTCTCGTATCCACCGTAc  <  1:101454/141‑1 (MQ=255)
                                                                                                                                          ttCATTAAATTTGTTAAGCTCTAAATTTTCAATTAATTTAGAATATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACCAAGGTACAATATCCTGAATTTCTCATATTTCTTTTTCTCGTATCCACCGTAc  <  2:143173/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TGTATGTTTTTCAATAACTTGCTCAACAGTTATTCGATTATGCTTTGCTACTTCTTCAATATCTGGTCCATATGTACCACCATATTGAACAGGTATTTTAATAATACGATTGGTTTGATTCACAGCATGAACATTTTTTTCATTAAATTTGTTAAGTTCTAAATTTTCAATTAATTTAGAAGATGTTATAGCTTGTTCATCAAAATATATTAGAACTGCTCGATACGAAGGGACAATATCTTGAATTTCTAATATTTCTTTTTCTCGTATCCACCGTAC  >  CP000730/1716618‑1716896

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