Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,335,195 G→T 5.5% L61I (CTC→ATC)  yjdF ← DUF2238 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,335,1950GT 5.5% 822.4 / 15.1 215L61I (CTC→ATC) yjdFDUF2238 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (112/91);  new base T (6/6);  total (118/97)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 225 total aligned reads displayed.

CCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGT  >  CP009273/4335049‑4335342
                                                                                                                                                  |                                                                                                                                                   
ccccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGTGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagc                                                                                                                                                    <  2:517703/148‑1 (MQ=255)
ccccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagc                                                                                                                                                    >  2:2378541/1‑148 (MQ=255)
ccccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagc                                                                                                                                                    >  2:2040426/1‑148 (MQ=255)
  ccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagcag                                                                                                                                                  >  1:2515751/1‑148 (MQ=255)
  ccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagcag                                                                                                                                                  <  2:1925793/148‑1 (MQ=255)
     ggAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAgcgg                                                                                                                                               >  2:1329684/1‑148 (MQ=255)
       aaaaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAgcggcg                                                                                                                                             <  1:3230016/148‑1 (MQ=255)
       aaaaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAgcggcg                                                                                                                                             <  1:1644278/148‑1 (MQ=255)
        aaaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGt                                                                                                                                            >  1:1342541/1‑148 (MQ=255)
         aaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGtt                                                                                                                                           <  2:1112366/148‑1 (MQ=255)
         aaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGtt                                                                                                                                           >  1:813540/1‑148 (MQ=255)
                 ccAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATa                                                                                                                                   >  1:3319649/1‑148 (MQ=255)
                 ccAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATa                                                                                                                                   >  2:1148152/1‑148 (MQ=255)
                     cTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGt                                                                                                                               >  2:301790/1‑148 (MQ=255)
                      ttGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTc                                                                                                                              <  2:1267330/148‑1 (MQ=255)
                       tGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCt                                                                                                                             >  2:3315822/1‑148 (MQ=255)
                         tCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTgg                                                                                                                           >  1:3374230/1‑148 (MQ=255)
                            tACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCgg                                                                                                                        >  1:586869/1‑148 (MQ=255)
                            tACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCgg                                                                                                                        >  1:2084999/1‑148 (MQ=255)
                                 aTTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCa                                                                                                                   >  2:1843350/1‑148 (MQ=255)
                                 aTTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCa                                                                                                                   >  2:1525994/1‑148 (MQ=255)
                                   tACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAg                                                                                                                 >  1:565606/1‑148 (MQ=255)
                                   tACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAg                                                                                                                 <  2:1484535/148‑1 (MQ=255)
                                    aCGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGc                                                                                                                >  2:1070940/1‑148 (MQ=255)
                                    aCGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGc                                                                                                                >  1:2703912/1‑148 (MQ=255)
                                     cGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCa                                                                                                               >  2:2897284/1‑148 (MQ=255)
                                     cGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCa                                                                                                               <  2:1144876/148‑1 (MQ=255)
                                      gACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACAATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAg                                                                                                              >  2:1911118/1‑148 (MQ=255)
                                       aCTCAACCCCAGCCATTCCTGCACCTCAAACCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGt                                                                                                             <  1:532147/148‑1 (MQ=255)
                                       aCTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGt                                                                                                             >  2:2809830/1‑148 (MQ=255)
                                                   ccATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGGGGCAAGCAGTAGCTGCACGACa                                                                                                 >  1:1517974/1‑148 (MQ=255)
                                                    cATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACaa                                                                                                <  2:2808516/148‑1 (MQ=255)
                                                      ttCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACaata                                                                                              <  2:3436715/148‑1 (MQ=255)
                                                       tCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACaataa                                                                                             >  1:2781538/1‑148 (MQ=255)
                                                          tGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCa                                                                                          <  1:3519766/148‑1 (MQ=255)
                                                          tGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCa                                                                                          <  1:2444465/148‑1 (MQ=255)
                                                           gCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCAc                                                                                         >  1:3477344/1‑148 (MQ=255)
                                                               cTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgt                                                                                     >  2:510633/1‑148 (MQ=255)
                                                               cTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgt                                                                                     >  2:642103/1‑148 (MQ=255)
                                                               cTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgt                                                                                     >  2:1511570/1‑148 (MQ=255)
                                                                   aaaCCACCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCa                                                                                 <  2:284649/148‑1 (MQ=255)
                                                                   aaaCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCa                                                                                 <  2:3374230/148‑1 (MQ=255)
                                                                    aaCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCAc                                                                                >  2:1118310/1‑148 (MQ=255)
                                                                    aaCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCAc                                                                                <  2:2252118/148‑1 (MQ=255)
                                                                      ccAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACtt                                                                              >  2:845994/1‑148 (MQ=255)
                                                                       cAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTc                                                                             <  1:1025622/148‑1 (MQ=255)
                                                                        aaCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTcc                                                                            <  1:535053/148‑1 (MQ=255)
                                                                          cGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCAt                                                                          <  2:664690/148‑1 (MQ=255)
                                                                            gggACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCa                                                                        >  1:1300356/1‑148 (MQ=255)
                                                                             ggACTTTCGCGTAGGTGTATTGTCCGCCGTCGATGATTATGATGGCGTGAAGGAATATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGGGTCACTTccttaag                                                                       <  2:726061/148‑6 (MQ=255)
                                                                             ggACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAg                                                                       >  1:780855/1‑148 (MQ=255)
                                                                             ggACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAg                                                                       >  2:1858298/1‑148 (MQ=255)
                                                                                 tttCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAg                                                                   >  1:1791503/1‑148 (MQ=255)
                                                                                  ttCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAgg                                                                  >  2:3306355/1‑148 (MQ=255)
                                                                                      cgTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAg                                                              <  1:301939/148‑1 (MQ=255)
                                                                                        tAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGtt                                                            <  2:1668027/148‑1 (MQ=255)
                                                                                        tAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGtt                                                            <  1:1078447/148‑1 (MQ=255)
                                                                                        tAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGtt                                                            <  1:2859030/148‑1 (MQ=255)
                                                                                         aGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCTGCGTTAACGGATAACGTCTGGCTGTGGCAAGCAGTAGCTGCACGACAATCATCACCGGTGTCACTTCCATCAGCCAGGTGAGttt                                                           <  2:1533453/148‑1 (MQ=255)
                                                                                         aGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGttt                                                           >  1:1631212/1‑148 (MQ=255)
                                                                                          ggTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTg                                                          >  1:2590095/1‑148 (MQ=255)
                                                                                          ggTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTg                                                          >  1:2098346/1‑148 (MQ=255)
                                                                                             gtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCa                                                       <  2:736782/148‑1 (MQ=255)
                                                                                              tATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCAt                                                      >  2:3538884/1‑148 (MQ=255)
                                                                                               aTTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATg                                                     >  1:856444/1‑148 (MQ=255)
                                                                                                 tGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATggg                                                   <  2:2462059/148‑1 (MQ=255)
                                                                                                  gTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGc                                                  <  1:2393877/148‑1 (MQ=255)
                                                                                                   tCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGcc                                                 >  1:1131539/1‑148 (MQ=255)
                                                                                                      gccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTGGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGaa                                              <  1:1223254/148‑1 (MQ=255)
                                                                                                      gccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGaa                                              <  1:2243558/148‑1 (MQ=255)
                                                                                                        cgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAAt                                            >  1:1959587/1‑148 (MQ=255)
                                                                                                         gACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATg                                           <  1:2313831/148‑1 (MQ=255)
                                                                                                         gACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATg                                           <  2:2084999/148‑1 (MQ=255)
                                                                                                            cATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCg                                        <  1:132054/148‑1 (MQ=255)
                                                                                                             aTCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCgg                                       >  2:2220210/1‑148 (MQ=255)
                                                                                                              tCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCGGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGt                                      <  1:2452746/148‑1 (MQ=255)
                                                                                                                     gatgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATc                               >  2:3142902/1‑148 (MQ=255)
                                                                                                                     gatgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATc                               >  1:3318907/1‑148 (MQ=255)
                                                                                                                      atgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCa                              >  2:1598576/1‑148 (MQ=255)
                                                                                                                       tgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAg                             >  2:2070448/1‑148 (MQ=255)
                                                                                                                        ggCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTAATGGGCCGAAATGCCGGTATAAATCAgg                            <  1:1762836/148‑1 (MQ=255)
                                                                                                                            tGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCTGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTa                        <  2:3043338/148‑1 (MQ=255)
                                                                                                                             gAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTaa                       >  2:2721482/1‑148 (MQ=255)
                                                                                                                              aaGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTGATGGGCCGAAATGCCGGTATAAATCAGGATTAAc                      >  1:3522726/1‑148 (MQ=255)
                                                                                                                               aGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACg                     <  2:1855428/148‑1 (MQ=255)
                                                                                                                                ggAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGt                    <  1:430724/148‑1 (MQ=255)
                                                                                                                                   aaaTGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAg                 <  2:862355/148‑1 (MQ=255)
                                                                                                                                    aaTGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGc                >  2:3120344/1‑148 (MQ=255)
                                                                                                                                    aaTGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGc                >  1:2954140/1‑148 (MQ=255)
                                                                                                                                      tGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGt              >  1:2374027/1‑148 (MQ=255)
                                                                                                                                       gAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTc             <  1:2037501/148‑1 (MQ=255)
                                                                                                                                         gCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAg           <  2:419728/148‑1 (MQ=255)
                                                                                                                                         gCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAg           >  1:695375/1‑148 (MQ=255)
                                                                                                                                         gCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAg           <  2:788338/148‑1 (MQ=255)
                                                                                                                                            tataGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgc        >  1:2644505/1‑148 (MQ=255)
                                                                                                                                             ataGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcg       >  2:571260/1‑148 (MQ=255)
                                                                                                                                             ataGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcg       <  2:2619834/148‑1 (MQ=255)
                                                                                                                                              taGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcgc      >  1:2539264/1‑148 (MQ=255)
                                                                                                                                                gagCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTg    >  2:2052660/1‑148 (MQ=255)
                                                                                                                                                  gcagcGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGgt  <  1:311554/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGT  >  CP009273/4335049‑4335342

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