Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,106,685 G→T 100% S51Y (TCT→TAT)  yghG ← secretin (GspDbeta) OM localization lipoprotein pilotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,106,6850GT97.8% 882.4 / 7.7 277S51Y (TCT→TAT) yghGsecretin (GspDbeta) OM localization lipoprotein pilotin
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base T (125/146);  total (131/146)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 279 total aligned reads displayed.

TCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTA  >  CP009273/3106539‑3106829
                                                                                                                                                  |                                                                                                                                                
tCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCatat                                                                                                                                                <  2:1478053/149‑1 (MQ=255)
   ccGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTa                                                                                                                                             >  2:2740233/1‑149 (MQ=255)
   ccGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTa                                                                                                                                             <  1:2707572/149‑1 (MQ=255)
   ccGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTa                                                                                                                                             <  2:1820154/149‑1 (MQ=255)
   ccGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTa                                                                                                                                             <  1:66846/149‑1 (MQ=255)
    cGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTaa                                                                                                                                            <  1:3770342/149‑1 (MQ=255)
     gTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAAt                                                                                                                                           <  1:569315/149‑1 (MQ=255)
     gTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAAt                                                                                                                                           >  2:3523567/1‑149 (MQ=255)
     gTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGGACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAAt                                                                                                                                           <  1:569282/149‑1 (MQ=255)
        gggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCgg                                                                                                                                        <  2:749467/149‑1 (MQ=255)
         ggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGc                                                                                                                                       >  1:2673415/1‑149 (MQ=255)
          gTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCa                                                                                                                                      <  2:3435358/149‑1 (MQ=255)
           tCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAg                                                                                                                                     >  2:1252888/1‑149 (MQ=255)
            cAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAgg                                                                                                                                    >  1:527760/1‑149 (MQ=255)
              gcgcACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGtt                                                                                                                                  >  1:995953/1‑149 (MQ=255)
                 cacaTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTg                                                                                                                               <  1:1244484/149‑1 (MQ=255)
                    aTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCt                                                                                                                            >  2:1650949/1‑149 (MQ=255)
                    aTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCt                                                                                                                            <  2:619650/149‑1 (MQ=255)
                    aTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCt                                                                                                                            <  2:2832993/149‑1 (MQ=255)
                      ttGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGa                                                                                                                          <  2:1415637/149‑1 (MQ=255)
                      ttGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGa                                                                                                                          <  2:1119316/149‑1 (MQ=255)
                            ttGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAtttt                                                                                                                    >  1:1004259/1‑149 (MQ=255)
                             tGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTg                                                                                                                   >  1:3406436/1‑149 (MQ=255)
                                 aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                               <  1:1358761/149‑1 (MQ=255)
                                 aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                               >  2:2149293/1‑149 (MQ=255)
                                  gCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCt                                                                                                              >  2:2341255/1‑149 (MQ=255)
                                    tGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGt                                                                                                            >  1:1138948/1‑149 (MQ=255)
                                    tGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGt                                                                                                            >  2:1186083/1‑149 (MQ=255)
                                     ggTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtt                                                                                                           >  2:245652/1‑149 (MQ=255)
                                     ggTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtt                                                                                                           >  1:3243364/1‑149 (MQ=255)
                                      gTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGttt                                                                                                          <  2:3780632/149‑1 (MQ=255)
                                        tAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGttttt                                                                                                        <  2:3383461/149‑1 (MQ=255)
                                         aaaaaGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtttttt                                                                                                       <  2:435306/149‑1 (MQ=255)
                                           aaaGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCg                                                                                                     >  2:2214337/1‑149 (MQ=255)
                                              ggcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCa                                                                                                  <  2:1433708/149‑1 (MQ=255)
                                                  tcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTaa                                                                                              >  1:2430147/1‑149 (MQ=255)
                                                    aggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTTCGCCAGCACTAAGGTATATACCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAAc                                                                                            >  2:227671/1‑149 (MQ=255)
                                                    aggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAAc                                                                                            <  2:546944/149‑1 (MQ=255)
                                                    aggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAAc                                                                                            <  1:1291008/149‑1 (MQ=255)
                                                      gcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                          <  2:490084/149‑1 (MQ=255)
                                                       cgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTgg                                                                                         <  1:870982/149‑1 (MQ=255)
                                                        gtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGc                                                                                        <  1:294592/149‑1 (MQ=255)
                                                        gtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGc                                                                                        >  2:3645075/1‑149 (MQ=255)
                                                         tcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCa                                                                                       <  2:1498671/149‑1 (MQ=255)
                                                         tcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCa                                                                                       >  1:2243161/1‑149 (MQ=255)
                                                          cTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAt                                                                                      >  1:1828494/1‑149 (MQ=255)
                                                           tGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAtt                                                                                     >  2:3906791/1‑149 (MQ=255)
                                                           tGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAtt                                                                                     >  2:1578919/1‑149 (MQ=255)
                                                            gCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAttt                                                                                    >  2:864661/1‑149 (MQ=255)
                                                            gCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAttt                                                                                    <  1:327214/149‑1 (MQ=255)
                                                             cGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAtttt                                                                                   <  1:2270462/149‑1 (MQ=255)
                                                                  tttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTa                                                                              <  2:2011114/149‑1 (MQ=255)
                                                                  tttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTa                                                                              >  1:3209230/1‑149 (MQ=255)
                                                                         aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctgg                                                                       >  1:3049299/1‑149 (MQ=255)
                                                                         aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctgg                                                                       >  2:1696022/1‑149 (MQ=255)
                                                                         aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctgg                                                                       >  1:3599889/1‑149 (MQ=255)
                                                                             gCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGcaca                                                                   >  2:3677752/1‑149 (MQ=255)
                                                                              cTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAg                                                                  <  2:3823585/149‑1 (MQ=255)
                                                                                 cGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCa                                                               >  2:3306851/1‑149 (MQ=255)
                                                                                   cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTACGGAATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTGGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACt                                                             <  2:637162/149‑1 (MQ=255)
                                                                                      aTAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaa                                                          >  2:3268656/1‑149 (MQ=255)
                                                                                       tAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaat                                                         >  1:1229164/1‑149 (MQ=255)
                                                                                       tAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGACATTTTCATTATGGCTGGCACAGCCACTtaat                                                         >  1:146434/1‑149 (MQ=255)
                                                                                         aTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaataa                                                       >  2:368692/1‑149 (MQ=255)
                                                                                         aTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaataa                                                       <  1:1577967/149‑1 (MQ=255)
                                                                                          tGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTaataat                                                      >  1:424718/1‑149 (MQ=255)
                                                                                           ggTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATc                                                     >  1:2659373/1‑149 (MQ=255)
                                                                                            gTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                    >  1:933591/1‑149 (MQ=255)
                                                                                            gTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                    <  1:4095427/149‑1 (MQ=255)
                                                                                             tCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCt                                                   <  1:3995289/149‑1 (MQ=255)
                                                                                             tCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCt                                                   <  2:2847203/149‑1 (MQ=255)
                                                                                                tttttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTa                                                <  1:3237172/149‑1 (MQ=255)
                                                                                                 ttttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTaa                                               <  1:1941041/149‑1 (MQ=255)
                                                                                                 ttttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTaa                                               <  2:695963/149‑1 (MQ=255)
                                                                                                 ttttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTaa                                               <  2:1302771/149‑1 (MQ=255)
                                                                                                 ttttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTaa                                               <  2:1302737/149‑1 (MQ=255)
                                                                                                       cGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTca                                         >  1:803144/1‑149 (MQ=255)
                                                                                                        gTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaa                                        >  1:1153651/1‑149 (MQ=255)
                                                                                                              ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAg                                  <  1:1040890/149‑1 (MQ=255)
                                                                                                                gCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCg                                >  2:1143263/1‑149 (MQ=255)
                                                                                                                gCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCg                                <  2:241646/149‑1 (MQ=255)
                                                                                                                gCTACTTTGCGCCAGCACTAAGGTATATCCCACATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCg                                >  2:1770889/1‑149 (MQ=255)
                                                                                                                   aCTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGata                             <  1:3167623/149‑1 (MQ=255)
                                                                                                                    cTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGatat                            >  2:409821/1‑149 (MQ=255)
                                                                                                                    cTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGatat                            >  2:409732/1‑149 (MQ=255)
                                                                                                                         cgcCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATa                       <  1:3884216/149‑1 (MQ=255)
                                                                                                                           ccAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAcc                     >  1:238632/1‑149 (MQ=255)
                                                                                                                           ccAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAcc                     <  2:1957634/149‑1 (MQ=255)
                                                                                                                            cAGCACTAAGGTATATCCCTCATATTTAATCGGCAGGTTTTTGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCAACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccc                    <  2:191188/149‑1 (MQ=255)
                                                                                                                              gCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctc                  >  1:2239009/1‑149 (MQ=255)
                                                                                                                                aCTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccc                >  1:961464/1‑149 (MQ=255)
                                                                                                                                  tAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg              >  1:2022113/1‑149 (MQ=255)
                                                                                                                                       tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt         <  2:2673415/149‑1 (MQ=255)
                                                                                                                                       tataTCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt         >  1:3012934/1‑149 (MQ=255)
                                                                                                                                        atatCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAttt        <  1:1439194/149‑1 (MQ=255)
                                                                                                                                          atCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGt      <  1:2186573/149‑1 (MQ=255)
                                                                                                                                           tCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt     <  2:3142324/149‑1 (MQ=255)
                                                                                                                                            cccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGttt    >  1:2160937/1‑149 (MQ=255)
                                                                                                                                            cccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAAAAATATCGATATTAATACCCTCCCTGGCATTTGttt    <  1:3425906/149‑1 (MQ=255)
                                                                                                                                             ccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt   <  1:2786334/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTA  >  CP009273/3106539‑3106829

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: