Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 684,444 G→C T9R (ACA→AGA)  insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273684,4440GC98.5% 193.3 / NA 67T9R (ACA→AGA) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base C (32/33);  minor base A (1/0);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAA  >  CP009273/684361‑684589
                                                                                   |                                                                                                                                                 
aggaaaGCTCTTGTCGTGTCTGACGGCGCTGACGTCTCAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTaa                                                                                  <  2:61736/146‑1 (MQ=11)
aGAAAATCTCTTTTCTGGTGTGGAGGAGTTTATTGCTTATTTGACTGTCGGCGGAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTaa                                                                                  >  2:626735/1‑149 (MQ=11)
aGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTaa                                                                                  <  2:1389196/149‑1 (MQ=14)
    actctcTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTTATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAtt                                                                              <  1:833379/147‑1 (MQ=16)
    aaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAtt                                                                              <  1:583365/149‑1 (MQ=21)
    aaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATCGTATTTGTTCTTGGACTTGTTCGCACCTTCCTTAGGTAACAtt                                                                              <  1:982339/149‑1 (MQ=11)
    aaTCTCTTTTCTGGGCTGACGGCGCTTACGGCTGAATTCCCTTTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAtt                                                                              <  1:519375/149‑1 (MQ=11)
     aTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAttt                                                                             <  1:1062069/149‑1 (MQ=25)
         ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGtt                                                                         >  1:19359/1‑149 (MQ=38)
          ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGttt                                                                        >  1:1525030/1‑149 (MQ=255)
          ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGAATCATGCTGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGttt                                                                        >  2:1688953/1‑149 (MQ=255)
           tttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTg                                                                       >  2:115043/1‑149 (MQ=255)
           tttCTGGTCTGACGGAGCTTCCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTg                                                                       <  1:936008/149‑1 (MQ=255)
            ttCCGGTCTGACGGCGCTTACGGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGCTTTGGTCGCACCTTCCTTAGGTAACATTTAGTTTgg                                                                      <  2:402194/149‑1 (MQ=14)
             tCTGGTCTGACGGCGCTTACTGCTGAATTCACGGGCCGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGc                                                                     <  1:269974/149‑1 (MQ=255)
              cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGCCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCt                                                                    <  2:1694202/149‑1 (MQ=255)
              cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCt                                                                    <  1:106532/149‑1 (MQ=255)
                 gTCTGACGGCGCTGACGGCTGAATTCACTGTCGGCGAATGTAATTTGATGACTCATGATGAACCCTCTTCGATGGCTCCAGATGACAAACATGAGCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTaaa                                                                 <  1:416337/149‑1 (MQ=14)
                  tCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAAt                                                                <  1:986158/149‑1 (MQ=255)
                    tGACGGCGCTTACTGCTTAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGATCGCCCCCGCCTTAGGTAACATTTAGTTTGGCTAAATGt                                                              <  2:402625/149‑1 (MQ=14)
                      aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTACATGTaa                                                            >  1:950614/1‑149 (MQ=255)
                       cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGACTCATGATGAACCATCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTACAATTTAGTTTGGCTCAATGTaac                                                           >  2:100612/1‑148 (MQ=14)
                       cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCATATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTaaa                                                           <  1:1769878/149‑1 (MQ=255)
                       cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTaaa                                                           >  2:32201/1‑149 (MQ=255)
                        ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGATGATGACTCATGATTCACCCTATTCTATGGTTTCAGATGACAAACATGATCTCATATCAAGAGCTTGTTCGCACCTTCCTTCGGTGACATTTCGTTTGGCTAAATGTAAAg                                                          >  2:115082/1‑149 (MQ=14)
                         gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTGGGTAACATTTAGTTTGGCTAAAGGTAAAgg                                                         >  2:943822/1‑148 (MQ=255)
                             ttACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGACAtt                                                     >  2:1745124/1‑149 (MQ=255)
                                 tgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTg                                                 >  2:1778882/1‑149 (MQ=255)
                                 tgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTg                                                 <  1:666260/149‑1 (MQ=255)
                                 tgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTg                                                 >  2:691474/1‑149 (MQ=255)
                                   ctgAATTCACTGTCGGCGAATGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCATATTACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGtt                                               <  1:1735026/149‑1 (MQ=255)
                                   ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGtt                                               <  1:482948/149‑1 (MQ=255)
                                   ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCGTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGtt                                               <  1:895052/149‑1 (MQ=255)
                                     gAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCGCCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGtttt                                             >  2:645866/1‑149 (MQ=255)
                                       aTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAATCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAt                                           >  2:1298587/1‑149 (MQ=255)
                                        ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTtttttt                                          >  2:979524/1‑149 (MQ=255)
                                        ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAtt                                          <  1:672273/149‑1 (MQ=255)
                                        ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCGAAATGTAAAGATATTGCTGTTTTAtt                                          >  1:433898/1‑149 (MQ=255)
                                        ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATCTCCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAtt                                          >  1:386606/1‑149 (MQ=255)
                                          cACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgt                                        >  1:939294/1‑149 (MQ=255)
                                          cACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATATCATATCCGGGACTTGGTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGAGAGTGCTGTTTTAttgt                                        >  1:781832/1‑149 (MQ=18)
                                            cTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgttt                                      <  1:1234100/149‑1 (MQ=255)
                                            cTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgttt                                      <  1:1958622/149‑1 (MQ=255)
                                            cTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCGATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTCAATGTAAAGATATTGCTGtgtgatgggtt                                      <  2:1223157/149‑11 (MQ=17)
                                            cTGTCGGCGAAGGGAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgttt                                      <  2:531787/149‑1 (MQ=255)
                                               tCGGCGAAGGTAAGTTGATGACTCATGATTAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCGGTTTTATtgtttgtt                                   >  1:1223138/1‑149 (MQ=255)
                                                 ggCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAACGGTATTGCTGTTTTATTGTTTGtttt                                 >  1:495085/1‑149 (MQ=255)
                                                  gCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGttttt                                >  1:1517659/1‑149 (MQ=255)
                                                  gCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAAAATGATCTCATATCAGGGACTTGTTCGCAACTTAATTAGGTAAAATTTAGTTTGGGTAAATGTAAAGATATTTCTGTTTTATTGTTTGttttt                                >  1:1708877/1‑149 (MQ=18)
                                                   cgacgctAAGTTGATGACTCATGATGACCCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTg                               <  1:150404/143‑1 (MQ=255)
                                                   cGAGGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTg                               <  2:522430/149‑1 (MQ=255)
                                                     aaGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCg                             >  1:1913187/1‑149 (MQ=255)
                                                      aGGTAAGGGGAAGACTCTAGATGCACCCTCTTCTATGTCTCCAGATGCCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAATGTTTGTTTTTGCga                            <  1:366575/149‑1 (MQ=21)
                                                        gTAAGTTGACGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTGAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCgaga                          <  1:440990/149‑1 (MQ=255)
                                                        gTAAGGTTATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCgaga                          <  2:1251354/149‑1 (MQ=255)
                                                               gatgaCTCATGATGAACCCTCTTCTGTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGcgccg                   <  1:1233897/149‑1 (MQ=255)
                                                               gatgaCTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTCAAGATATTGCTGTTTTATTGTTTGTTTTTTCGCGCTGcgccc                   >  2:841149/1‑148 (MQ=21)
                                                               gatgaCTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGcgccg                   >  2:278439/1‑149 (MQ=255)
                                                                     tCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGCCTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCAt             >  1:539697/1‑149 (MQ=255)
                                                                        tgatgaATCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTcc          <  1:1411822/149‑1 (MQ=255)
                                                                        tgatgaACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGATATGCGCCGCACCATTcc          >  2:1033789/1‑149 (MQ=255)
                                                                          atgaACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGa        >  2:344507/1‑149 (MQ=255)
                                                                           ggatcccTCTTCTATGTCTCCAGATGACAAACATGATCTCATAGCAGGGACTTGTTCCCACCTTCCTTAGGTAACATTTAGTTTGGTTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGaa       <  2:748681/145‑1 (MQ=37)
                                                                            gaACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTGTTTGCGAGATGCGCCGCACCATTCCGAAg      >  2:1465341/1‑149 (MQ=255)
                                                                                ccTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGCCTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCaaa  >  2:63285/1‑149 (MQ=255)
                                                                                ccTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCaaa  >  1:1834819/1‑149 (MQ=255)
                                                                                ccTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCaaa  >  2:375567/1‑149 (MQ=255)
                                                                                   |                                                                                                                                                 
AGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAAA  >  CP009273/684361‑684589

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: