Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,456,014 C→A 5.5% H147N (CAT→AAT)  gspC → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,456,0140CA 5.5% 195.5 / 3.6 74H147N (CAT→AAT) gspCgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base A (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

CTCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGAACCCTCGTA  >  NC_000913/3455867‑3456161
                                                                                                                                                   |                                                                                                                                                   
cTCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGCGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGCGTTTGAAACTCACGGGcc                                                                                                                                                    >  2:69704/1‑148 (MQ=255)
 tCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAATAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACAAGCGTTGTGTTTGAAACTCACGGGCAt                                                                                                                                                   >  1:657248/1‑149 (MQ=255)
 tCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCAt                                                                                                                                                   >  2:920685/1‑149 (MQ=255)
 tCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCAt                                                                                                                                                   >  2:359882/1‑149 (MQ=255)
 tCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCAt                                                                                                                                                   >  1:741198/1‑149 (MQ=255)
 tCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCAt                                                                                                                                                   >  1:707212/1‑149 (MQ=255)
 tCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCAt                                                                                                                                                   >  1:397476/1‑149 (MQ=255)
  cacaTCGAATGATGAGACTTCTTTTGTCTTAAATAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCAtt                                                                                                                                                  <  1:98829/149‑1 (MQ=255)
   acaTCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGACACTCACGTGCATTa                                                                                                                                                 >  1:344870/1‑149 (MQ=255)
      tCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGa                                                                                                                                              <  1:247163/149‑1 (MQ=255)
          atgatgAGTCTTCTTTTGTCTGCATTTATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGaaaaa                                                                                                                                          <  2:861598/149‑1 (MQ=255)
              tgaGACTTCTTTTTTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAGATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTAAGAAAAAGTaa                                                                                                                                      <  2:441856/149‑1 (MQ=255)
              tgaGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTaa                                                                                                                                      >  2:578935/1‑149 (MQ=255)
                 gaCTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAActc                                                                                                                                   >  2:222630/1‑149 (MQ=255)
                    ttcttTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCtcc                                                                                                                                >  2:270607/1‑149 (MQ=255)
                          cggcTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCTGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAg                                                                                                                          <  2:857840/146‑1 (MQ=255)
                                       gAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGATACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGAcat                                                                                                             <  1:790593/149‑1 (MQ=255)
                                       gAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTGGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCA‑‑GATGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGAcatca                                                                                                           <  2:1065952/149‑1 (MQ=255)
                                             aGCGAACAACAGCGGTGCTCCCTCAACGAAGCCCGTGGGTCAGCGCCGGGAACGTTTAGTAGAAAAATAAATAAACCCAGCGTTGTGTTTGAAACTCACAGGCATTACGAAAAAGTATCTCTCCATCCAGGGCTACCTGACATCATCaa                                                                                                       <  2:628873/149‑1 (MQ=255)
                                             aGCGAACAAAAGCTATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGTGTTGTGTTTGAAGCTCACGTGCGTTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCaa                                                                                                       <  2:192598/149‑1 (MQ=255)
                                                 aaCAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCCTCAagcag                                                                                                   >  1:21954/1‑149 (MQ=255)
                                                 aaCAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAagcag                                                                                                   >  2:646524/1‑149 (MQ=255)
                                                  aCAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCTATCCAGGGCTACCTGACATCATCAagcagc                                                                                                  <  1:600411/149‑1 (MQ=255)
                                                    aaaaGCGAAACTCCCTCAACGAAGCCCTTGAGTGAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACACCGTCAAGCAGCCt                                                                                                >  1:841723/1‑149 (MQ=255)
                                                      acctgATCCTCCAGCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACCCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGa                                                                                              <  2:805074/145‑1 (MQ=255)
                                                      aaGCGATACTCCCTCAACGAAACCCTTGAGTCAGAGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTGGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGa                                                                                              >  2:872974/1‑149 (MQ=255)
                                                           attcTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAACCTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAg                                                                                         <  2:516299/146‑1 (MQ=255)
                                                              cTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGaaa                                                                                      <  2:1011081/149‑1 (MQ=255)
                                                                 ccTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCCTCGAGCAGCCGTATTCAGAAAGTc                                                                                   >  2:467356/1‑149 (MQ=255)
                                                                 ccTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGAATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTc                                                                                   >  2:196868/1‑149 (MQ=255)
                                                                 ccTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGAATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTc                                                                                   >  2:196846/1‑149 (MQ=255)
                                                                  cTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCa                                                                                  <  2:397476/149‑1 (MQ=255)
                                                                      aCGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATACAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAAt                                                                              <  1:670447/149‑1 (MQ=255)
                                                                        gAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGAAGCCTGATTCAGAAAGACAAAATGt                                                                            >  2:602614/1‑149 (MQ=255)
                                                                            gacgtGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGAACTc                                                                        <  2:707212/145‑1 (MQ=255)
                                                                             ccTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCg                                                                       >  1:507907/1‑149 (MQ=255)
                                                                              cTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCTAGGTCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGc                                                                      >  2:323204/1‑149 (MQ=255)
                                                                               tggAGTCAGCGCCGGGAACGTTTATCAGAAAAGCAGAGGAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCt                                                                     <  1:601488/147‑1 (MQ=255)
                                                                                tGAGTCAGCGCCGGGAACGTTTAGTAGAAGAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTg                                                                    <  1:80944/149‑1 (MQ=255)
                                                                                    ttaGCGGCGGGAACGTTTATTAGAAAAGTGAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCGTTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGatta                                                                >  2:596449/3‑149 (MQ=255)
                                                                                    tCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTCCCCGACATCATCAACCACCCTGATCCAGAACGTCAAAATCAACTCCCTGatta                                                                >  2:130819/1‑149 (MQ=255)
                                                                                    tCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGatta                                                                >  1:940102/1‑149 (MQ=255)
                                                                                          ccGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACCCACGACCGCGACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGattatattat                                                          <  1:1128807/149‑1 (MQ=255)
                                                                                           cGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCCCTGATGATATTATc                                                         >  2:943386/1‑149 (MQ=255)
                                                                                             ggAACGTTTATTCGAAAAATGAATGAGACCAGCGTTGTGTTTGAGACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACACCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCga                                                       <  2:770867/149‑2 (MQ=255)
                                                                                              gAACGTTTTATAGAAAAATAAATAAAGCCAGCGTTGTGTTTGAAACTCACGGGCAATACGAAAAAGTAACTCTCCATCCAGGGCTCCCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGcc                                                      <  2:581481/149‑1 (MQ=255)
                                                                                                aCGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCAc                                                    <  1:683132/149‑1 (MQ=255)
                                                                                                  gTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGc                                                  <  1:578935/149‑1 (MQ=255)
                                                                                                        taGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATAc                                            <  1:1074830/149‑1 (MQ=255)
                                                                                                         aGAAAAATAAATCACACCAGCTTTGTGTTTGAAGCTCACGGGCATTACGAAAAAGTAACGCTCCATCCAGGGCTACCGGACATCATCGAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACg                                           <  2:885688/149‑1 (MQ=255)
                                                                                                          gAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCACTGATTATCTTATCGCCACGCCTATACGt                                          >  2:916351/1‑149 (MQ=255)
                                                                                                              aataaataAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATg                                      <  1:359882/149‑1 (MQ=255)
                                                                                                              aataaacaAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATg                                      <  1:935783/149‑1 (MQ=255)
                                                                                                              aataaaaaAAGCCAGCGTCGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATg                                      <  2:433344/149‑1 (MQ=255)
                                                                                                                taaataaaACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGga                                    <  2:663100/149‑1 (MQ=255)
                                                                                                                      aaaCCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGTAGAACAg                              >  2:1068305/1‑149 (MQ=255)
                                                                                                                      aaaCCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAg                              >  1:471659/1‑149 (MQ=255)
                                                                                                                       aaCCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGTCATTTTCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGa                             <  1:270607/149‑1 (MQ=255)
                                                                                                                             cGTTGTGTTTGAAACCCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACg                       <  1:855943/149‑1 (MQ=255)
                                                                                                                              ggtgtgTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCGAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACgg                      <  2:84118/147‑1 (MQ=255)
                                                                                                                              gTTGTGTTTGAAACTCACGGGAATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACgg                      <  1:196868/149‑1 (MQ=255)
                                                                                                                              gTTGTGTTTGAAACTCACGGGAATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACgg                      <  1:196846/149‑1 (MQ=255)
                                                                                                                                 gtgtTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCAGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACACACTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCt                   >  2:635964/1‑149 (MQ=255)
                                                                                                                                    tttGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGAGATGGAGAACAGATATACGGGCTTCg                >  2:8215/1‑149 (MQ=255)
                                                                                                                                    tttGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATAATATTATCGCCACGCCTATACGAGATGGAGAACAGATATACGGGCTTCg                >  1:161789/1‑149 (MQ=255)
                                                                                                                                    tgtGAAACTCACGGGCATTACCAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACTCCTATACGTGATGGAGAACAGATATACGGGCTTCg                <  2:695392/147‑1 (MQ=255)
                                                                                                                                     ttGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGa               <  1:646524/149‑1 (MQ=255)
                                                                                                                                      tGAAACTCACGGGCATTACGAAAAAGTAATTCTCCATCCAGGGCTATCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGAt              <  1:368779/149‑1 (MQ=255)
                                                                                                                                       gAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAAAAGAAATACGGGCTTCGAtt             >  2:205790/1‑149 (MQ=255)
                                                                                                                                          aCTCATGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCCAGTAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGaa          <  2:344870/149‑1 (MQ=255)
                                                                                                                                           cTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCCCGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGAAc         >  2:895017/1‑149 (MQ=255)
                                                                                                                                            tCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGAAcc        >  2:1076634/1‑149 (MQ=255)
                                                                                                                                               ccggCCTTACGAAAAAGTAACCCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGAACACTc     <  2:877837/147‑1 (MQ=255)
                                                                                                                                               cGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGAACCCTc     <  1:872974/149‑1 (MQ=255)
                                                                                                                                                  gCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGAACCCTCGTa  >  1:71671/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTCACATCGAATGATGAGACTTCTTTTGTCTTAATTAATGAGGGGAGCGAACAAAAGCGATACTCCCTCAACGAAGCCCTTGAGTCAGCGCCGGGAACGTTTATTAGAAAAATAAATAAAACCAGCGTTGTGTTTGAAACTCACGGGCATTACGAAAAAGTAACTCTCCATCCAGGGCTACCTGACATCATCAAGCAGCCTGATTCAGAAAGTCAAAATGTACTCGCTGATTATATTATCGCCACGCCTATACGTGATGGAGAACAGATATACGGGCTTCGATTGAACCCTCGTA  >  NC_000913/3455867‑3456161

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