Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 580,366 G→T 100% R242S (CGC→AGC)  ompT ← DLP12 prophage; outer membrane protease VII (outer membrane protein 3b)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,3660GT96.3% 836.5 / 20.0 273R242S (CGC→AGC) ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base T (133/130);  total (137/136)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 274 total aligned reads displayed.

TATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTG  >  CP009273/580222‑580510
                                                                                                                                                |                                                                                                                                                
   tgtgATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGt                                                                                                                                            >  1:2901536/1‑148 (MQ=255)
     tgATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGa                                                                                                                                          <  1:1219419/148‑1 (MQ=255)
       aTCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGAtt                                                                                                                                        <  1:1383378/148‑1 (MQ=255)
        tCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTc                                                                                                                                       <  1:1882680/148‑1 (MQ=255)
         cATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCt                                                                                                                                      >  1:2324231/1‑148 (MQ=255)
            aaaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtttt                                                                                                                                   <  2:1993184/148‑1 (MQ=255)
            aaaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtttt                                                                                                                                   >  1:3452051/1‑148 (MQ=255)
              aGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                 >  1:1005234/1‑148 (MQ=255)
              aGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                 >  1:705041/1‑148 (MQ=255)
                tGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTccc                                                                                                                               <  2:151061/148‑1 (MQ=255)
                 gAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCg                                                                                                                              <  1:1601570/148‑1 (MQ=255)
                   aGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCggg                                                                                                                            <  2:1822233/148‑1 (MQ=255)
                     tATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTc                                                                                                                          >  1:968356/1‑148 (MQ=255)
                      aTTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCa                                                                                                                         <  1:2169710/148‑1 (MQ=255)
                      aTTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCa                                                                                                                         >  1:2423396/1‑148 (MQ=255)
                        tACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATa                                                                                                                       >  1:3707478/1‑148 (MQ=255)
                         aCCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAg                                                                                                                      >  2:3692036/1‑148 (MQ=255)
                         aCCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAg                                                                                                                      >  2:3598134/1‑148 (MQ=255)
                          ccTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgt                                                                                                                     >  2:683100/1‑148 (MQ=255)
                          ccTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAgt                                                                                                                     >  1:3734512/1‑148 (MQ=255)
                           cTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtg                                                                                                                    <  2:3144848/148‑1 (MQ=255)
                            ttttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtgt                                                                                                                   >  2:3280489/1‑148 (MQ=255)
                             tttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGtt                                                                                                                  <  1:136202/148‑1 (MQ=255)
                              ttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTc                                                                                                                 >  1:2704295/1‑148 (MQ=255)
                              ttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTc                                                                                                                 >  1:3749026/1‑148 (MQ=255)
                                ttattaGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                               >  2:779378/1‑148 (MQ=255)
                                   ttCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGt                                                                                                            <  2:1275618/148‑1 (MQ=255)
                                   ttCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGt                                                                                                            <  1:1996196/148‑1 (MQ=255)
                                    tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           <  1:1055867/148‑1 (MQ=255)
                                    tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           <  1:472774/148‑1 (MQ=255)
                                    tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           >  1:3491603/1‑148 (MQ=255)
                                         aCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAg                                                                                                      >  2:2915924/1‑148 (MQ=255)
                                         aCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAg                                                                                                      >  2:45165/1‑148 (MQ=255)
                                          cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAga                                                                                                     >  2:1802669/1‑148 (MQ=255)
                                            cGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                   >  2:3161834/1‑148 (MQ=255)
                                             gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatga                                                                                                  <  2:944671/148‑1 (MQ=255)
                                             gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatga                                                                                                  <  1:1138700/148‑1 (MQ=255)
                                              aTTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatgat                                                                                                 <  1:1507445/148‑1 (MQ=255)
                                              aTTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatgat                                                                                                 >  2:2378889/1‑148 (MQ=255)
                                                    tGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAcc                                                                                           <  1:3261465/148‑1 (MQ=255)
                                                      cgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCa                                                                                         <  1:3671003/148‑1 (MQ=255)
                                                          ttCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCg                                                                                     >  1:3512312/1‑148 (MQ=255)
                                                            cAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCt                                                                                   <  2:194675/148‑1 (MQ=255)
                                                               cATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTa                                                                                >  1:413922/1‑148 (MQ=255)
                                                                     cTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                          <  2:3620188/148‑1 (MQ=255)
                                                                      ttttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAt                                                                         <  1:57397/148‑1 (MQ=255)
                                                                      ttttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGCGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAt                                                                         <  1:1970063/148‑1 (MQ=255)
                                                                         tGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtgtg                                                                      <  1:3027044/148‑1 (MQ=255)
                                                                           cGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGcc                                                                    >  2:2619108/1‑148 (MQ=255)
                                                                           cGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGcc                                                                    >  2:1021671/1‑148 (MQ=255)
                                                                            gTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCa                                                                   <  2:297499/148‑1 (MQ=255)
                                                                                 gtgtGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAg                                                              <  1:3271982/148‑1 (MQ=255)
                                                                                 gtgtGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAg                                                              <  1:1314092/148‑1 (MQ=255)
                                                                                   gtgACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATTTGCCACCGAGtt                                                            >  1:1003786/1‑148 (MQ=255)
                                                                                       cGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaaa                                                        <  2:3052833/148‑1 (MQ=255)
                                                                                       cGTAATAACCTGCATTGACCGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaaa                                                        >  2:1262547/1‑148 (MQ=255)
                                                                                         taataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCa                                                                                                                <  1:14128/90‑1 (MQ=255)
                                                                                         taataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAAt                                                      <  1:3569018/148‑1 (MQ=255)
                                                                                          aataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATc                                                     >  1:709738/1‑148 (MQ=255)
                                                                                          aataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATc                                                     <  2:3749026/148‑1 (MQ=255)
                                                                                          aataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATc                                                     <  2:2704295/148‑1 (MQ=255)
                                                                                           ataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCt                                                    <  1:1685104/148‑1 (MQ=255)
                                                                                           ataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCt                                                    >  2:1426558/1‑148 (MQ=255)
                                                                                           ataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCt                                                    >  1:182473/1‑148 (MQ=255)
                                                                                             aaCCTGCATTGACTGCAACAGAATATTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTc                                                  <  2:2293673/148‑1 (MQ=255)
                                                                                              aCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCa                                                 <  1:201497/148‑1 (MQ=255)
                                                                                              aCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCa                                                 >  2:2198613/1‑148 (MQ=255)
                                                                                              aCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCa                                                 >  2:271008/1‑148 (MQ=255)
                                                                                                  gCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAAc                                             >  1:1444739/1‑148 (MQ=255)
                                                                                                   cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACg                                            <  2:1868393/148‑1 (MQ=255)
                                                                                                   cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACg                                            <  1:3220775/148‑1 (MQ=255)
                                                                                                    aTTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGa                                           <  1:3123283/148‑1 (MQ=255)
                                                                                                       gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATaa                                        <  1:3418540/148‑1 (MQ=255)
                                                                                                       gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATaa                                        <  2:3247629/148‑1 (MQ=255)
                                                                                                       gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATaa                                        <  1:2869474/148‑1 (MQ=255)
                                                                                                           gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATATTAATAACGATAACTTc                                    >  1:1534851/1‑148 (MQ=255)
                                                                                                             aaCAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCa                                  <  1:2727958/148‑1 (MQ=255)
                                                                                                                 gAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCa                              <  2:3437876/148‑1 (MQ=255)
                                                                                                                    tAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGc                           >  2:613049/1‑148 (MQ=255)
                                                                                                                       tAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCaa                        <  1:873097/148‑1 (MQ=255)
                                                                                                                       tAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCaa                        <  2:1734813/148‑1 (MQ=255)
                                                                                                                       tAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATATTCATAACGATAACTTCCAGTCAAGCCaa                        >  2:1369834/1‑148 (MQ=255)
                                                                                                                        aaTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCGGTATTTAAATGTGCCACCGAGTTCAAAATCTTAATAACGATAAATTACAGTCAAGCCaaa                       >  1:2256760/1‑147 (MQ=255)
                                                                                                                         aTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATg                      <  2:2423396/148‑1 (MQ=255)
                                                                                                                          ttttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGt                     <  2:3029976/148‑1 (MQ=255)
                                                                                                                           tttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAAATTCCAGTCAAGCCAATGTa                    >  1:2963022/1‑148 (MQ=255)
                                                                                                                            ttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAg                   <  1:667766/148‑1 (MQ=255)
                                                                                                                            ttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGGCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAg                   <  2:3193210/148‑1 (MQ=255)
                                                                                                                               gTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGc                >  2:1169162/1‑148 (MQ=255)
                                                                                                                                tCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCa               <  2:530530/148‑1 (MQ=255)
                                                                                                                                  tttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATg                      >  2:1528128/1‑139 (MQ=255)
                                                                                                                                   ttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt            >  1:793143/1‑148 (MQ=255)
                                                                                                                                   ttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt            >  2:2526008/1‑148 (MQ=255)
                                                                                                                                   ttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt            >  1:2478486/1‑148 (MQ=255)
                                                                                                                                   ttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt            <  1:3280489/148‑1 (MQ=255)
                                                                                                                                   ttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt            >  1:2468099/1‑148 (MQ=255)
                                                                                                                                      aCCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTaa         >  2:3790734/1‑148 (MQ=255)
                                                                                                                                       ccTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTaaa        <  1:1785233/148‑1 (MQ=255)
                                                                                                                                          tACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg     >  2:1599440/1‑148 (MQ=255)
                                                                                                                                           aCTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgt    <  2:633251/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTG  >  CP009273/580222‑580510

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: