Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,295,828 T→C 5.3% I165V (ATT→GTT)  ccmC ← heme export ABC transporter permease; CcmE‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,295,8280TC 5.3% 235.0 / 2.5 75I165V (ATT→GTT) ccmCheme export ABC transporter permease; CcmE‑interacting protein
Reads supporting (aligned to +/- strand):  ref base T (38/33);  new base C (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

AGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTT  >  NC_000913/2295682‑2295976
                                                                                                                                                  |                                                                                                                                                    
agccAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCTGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCa                                                                                                                                                    <  8:262763/146‑1 (MQ=255)
aGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGTGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAAt                                                                                                                                                      >  4:185999/1‑147 (MQ=255)
      gAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCa                                                                                                                                              <  6:212929/149‑1 (MQ=255)
        gCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAgg                                                                                                                                            <  7:128470/149‑1 (MQ=255)
            gAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACATTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATAc                                                                                                                                        <  5:296549/149‑1 (MQ=255)
             aaGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATAcc                                                                                                                                       >  7:214568/1‑149 (MQ=255)
                ccAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGc                                                                                                                                    >  7:253189/1‑149 (MQ=255)
                  aaaaaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCg                                                                                                                                  >  8:266824/1‑149 (MQ=255)
                  aaaaaTCGACCAGCGCAGCGGCGAACGCATCGCCCGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCg                                                                                                                                  >  6:319253/1‑149 (MQ=255)
                     aaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAACCAGCACCAGGATACCTGCCGCa                                                                                                                                <  8:284354/148‑1 (MQ=255)
                     aaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAACCAGCACCAGGATACCTGCCGCa                                                                                                                                >  1:228631/1‑148 (MQ=255)
                     aaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAACCAGCACCAGGATACCTGCCGCa                                                                                                                                <  2:228631/148‑1 (MQ=255)
                     aaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAACCAGCACCAGGATACCTGCCGCa                                                                                                                                >  7:284354/1‑148 (MQ=255)
                      aTCGACCAGCGCAGCGGCGAACGCATCGCCGGATTGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACg                                                                                                                              >  3:266028/1‑149 (MQ=255)
                      aTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGAGCCCTGCTGCAGGGTGTTCCACCACGCCCCGGAGTAATGAATAATCGGCAGATTCCCCACGCGAGTGAGCAGCAGGATCCCTGCCggccg                                                                                                                              >  2:312282/1‑145 (MQ=255)
                               cgcagcgGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGc                                                                                                                        >  3:186626/1‑146 (MQ=255)
                                   gcggcgAACGCATCGCCGGATTGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGAcgg                                                                                                                 <  4:266028/149‑1 (MQ=255)
                                            gCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGgtcgtc                                                                                                         >  8:222939/1‑148 (MQ=255)
                                              aTCGCAGGATCGATACTTTGCTGCATCAGAGTTGATCACTGATGCAGTGTGTTCCACCACTACACGGAGTAATGAATAATCGGCAGATTCACCACTCCAATCAGCACCAGGATACCTGCCGCACTGCCCGCCAGACGGCGGtcgtcg                                                                                                        >  8:320255/1‑147 (MQ=255)
                                               tCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAg                                                                                                     <  6:23545/149‑1 (MQ=255)
                                                  ccGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCg                                                                                                  >  1:127982/1‑149 (MQ=255)
                                                  ccGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCg                                                                                                  <  1:65835/149‑1 (MQ=255)
                                                   cGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGAcggcg                                                                                                               >  4:312155/1‑135 (MQ=255)
                                                   cGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGAcggcg                                                                                                               <  3:312155/135‑1 (MQ=255)
                                                   cGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGt                                                                                                 <  7:222939/149‑1 (MQ=255)
                                                   cGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGt                                                                                                 >  7:235085/1‑149 (MQ=255)
                                                       tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCca                                                                                             >  4:311030/1‑149 (MQ=255)
                                                         gATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCca                                                                                             >  2:132027/1‑147 (MQ=255)
                                                             cTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAggg                                                                                        <  1:132027/148‑1 (MQ=255)
                                                             cTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAggg                                                                                        <  4:186626/148‑1 (MQ=255)
                                                             cTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGc                                                                                       <  7:266824/149‑1 (MQ=255)
                                                                 gctgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGTCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAAt                                                                                    <  5:126642/148‑1 (MQ=255)
                                                                 gcggcATCCGCGGTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAAt                                                                                    <  1:312282/145‑1 (MQ=255)
                                                                  ctgcATCCGCGTTGATCCCTGATGCAGGGTGTGCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATca                                                                                  <  6:147584/149‑1 (MQ=255)
                                                                              tGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCacac                                                                               <  4:44090/140‑1 (MQ=255)
                                                                              tGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCacac                                                                               >  3:44090/1‑140 (MQ=255)
                                                                                 tCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaa                                                                    >  6:40784/1‑148 (MQ=255)
                                                                                  cccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaaaa                                                                  >  4:29655/1‑149 (MQ=255)
                                                                                   ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaaaa                                                                  >  4:197088/1‑148 (MQ=255)
                                                                                   ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaaaa                                                                  >  1:121645/1‑148 (MQ=255)
                                                                                   ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaaaa                                                                  >  7:345287/1‑148 (MQ=255)
                                                                                   ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACa                                                                      <  6:325751/144‑1 (MQ=255)
                                                                                   ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACa                                                                      >  5:325751/1‑144 (MQ=255)
                                                                                    cTGATGCAGTGTGTTCCACCAATCCACTGATTAATTAATAATCGGCAGATTCACCACGCCAATCAGAACAAGGATACATGCCTCACGGCCCGCCAGACGGCGGTCGTCGAATGCGTGCCACATGGCAATCACACCCACATACAAAAAca                                                                <  7:320255/149‑1 (MQ=255)
                                                                                     tGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaaacccg                                                               >  4:352762/1‑145 (MQ=255)
                                                                                       aTGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcag                                                               <  8:235085/147‑1 (MQ=255)
                                                                                       aTGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcag                                                               <  5:319253/147‑1 (MQ=255)
                                                                                       aTGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcag                                                               <  2:127982/147‑1 (MQ=255)
                                                                                           aGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAc                                                         >  6:268626/1‑149 (MQ=255)
                                                                                            gggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAcc                                                        <  3:197088/149‑1 (MQ=255)
                                                                                               tgtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAg                                                                                                     >  1:4700/1‑101 (MQ=255)
                                                                                               tgtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAg                                                                                                     <  2:4700/101‑1 (MQ=255)
                                                                                                 ttCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGtt                                                    >  7:77997/1‑148 (MQ=255)
                                                                                                       cacTCCACGGAGTAATGAATAATCGGCAGATTCACTACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACAccc                                                                             <  4:21882/117‑1 (MQ=255)
                                                                                                       cacTCCACGGAGTAATGAATAATCGGCAGATTCACGACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACAccc                                                                             >  3:21882/1‑117 (MQ=255)
                                                                                                       cacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCcac                                                                           <  6:302541/119‑1 (MQ=255)
                                                                                                          tCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCa                                                  <  8:256496/141‑1 (MQ=255)
                                                                                                          tCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCa                                                  >  7:256496/1‑141 (MQ=255)
                                                                                                            cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGAc                                        <  8:253189/149‑1 (MQ=255)
                                                                                                            cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAACAGCAGCACCAGTTCAGAAGTCAGACg                                       >  7:135839/1‑149 (MQ=255)
                                                                                                            cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTAGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGAc                                        >  7:163309/1‑149 (MQ=255)
                                                                                                              cGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGt                                      <  8:163309/149‑1 (MQ=255)
                                                                                                                gAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAGAGCAGCACCAGTTCAGAAGTCAGACGTGc                                    <  5:40784/149‑1 (MQ=255)
                                                                                                                     aTGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc                                >  5:304925/1‑148 (MQ=255)
                                                                                                                       gAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCa                              >  7:211641/1‑148 (MQ=255)
                                                                                                                        aataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATa                            <  3:185999/149‑1 (MQ=255)
                                                                                                                         ataatCGGCAGATTCACCACGCCAATCAGCCCCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAc                           >  2:311426/1‑149 (MQ=255)
                                                                                                                          taatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAcc                          >  7:371042/1‑149 (MQ=255)
                                                                                                                              cGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACccac                       >  2:205123/1‑148 (MQ=255)
                                                                                                                                gCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCCTACCCACCAg                    >  6:33891/1‑149 (MQ=255)
                                                                                                                                     ttCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGcc               <  3:29655/149‑1 (MQ=255)
                                                                                                                                          cacGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCcaca          >  3:242251/1‑149 (MQ=255)
                                                                                                                                             gCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCg       >  1:202621/1‑149 (MQ=255)
                                                                                                                                               cAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACagcagc                                                           >  8:197094/1‑95 (MQ=255)
                                                                                                                                               cAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACagcagc                                                           <  7:197094/95‑1 (MQ=255)
                                                                                                                                                aaTCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGtt    <  8:371042/149‑1 (MQ=255)
                                                                                                                                                  tCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGtttt  <  8:106496/149‑1 (MQ=255)
                                                                                                                                                  tCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGtttt  >  1:179203/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTT  >  NC_000913/2295682‑2295976

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