Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 580,366 G→T 85.8% 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,3660GT85.8% 415.4 / 71.6 183R242S (CGC→AGC) ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (12/14);  new base T (83/74);  total (95/88)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 185 total aligned reads displayed.

ATTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTG  >  CP009273/580220‑580510
                                                                                                                                                  |                                                                                                                                                
attattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGa                                                                                                                                                 >  2:1457053/1‑148 (MQ=255)
  tattGTGACCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCtata                                                                                                                                               >  2:1099868/1‑148 (MQ=255)
   attGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATaa                                                                                                                                              <  2:2172089/148‑1 (MQ=255)
    ttGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAg                                                                                                                                             <  2:39645/148‑1 (MQ=255)
     tgtgATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGt                                                                                                                                            >  1:895049/1‑148 (MQ=255)
     tgtgATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGt                                                                                                                                            <  1:17487/148‑1 (MQ=255)
         aTCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGAtt                                                                                                                                        >  2:967220/1‑148 (MQ=255)
           cATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCt                                                                                                                                      <  2:264850/148‑1 (MQ=255)
            aTAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtt                                                                                                                                     <  1:1022976/148‑1 (MQ=255)
            aTAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtt                                                                                                                                     <  1:1022945/148‑1 (MQ=255)
             tAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCttt                                                                                                                                    >  1:772588/1‑148 (MQ=255)
                aGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                 >  1:747844/1‑148 (MQ=255)
                aGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACCTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                 <  1:1951964/148‑1 (MQ=255)
                 gTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTcc                                                                                                                                <  1:1497216/148‑1 (MQ=255)
                    aaGTATTAACTTTTTTATTCGTAACCCGATTCAATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGAAGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGATATAAGTGATACTTTTTccagt                                                                                                                             >  2:722795/1‑145 (MQ=255)
                      gTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGt                                                                                                                           >  1:2545586/1‑148 (MQ=255)
                      gTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGt                                                                                                                           >  2:863798/1‑148 (MQ=255)
                        aTTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCa                                                                                                                         >  2:340474/1‑148 (MQ=255)
                         ttACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCAt                                                                                                                        >  2:1229249/1‑148 (MQ=255)
                          tACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATa                                                                                                                       >  1:146430/1‑148 (MQ=255)
                          tACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATa                                                                                                                       >  2:2354744/1‑148 (MQ=255)
                              ttttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtgt                                                                                                                   >  2:1369622/1‑148 (MQ=255)
                              ttttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtgt                                                                                                                   >  1:1203886/1‑148 (MQ=255)
                              ttttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtgt                                                                                                                   >  1:2858618/1‑148 (MQ=255)
                                 tttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCa                                                                                                                >  1:849016/1‑148 (MQ=255)
                                  ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCTGGTCATAGTGTTCAt                                                                                                               >  2:892542/1‑148 (MQ=255)
                                  ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                               <  2:3044077/148‑1 (MQ=255)
                                  ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                               >  1:3284636/1‑148 (MQ=255)
                                      tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           <  2:1538261/148‑1 (MQ=255)
                                      tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           >  2:270452/1‑148 (MQ=255)
                                      tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           <  2:2859767/148‑1 (MQ=255)
                                         tAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATc                                                                                                        <  1:2603191/148‑1 (MQ=255)
                                           aCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAg                                                                                                      <  2:1126725/148‑1 (MQ=255)
                                           aCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAg                                                                                                      <  2:3268711/148‑1 (MQ=255)
                                            cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAga                                                                                                     <  2:133093/148‑1 (MQ=255)
                                            cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAga                                                                                                     >  1:572139/1‑148 (MQ=255)
                                               gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatga                                                                                                  <  1:2239489/148‑1 (MQ=255)
                                                 ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                >  2:1463444/1‑148 (MQ=255)
                                                 ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                >  2:839637/1‑148 (MQ=255)
                                                 ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                >  2:839586/1‑148 (MQ=255)
                                                 ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                >  1:3141080/1‑148 (MQ=255)
                                                 ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                >  1:3141040/1‑148 (MQ=255)
                                                     aTGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAc                                                                                            >  1:3175768/1‑148 (MQ=255)
                                                       gcgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAccc                                                                                          >  2:1525587/1‑148 (MQ=255)
                                                        cgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCa                                                                                         >  2:2397984/1‑148 (MQ=255)
                                                            ttCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCg                                                                                     <  1:2952676/148‑1 (MQ=255)
                                                             tCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGc                                                                                    >  1:479993/1‑148 (MQ=255)
                                                             tCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGc                                                                                    >  2:2697680/1‑148 (MQ=255)
                                                              cAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCt                                                                                   >  2:2152920/1‑148 (MQ=255)
                                                                 cATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTa                                                                                >  1:235619/1‑148 (MQ=255)
                                                                    aaaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAttt                                                                             <  2:2514239/148‑1 (MQ=255)
                                                                      aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                           <  2:953024/148‑1 (MQ=255)
                                                                       cTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                          <  1:1219079/148‑1 (MQ=255)
                                                                          ttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtgt                                                                       <  1:86193/148‑1 (MQ=255)
                                                                              gTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCa                                                                   <  2:572139/148‑1 (MQ=255)
                                                                               ttAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCAc                                                                  >  1:1951212/1‑148 (MQ=255)
                                                                                tAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCAcc                                                                 >  2:3191719/1‑148 (MQ=255)
                                                                                 aGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCg                                                                <  2:3088079/148‑1 (MQ=255)
                                                                                        aCGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaa                                                         >  1:2714929/1‑148 (MQ=255)
                                                                                              taaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCtt                                                   <  1:3047104/148‑1 (MQ=255)
                                                                                                 ccTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCAt                                                <  1:1565896/148‑1 (MQ=255)
                                                                                                 ccTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCAt                                                >  1:395426/1‑148 (MQ=255)
                                                                                                     cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACg                                            >  1:207210/1‑148 (MQ=255)
                                                                                                     cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACg                                            <  2:747844/148‑1 (MQ=255)
                                                                                                       ttGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGAt                                          >  2:2013483/1‑148 (MQ=255)
                                                                                                       ttGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGAt                                          >  1:1361609/1‑148 (MQ=255)
                                                                                                         gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATaa                                        <  1:341863/148‑1 (MQ=255)
                                                                                                            tGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACtt                                     <  1:839586/148‑1 (MQ=255)
                                                                                                             gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                    <  2:1010833/148‑1 (MQ=255)
                                                                                                              cAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTcc                                   >  1:3028989/1‑148 (MQ=255)
                                                                                                                 cAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGt                                <  2:1203886/148‑1 (MQ=255)
                                                                                                                   gAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCa                              <  2:3175768/148‑1 (MQ=255)
                                                                                                                     aTAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAg                            >  1:924110/1‑148 (MQ=255)
                                                                                                                     aTAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAg                            >  1:2977153/1‑148 (MQ=255)
                                                                                                                      tAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGc                           <  1:543220/148‑1 (MQ=255)
                                                                                                                       aGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTACAGTCAAGcc                          >  2:746067/1‑148 (MQ=255)
                                                                                                                        gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCa                         <  1:1369622/148‑1 (MQ=255)
                                                                                                                        gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCa                         <  1:1820157/148‑1 (MQ=255)
                                                                                                                        gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTACAGTCAAGCCa                         >  2:1833063/1‑148 (MQ=255)
                                                                                                                         tAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCaa                        <  1:2396294/148‑1 (MQ=255)
                                                                                                                          aaTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAAt                       <  2:2967766/148‑1 (MQ=255)
                                                                                                                           aTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATg                      >  1:2653553/1‑148 (MQ=255)
                                                                                                                            ttttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGt                     <  2:2894216/148‑1 (MQ=255)
                                                                                                                            ttttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGt                     <  2:2620337/148‑1 (MQ=255)
                                                                                                                            ttttGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGt                     <  1:1457053/148‑1 (MQ=255)
                                                                                                                             tttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTa                    <  1:186771/148‑1 (MQ=255)
                                                                                                                              ttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAg                   <  1:203784/148‑1 (MQ=255)
                                                                                                                                ggTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAggg                 >  2:518425/1‑148 (MQ=255)
                                                                                                                                ggTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAggg                 <  1:2686477/148‑1 (MQ=255)
                                                                                                                                ggTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAggg                 <  1:2397984/148‑1 (MQ=255)
                                                                                                                                  tCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCa               <  2:1981135/148‑1 (MQ=255)
                                                                                                                                  tCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCa               <  2:1709398/148‑1 (MQ=255)
                                                                                                                                   cTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATTTTCCACCGAGTTCAAAATCTTCATAAAGATAAATTACAGTCAAGCCAATTTAGGGCAt              >  1:678342/1‑148 (MQ=255)
                                                                                                                                   cTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAt              >  2:1002215/1‑148 (MQ=255)
                                                                                                                                   cTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAAAGATAAAATCCAGTCAAGCCCAAGTAGGGCAt              >  1:8673/1‑148 (MQ=255)
                                                                                                                                    tttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAATCCAATGTAGGGCAtt             >  1:2119243/1‑148 (MQ=255)
                                                                                                                                        aCCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGACAATGTACGGCATTTTaa         >  2:2399679/1‑148 (MQ=255)
                                                                                                                                          cTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTaaaa       >  2:1406359/1‑148 (MQ=255)
                                                                                                                                            tACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg     >  1:2456616/1‑148 (MQ=255)
                                                                                                                                            tACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg     <  2:895049/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTG  >  CP009273/580220‑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: